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1.0  INTRODUCTION 

This  report  describes  the  second  year  of  work  on  contract  number  N00014- 
79-F-0005.  The  work  has  the  following  general  goals: 

Demonstration  of  a  superconducting  300  megasample  per  second 
A/D  converter  within  the  first  three  years. 

Demonstration  of  a  superconducting  A/D  converter  operating 
well  Into  the  glgasample  per  second  region  within  five  years. 

During  this  second  year  the  goal  for  the  third  year  was  essentially 
accomplished  and  described  In  a  publication  In  Applied  Physics  Letters 
on  an  A/D  converter  operating  between  200  and  700  megasampl es  per  second. 
This  remarkable  success.  In  a  totally  new  technology,  was  the  result  of 
advances  In  three  areas: 

Early  In  the  reporting  period,  the  catalog  of  processes 
necessary  for  fabricating  superconducting  Integrated  circuits 
was  completed,  and  the  first  circuits  were  successfully 
demonstrated.  This  work  had  been  the  primary  subject  of 
the  first  year’s  contract. 

Simulations  of  the  threshold  curves  of  two-  and  three- 
junction  Interferometers  enabled  understanding  of 
experimental  measurements  of  the  first  circuits. 

A  simple  A/D  converter  design  was  Invented.  The  design 
greatly  reduced  the  complexity  of  the  expected  chip  layout 
and  testing,  making  possible  an  early  demonstration  of  a 
4-blt  A/D  converter. 

In  Section  2  we  summarize  the  main  portion  of  the  work.  Details  are 
presented  as  reprints  In  the  Appendix.  The  work  described  Is  the  demon¬ 
stration  of  a  very  high  speed  superconducting  A/D  converter,  the  calculation 
of  threshold  curves  of  Interferometers,  theoretical  analysis  of  an 
alternate  design  of  an  A/0  converter,  and  experimental  and  theoretical 
work  on  strlpllnes. 

In  Section  3  we  describe  other  work,  supported  by  0NR,  but  under  a 
different  contract.  Much  of  this  work  Is  closely  related  to  the  A/D 
converter  work.  He  discuss  the  demonstration  of  a  9  ps  risetime  sampler 
and  a  demonstration  of  the  quietest  known  microwave  mixer  at  36  GHz. 


Section  4  contains  a  discussion  of  those  feu  elements  of  the  project 
which  have  not  been  submitted  for  publication.  These  are  a  complete 
description  of  the  present  fabrication  process,  the  details  of  the 
process  for  fabrication  of  resistors,  the  fabrication  of  several  test 
circuits,  and  a  few  comments  on  factors  limiting  yield. 


2.0  CONTRACT  WORK  PUBLISHED  IN  THE  LITERATURE 

The  work  which  directly  resulted  from  the  ONR  contract  on  high  speed  A/0 
conversion  will  result  In  five  publications.  These  papers  are  Included  In 
the  Appendix  and  are  described  briefly  below: 

2.1  Multiple-Quantum  Interference  Superconducting  A/0  Converter 
R.  E.  Harris,  C.  A.  Hamilton,  and  F.  L.  Lloyd 

Appl .  Phys.  Letters,  July  15,  1979. 

This  paper  describes  the  essential  work  under  the  contract.  A  novel 
design  was  Invented  for  a  very  simple,  fully  parallel  A/D  converter. 

The  design  makes  possible  an  N-blt  converter  using  only  N  gates.  The 
Gray  code  output  of  the  converter  results  naturally  from  the  design. 
Sampling  rates  for  the  two  bits  tested  were  2  x  10°  per  second.  Con¬ 
ventional  sampling  measurements  on  one  bit  suggested  the  entire  device 
will  operate  at  7  x  10s  per  second.  The  major  limitation  appears  to  be 
cross-talk  between  the  leads  to  the  chip.  It  1$  expected  thatgfuture 
versions  of  the  devices  will  operate  at  speeds  well  Into  the  10  samples 
per  second  region. 

2.2  Analysis  of  Threshold  Curves  for  Superconducting  Interferometers 
R.  L.  Peterson  and  C.  A.  Hamilton 

J.  Appl.  Physics,  50,  8135  (1979). 

Theoretical  work  In  support  of  the  experimental  demonstration  of  a  very 
high  speed  A/D  converter  was  concentrated  on  the  analysis  of  threshold 
curves  of  superconducting  Interferometers.  These  devices  are  the  gates 
used  In  the  A/D  converter.  The  details  of  their  Interference  patterns 
provide  the  basis  for  the  operation  of  the  A/D  converter. 

In  the  absence  of  perfect  fabrication  and  testing  In  a  zero-field,  zero- 
transient  environment,  superconducting  three-junction  Interferometers  do 
not  behave  In  the  Ideal  manner  described  previously  In  the  literature. 

This  theoretical  work  provides  a  description  of  non-ideal  behavior  due 
to  asymmetries  In  junction  critical  currents  and  loop  Inductances,  and 
due  to  trapped  flux.  The  work  provides  the  basis  for  understanding  how 
nearly  Identical  critical  currents  and  loop  Inductances  must  be  In  order 
to  achieve  Interferometers  having  the  required  characteristics.  It 
further  provides  motivation  for  using  two- junction  Interferometers  In 
the  future,  rather  than  three- junction  ones  as  In  the  past. 
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2.3  Analog-to-DIgltal  Conversion  with  a  SQUID:  Conditions  for  a 
Countable  Pulse  Train 

R.  L.  Peterson 

J.  Appl.  Phys.  50,  4231  (1979). 

This  work  provides  a  theoretical  analysis  of  a  different  type  of  super¬ 
conducting  A/0  converter  than  the  one  demonstrated  under  the  contract. 

This  one  Is  based  on  a  Josephson  junction  In  a  superconducting  loop  (single 
junction  SQUID).  The  device  generates  pulses  as  the  magnetic  flux  through 
the  loop  Is  raised  and  lowered  by  an  applied  field.  Counting  the  pulses 
provides  a  measure  of  the  analog  value  of  the  magnetic  flux.  The  work 
determines  conditions  on  circuit  parameters  required  for  well-resolved 
pulses.  The  capacitance  Is  shown  to  be  especially  critical. 

2.4  Attenuation  In  Superconducting  Strlpllnes 
R.  L.  Kautz 

IEEE  Trans.  Mag.  MAG-15.  566  (1979). 

Superconducting  strlpllnes  provide. an  Important  part  of  the  advantage  of 
superconducting  electronic  technology.  They  are  essentially  lossless 
and  dlsperslonless,  providing  the  capability  of  transmitting  signals  with 
frequency  components  In  the  THz  (1 0,c)  region  over  relatively  long  distances. 
This  work  provides  an  experimental  measure  of  the  losses  In  such  lines.  At 
4  K  the  attenuation  between  50  and  500  MHz  Is  proportional  to  frequency  and 
probably  results  from  dielectric  losses.  Extrapolation  of  the  dielectric 
losses  to  10  GHz  yields  an  attenuation  of  10  dB/m,  leading  to  losses  accept¬ 
able  for  most  microcircuit  applications  where  strlpllne  lengths  are  1  cm  or 
less. 

2.5  Miniaturization  of  Normal -state  and  Superconducting  Strlpllnes 
R.  L.  Kautz 

J.  Res.  Nat'l .  Bur.  Stds.  84,  247  (1979). 

Theoretical  work  has  been  performed  to  quantitatively  determine  the  advantages 
offered  by  superconducting  strlpllnes,  as  compared  with  normal  ones.  The  work 
employs  all  of  the  relevant  theory  to  characterize  both  superconducting  and 
normal  strlpllnes.  For  normal  conductors,  the  Reuter- Sondhelmer  theory  Is 
applied  In  order  to  account  for  the  effects  of  finite  film  thickness  and  mean 
free  path.  For  superconductors  the  Mattls- Bardeen  theory  Is  used  In  order  to 
Include  effects  due  to  the  energy  gap.  Comparisons  are  made  for  three  example 
conductors:  copper  at  295  K  and  4.2  K  and  niobium  at  4.2  K.  Simulations  of 
pulse  transmission  on  lines  having  0.2  micron  dielectric  thickness  Indicate 
limitations  of  the  copper  lines  for  transmitting  50  ps  pulses  (half-width  at 
half-maximum)  more  than  0.1  cm  for  295  K  and  1.0  cm  for  4.2  K.  On  the  other 
hand,  superconducting  lines  will  carry  1  ps  pulses  10  cm  with  little  distortion. 
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3.0  RELATED  PUBLISHED  WORK  SUPPORTED  UNDER  OTHER  ONR  CONTRACTS 

Additional  work  has  been  supported  by  the  Office  of  Naval  Research  under 
contract  number  N00014-79-F-0020  and,  of  course,  by  the  National  Bureau  of 
Standards.  Much  of  this  work  Is  closely  related  to  the  ONR  A/D  converter 
work  for  two  reasons.  First  It  demonstrates  the  technical  capability  which 
the  Cryoelectronics  Group  can  bring  to  bear  on  the  A/D  converter  problem, 
and  second.  It  provides  tools  which  can  be  used  In  future  A/D  converter 
work.  Two  relevant  papers  are  discussed  below  and  also  Included  In  the 
Appendix. 

3.1  A  Superconducting  Sampler  for  Josephson  Logic  Circuits 

C.  A.  Hamilton,  F.  L.  Lloyd,  R.  L.  Peterson,  and  J.  R.  Andrews 
Appl.  Phys.  Letters,  Nov.  1,  1979 

A  superconducting  on-chip  sampler  has  been  developed  which  measures  rlsetlmes 
as  short  as  9  ps.  The  sampler  uses  a  single  biased  Josephson  junction  on  the 
chip  containing  the  signal  to  be  measured.  The  junction  converts  a  given 
signal  level  to  a  time  delay  which  Is  measured  outside  the  cryogenic  region 
by  a  conventional  sampling  oscilloscope.  Such  samplers  will  be  used  In  the 
future  to  fully  characterize  the  speed  of  A/D  converters. 

3.2  Quasiparticle  Heterodyne  Mixing  In  SIS  Tunnel  Junctions 
P.  L.  Richards,  T.  M.  Shen,  R.  E.  Harris,  and  F.  L.  Lloyd 
Appl.  Phys.  Letters,  34,  345  (1979). 

The  fabrication  technology  which  has  been  developed  was  used  to  fabricate 
microwave  mixers  which  were  evaluated  at  36  GHz  and  found  to  be  the  quietest 
mixers  yet  developed  for  that  frequency.  The  work  was  done  In  collaboration 
with  the  University  of  California  at  Berkeley  where  the  testing  was  done. 
Although  the  mixers  were  superconductor-lnsulator-superconductor  junctions 
which  exhibited  the  Josephson  effect.  It  was  the  energy  gap  structure  which 
provided  the  basis  for  the  mixing.  The  project  Is  continuing  In  the  direction 
of  multiple  junction  series  arrays  and  higher  frequencies. 


4.0  UNPUBLISHED  RESULTS 

This  Section  describes  some  unpublished  results  which  may  be  useful  to 
others  In  the  field.  All  of  the  work  was  either  on  a  procedure  necessary 
to  accomplish  other  goals,  or  was  on  a  device  built  and  tested  In  order 
to  learn  the  technology. 

7  4.1  Special  Equipment 

-#■ 

A  variety  of  special  pieces  of  equipment  are  now  In  use  and  are  described 
-  below.  Many  of  these  Items  were  specially  designed  and  constructed  for 

this  work  and  have  played  a  vital  role  In  the  accomplishments  to  date. 

■*  Of  particular  Interest,  several  Items  of  test  equipment  were  required  for 

the  characterization  of  our  superconducting  A/D  converters. 
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3.0  RELATED  PUBLISHED  WORK  SUPPORTED  UNDER  OTHER  ONR  CONTRACTS 


Additional  work  has  been  supported  by  the  Office  of  Naval  Research  under 
contract  number  N00014-79-F-0020  and,  of  course,  by  the  National  Bureau  of 
Standards.  Much  of  this  work  Is  closely  related  to  the  ONR  A/D  converter 
work  for  two  reasons.  First  It  demonstrates  the  technical  capability  which 
the  Cryoelectronics  Group  can  bring  to  bear  on  the  A/D  converter  problem, 
and  second.  It  provides  tools  which  can  be  used  In  future  A/D  converter 
work.  Two  relevant  papers  are  discussed  below  and  also  Included  In  the 
Appendix. 

3.1  A  Superconducting  Sampler  for  Josephson  Logic  Circuits 

C.  A.  Hamilton,  F.  L.  Lloyd,  R.  L.  Peterson,  and  J.  R.  Andrews 
Appl.  Phys.  Letters,  Nov.  1,  1979 

A  superconducting  on-chip  sampler  has  been  developed  which  measures  rlsetlmes 
as  short  as  9  ps.  The  sampler  uses  a  single  biased  Josephson  junction  on  the 
chip  containing  the  signal  to  be  measured.  The  junction  converts  a  given 
signal  level  to  a  time  delay  which  Is  measured  outside  the  cryogenic  region 
by  a  conventional  sampling  oscilloscope.  Such  samplers  will  be  used  in  the 
future  to  fully  characterize  the  speed  of  A/D  converters. 

3.2  Quasiparticle  Heterodyne  Mixing  In  SIS  Tunnel  Junctions 
P.  L.  Richards,  T.  M.  Shen,  R.  E.  Harris,  and  F.  L.  Lloyd 
Appl.  Phys.  Letters,  34,  345  0979). 

The  fabrication  technology  which  has  been  developed  was  used  to  fabricate 
microwave  mixers  which  were  evaluated  at  36  GHz  and  found  to  be  the  quietest 
mixers  yet  developed  for  that  frequency.  The  work  was  done  In  collaboration 
with  the  University  of  California  at  Berkeley  where  the  testing  was  done. 
Although  the  mixers  were  superconductor- Insulator-superconductor  junctions 
which  exhibited  the  Josephson  effect.  It  was  the  energy  gap  structure  which 
provided  the  basis  for  the  mixing.  The  project  Is  continuing  In  the  direction 
of  multiple  junction  series  arrays  and  higher  frequencies. 


4.0  UNPUBLISHED  RESULTS 

This  Section  describes  some  unpublished  results  which  may  be  useful  to 
others  In  the  field.  All  of  the  work  was  either  on  a  procedure  necessary 
to  accomplish  other  goals,  or  was  on  a  device  built  and  tested  In  order 
to  learn  the  technology. 

4.1  Special  Equipment 

A  variety  of  special  pieces  of  equipment  are  now  In  use  and  are  described 
below.  Many  of  these  Items  were  specially  designed  and  constructed  for 
this  work  and  have  played  a  vital  role  In  the  accomplishments  to  date. 

Of  particular  Interest,  several  Items  of  test  equipment  were  required  for 
the  characterization  of  our  superconducting  A/D  converters. 


4.1.1  High  Speed  Probe  -  The  most  crucial  Item  has  been  a  probe  for  mounting 
the  sample  Integrated  circuit  chips  and  delivering  to  and  from  them  high 
speed  signals  over  a  large  number  of  lines.  The  basic  design  of  this  device 
was  reported  last  year.  Since  then  five  more  probes  have  been  built  with 
numerous  changes  designed  to  Improve  the  three  most  Important  requirements: 
contact  reliability,  large  bandwidth,  and  absence  of  crosstalk  between  lines. 
Contact  reliability  has  been  Improved  by  developing  a  special  printed  circuit 
board  on  which  the  contact  fingers  are  formed  from  beryllium  copper  and  thus 
have  Individual  compliance  with  every  chip  pad.  Bandwidth  has  been  Improved 
by  better  matching  the  contact  finger  strlpllnes  to  the  50  ohm  coaxial 
cables  which  carry  signals  to  room  temperature.  Crosstalk  between  power 
supply  and  output  signals  has  been  Identified  as  a  major  problem  In  high 
speed  testing.  The  crosstalk  In  our  probe  Is  largely  a  result  of  excess 
Inductance  In  making  a  common  ground  contact  to  the  ground  plane  on  the 
chip.  Crosstalk  has  been  reduced  by  about  a  factor  of  four  by  a  modification 
which  reduces  this  Inductance.  Further  reduction  is  required  and  will  be 
accomplished  by  a  new  probe  design  in  which  each  signal  line  has  Its  own 
individual  ground  line  to  the  coaxial  shield. 

4.1.2  A/D  Converter  Test  Electronics  -  Another  Item  of  test  equipment 
designed  and  fabricated  this  year  is  an  electronics  box  for  Interpreting 
A/D  converter  output.  This  box  allows  a  direct  plot  of  the  input-output 
characteristics  as  shown  In  Fig.  2  of  Appendix  A.  It  consists  of  (1)  four 
low-level  voltage  comparators  to  convert  Josephson  logic  levels  to  CHOS 
levels,  C21  a  Gray  code  to  binary  code  converter,  and  (3)  a  D/A  converter. 

This  device  has  proved  to  be  Invaluable  In  setting  the  drive  conditions  for 
accurate  A/D  converter  operation. 

4.1.3  Deposition  System  Automation  -  The  vacuum  system  In  which  lead  alloy 
films  are  deposited  and  in  which  oxidation  Is  done,  has  been  significantly 
automated.  The  automation  eliminates  many  possibilities  for  human  error  and 
speeds  the  process.  The  greater  speed  of  the  film  thickness  monitor  now  In 
use  also  Increases  the  degree  of  control  over  film  thicknesses.  The  system 
performs  the  following  functions:  the  proper  source  Is  moved  directly  under 
the  substrate,  the  source  Is  heated  and  the  evaporation  rate  Is  controlled,  the 
shutter  Is  automatically  opened  and  closed,  and  the  process  Is  repeated  for 

up  to  a  total  of  four  layers. 

4.1.4  Move  Into  Clean  Room  -  All  of  the  fabrication  equipment  has  been  moved 
from  an  ordinary  laboratory  Into  a  Class  300  clean  room.  The  move  not  only 
eliminates  dust,  but  also  centralizes  the  optical  mask  making,  fabrication, 
and  soon-to-be-installed  electron  beam  fabricator.  The  move  has  required 
significant  time  during  this  contract  year,  but  overall  It  should  prove  to  be 
extremely  beneficial . 

4.2.  Present  Fabrication  Process 

The  present  fabrication  process  Is  slightly  different  from  that  described 
In  the  last  Annual  Report.  Although  it  Is  much  like  one  In  the  literature 
(1),  the  films  are  about  20%  thinner.  The  process  Is  described  here  In  Its 
entirety  and  summarized  In  Table  I. 
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4.2.1  Substrates  -  The  substrates  are  silicon  wafers  which  are  Insulated  with 
S10(21.  These  wafers  have  a  2  inch  diameter  and  0.010  inch  thickness.  They 
are  purchased  from  an  outside  vendor  with  the  S10(2)  already  formed. 

4.2.2  Layer  1 :  Niobium  Ground  Plane  -  Niobium  is  deposited  onto  the  wafer  by 
e-beam  evaporation  in  an  ultra-high  vacuum  system.  The  films  are  nominally 
4000  Angstroms  thick.  They  are  subsequently  patterned  by  selective  etching. 

The  parts  which  are  not  to  be  etched  are  protected  by  a  patterned  layer  of 
Shipley  AZ-1350J  or  AZ-1370  photoresist.  The  etchant  Is  a  mixture  of  1  part 
hydrofluoric  acid,  9  parts  nitric  acid,  and  20  parts  water,  all  measured  by 
volume.  The  etching  Is  done  at  room  temperature. 

4.2.3  Layer  2:  Nb(2)0(5)  Ground  Plane  Insulator  -  The  Nb(2)0(5)  insulation 
is  formed  by  anodizing  the  Nb  ground  plane.  Some  areas  of  the  ground  plane 
must  remain  uninsulated,  however,  to  permit  electrical  contact.  One  might 
think  of  using  Shipley  1350  Series  resist  to  protect  these  regions  from 
anodization.  However,  the  anodizing  solution  is  a  solvent  for  this  resist. 
Therefore,  the  areas  of  niobium  which  are  to  remain  uninsulated  are  coated 
with  1000  Angstroms  of  aluminum  patterned  using  the  lift-off  process.  The 
anodization  is  carried  out  in  a  solution  (2)  of  18  g  or  ammonium  pentaborate  in 
200  ml  of  ethylene  glycol  at  70  C.  The  thickness  of  the  Nb( 2 )0( 5 )  is  approximately 
proportional  to  the  anodization  voltage,  the  ratio  being  in  the  range  22  to  24 
Angstroms/volt.  Layers  of  700  to  2000  Angstrom  thickness  have  been  made.  The 
protective  aluminum  layer  is  removed  in  50%  H(3)P0(4)  at  60  C. 

It  is  our  impression  that  the  present  IBM  process  employs  photoresist  to  prevent 
regions  from  being  anodized.  Such  a  process  might  be  faster  than  the  one 
discussed  here. 

4.2.4  Layers  3,  6,  9,  and  1_1_:  SiO  -  Silicon  monoxide  is  evaporated  in  a  high 
vacuum  system,  with  a  pressure  during  evaporation  of  about  1  x  10”°  mm  Hg. 

The  material  used  is  a  powder  which  is  evaporated  from  a  baffled  source 
designed  especially  for  SiO.  Both  powder  and  source  are  commercially 
available.  The  S10  is  evaporated  at  a  rate  of  about  20  Angstroms/s. 

The  thicknesses  used  are  layer  3,  2500;  layer  6,  3660;  layer  9,  7200; 
and  layer  11  (rarely  used),  about  10000  Angstroms. 

4.2.5  Layer  4:  Auln(2)  Resistors  -  The  material  used  for  resistors  is 
Auln(2),  an  intermetal  lie  which  is  present  In  the  base  electrode  of  the 
junction.  Resistors  consist  of  46  wt.  percent  Au  and  54  wt.  percent  In. 

With  total  film  thickness  of  437  Angstroms,  the  resistivity  is  1  ohm  per 
square  within  a  few  percent.  The  layers  are  deposited  with  the  following 
thicknesses:  In,  165;  Au,  107;  In,  165  Angstroms.  No  difficulties  have 
been  experienced  with  these  resistors.  They  have  shown  negligible 
resistance  change  with  up  to  20  thermal  cycles  and  room  temperature 
storage  for  one  month. 
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Layer 

Material 

Thickness 

(Angstroms) 

Function 

1 

Niobium 

4000 

Ground  plane.  Defines  magnetic  field 
configuration  near  junctions  and  forms 
striplines  from  conductors  above  it. 

2 

Nb(2)0(5) 

700 

Insulator  above  ground  plane. 

3 

SiO 

2500 

Additional  insulator  above  ground  plane. 
Low  dielectric  constant  and  readily 
variable  thickness  permit  good  control 
of  stripline  impedance. 

4 

Auln(2] 

437 

Resistors  for  terminating  striplines, 
voltage  dividers,  damping  interferometers, 
etc. 

5 

Pb  alloy 

Pb  (88%) 

In  (8%) 

Au  (4%) 

3280 

Junction  base  electrode  and  stripline 
conductor. 

6 

SIO 

3660 

Forms  window  through  which  barrier  is 
formed.  Minimizes  anomalous  effects  at 
Junction  edges. 

7 

Pb(x)Q(yl 

20-50 

Insulating  tunneling  barrier  formed  by 
oxidizing  base  electrode. 

8 

Pb  alloy 

Pb  (98% 1 

Au  (2%l 

404Q 

Junction  upper  electrode. 

9 

SiO 

720Q 

Insulator  between  junctions  and  control 
leads. 

IQ 

Pb  alloy 

8Q02 

Control  leads  for  junctions  and  additional 
stripline  conductors. 

11 

SiO 

10QQQ 

Protective  layer  for  entire  circuit. 
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4.2.6  Layer  5:  Pfa-In-Au  Base  Electrode  -  The  base  electrode  Is  formed 
by  a  sequential  evaporation  of  In,  Au,  and  Pb  having  thickness  400,  80, 
and  2800  Angstroms  respectively.  During  the  baking  of  photoresist  for 
the  next  layers,  these  films  intermix  to  form  a  homogeneous  alloy.  The 
heat  treatment  Is  at  70  C  for  25  minutes. 

4.2.7  Layer  7j_  Tunneling  Barrier  -  The  tunneling  barrier  is  formed  by 
oxidizing  the  base  electrode  through  openings  in  the  photoresist  pattern 
for  layer  8  and,  usually,  the  S10  of  layer  6.  The  oxidation  process  Is 
described  in  three  publications  (3-5)  and  other  references  therein.  It 
Involves  using  rf  sputtering  In  oxygen  to  remove  the  surface  of  layer  5. 

Since  the  surface  continually  reoxidizes,  an  equilibrium  between  the 
rate  of  removal  of  material  and  the  rate  of  reoxidation  occurs.  The 
equilibrium  takes  place  at  a  specific  oxide  thickness  essentially  In¬ 
dependent  of  the  time  taken  (beyond  a  few  minutes).  The  oxygen  pressure 
and  the  cathode  self-bias  voltage  are  carefully  controlled  to  produce 
the  desired  thickness. 

4.2.8  Layer  8:  Upper  Electrode  -  The  upper  junction  electrode  Is  formed 
through  a  sequential  deposition  of  Pb,  Au,  and  Pb,  having  thickness 
2400,  40,  and  1600  Angstroms  respectively.  This  deposition  immediately 
follows  the  formation  of  the  oxide  tunneling  barrier  (without  opening 
the  vacuum  system).  Homogenization  of  the  alloy  Is  promoted  through 
baking  of  the  photoresist  In  subsequent  steps. 

4.2.9  Layer  10:  Pb-In-Au  Control  Lines  -  This  layer  is  formed  from  the 
same  alloy  as  the  base  electrode.  However,  It  is  thicker  to  provide 
good  edge  coverage  over  the  Intervening  layers.  The  presence  of  indium 
In  this  electrode  promotes  Its  connection  to  other  lead  layers.  If 
Indium  were  not  present,  connections  to  other  layers  would  be  Josephson 
junctions,  rather  than  the  superconducting  shorts  which  are  desired. 

The  thicknesses  used  are  In,  976,  Au,  196,  and  Pb,  6830  Angstroms  respectively. 

4.2.10  Design  Rules  Usually  Used  -  Detailed  design  rules  have  not  been 
worked  out.  However,  several  practices  seem  plausible  and  are  usually 
followed. 

4.2.10.1  Step  Coverage  -  It  has  been  found  that  step  coverage  Is  adequate 
If  every  layer  Is  at  least  4Q0  Angstroms  thicker  than  steps  it  must 
cover.  This  procedure  provides  adequate  yield  for  the  circuits  we  have 
fabricated  to  date. 

4.2.10.2  Metal  Layer  Interconnections  -  It  is  thought  that  Interconnections 
between  metal  layers  are  superconducting  only  if  there  is  Indium  In  the 
upper  layer.  The  Indium  combines  with  any  residual  oxide  at  the  interface, 
producing  a  short.  Even  with  careful  sputter  cleaning  of  the  Interface 
prior  to  deposition  of  the  upper  layer,  shorts  are  not  formed  without 

the  upper  layer  Indium.  Thus  there  can  be  no  Interconnections  (except 
Josephson  junctions)  between  the  counterelectrode  and  the  base  electrode. 
Connections  between  the  counterelectrode  and  resistors  are  similarly 
avoided.  Thus  connections  may  be  formed  between  resistor  and  base  or 
control  electrode,  and  they  may  be  formed  between  counterelectrode  and 
control  electrode.  Connections  to  the  ground  plane  are  from  the  base 
electrode,  and  Infrequently  from  the  control  layer.  No  other  Interconnections 
are  permitted. 
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4.3.  Process  for  Fabricating  Resistors  at  One  Ohm/Square 

It  Is  useful  to  be  able  to  deposit  resistors  of  sufficiently  high 
resistivity  so  that  the  length-to-wldth  ratio  of  most  resistors  can  be 
fairly  low.  Since  the  usual  range  of  resistance  we  expect  will  be 
between  1  and  100  ohms,  a  resistivity  of  1  ohm/square  Is  satisfactory. 

The  one  reference  available  (6)  at  the  time  the  work  was  done  describes 
resistors  having  thickness  of  2850  Angstroms  and  resistivity  of  0.045 
ohms/  square.  The  material  Is  Auln(2)  deposited  In  a  layered  form  In 
the  order  In-Au-In.  We  have  tried  this  prescription  and  achieved  a 
resistivity  of  0.03  ohms/square.  In  order  to  Increase  the  resistivity 
we  decrease  the  film  thickness  while  maintaining  the  ratio  of  In  to  Au. 
Achieved  values  of  resistivity  as  a  function  of  thickness  are  shown  In 
Fig.  1.  The  curve  Is  not  quite  smooth,  probably  because  of  errors  due 
to  our  quartz  crystal  film  thickness  monitor.  A  sample  error  bar  Is 
shown  for  the  left-most  point.  No  changes  In  resistivity  greater  than  2 
percent  with  room  temperature  storage  or  thermal  cycling  have  been 
observed  during  a  20  day  storage  period  and  up  to  23  thermal  cycles.  It 
appears  possible  to  fabricate  resistors  having  values  within  ten  percent 
of  that  desired.  Often  results  are  within  a  few  percent.  Indeed  the 
resistors  are  the  most  repeatable  part  of  our  fabrication  process. 

A  new  film  thickness  monitor  has  been  Installed  in  order  to  sample  the 
film  thickness  more  often  (every  0.2  s  compared  to  every  1  s  before). 

In  this  way  better  control  of  thickness  is  expected. 

4.4  Test  Circuits  Including  7-0R  String 

The  following  test  circuits  were  fabricated  and  evaluated: 

1.  Shorts  from  the  base  electrode  to  the  ground  plane.  Areas 
of  5Q  x  5Q  microns  were  found  to  have  critical  currents  In  excess 
of  40  mA.  The  standard  size  ground  plane  short  now  used  Is  20  x  20 
microns  and  has  a  critical  current  known  to  be  greater  than  5  mA. 
Exact  values  of  the  critical  current  have  not  been  determined  as  such 
measurements  have  not  been  required. 

2.  Numerous  transmission  lines.  It  is  important  to  be  able  to 
predict  the  impedance  of  such  lines.  Impedance  was  calculated 
from 

Z  *  (377/W)*CtSl  +S2+L1  +12  ]*(.S1  /el  +S2/e2 )  )**(! /2 ) 


where  W  Is  the  width  of  the  line,  Sn  Is  the  thickness  of  the  nth 
dielectric  layer,  Ln  Is  the  penetration  depth  of  the  nth  metal 
film,  and  en  Is  the  dielectric  constant  of  the  nth  layer.  This 
formula  assumes  metal  films  which  are  thick  compared  to  their 
penetration  depths  and  wide  compared  to  the  dielectric  thickness. 

As  an  example  a  12  micron  wide,  lead  alloy  line  was  placed  over  a 
niobium  ground  plane  covered  with  1173  Angstroms  of  niobium  oxide. 
The  dielectric  constant  was  assumed  to  be  29,  the  penetration  depth 
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for  the  lead  alloy  was  1370  Angstroms,  and  for  the  niobium  was  860 
Angstroms.  The  calculated  Impedance  was  1.17  ohms  while  that 
measured  was  1.55  ohms,  a  difference  of  about  25X.  However,  when 
an  additional  layer  of  S10,  5000  Angstroms  thick,  was  added  to  the 
niobium  oxide,  the  calculated  Impedance  was  9.30  ohms  while  that 
measured  was  9.52  ohms.  The  assumed  S10  dielectric  constant  was  5. 

The  difference  In  this  case  Is  less  than  3  percent.  These  results 
suggest  an  error  In  the  niobium  oxide  thickness. 

Currently  a  more  sophisticated  formula  Is  used  to  calculate  Impedance. 
The  formula  was  obtained  from  Wen  Chang  of  IBM  Research,  and  Is 
about  to  be  published. 

3.  Junction  matched  to  external  transmission  lines  with  50  ohm 
resistors  In  order  to  evaluate  risetime  of  probe.  The  probe  rise¬ 
time  was  found  to  be  about  50  ps  for  a  signal  travelling  on  to 

or  off  a  chip. 

4.  Bridge-type  Interferometer  having  two  control  lines.  This 
device  was  the  first  fabricated  which  permitted  logic  to  be  per¬ 
formed.  The  functions  AND,  OR,  exclusive  OR,  and  NOT  were  per¬ 
formed  at  25  MHz  and  provided  our  first  experience  with  the  dif¬ 
ficulties  of  high  speed  testing. 

5.0  ANTICIPATED  FUTURE  PROCESS  CHANGES 

1.  A  new  upper  electrode  alloy  has  appeared  In  the  literature 
(7,8).  The  alloy  Is  29  wt.  percent  bismuth  In  lead,  deposited  by 
evaporation  from  a  premixed  source.  This  material  Is  reputed  to 
provide  reduced  excess  current,  greater  yield,  and  greater  cor¬ 
rosion  resistance.  Initial  experiments  Involving  this  alloy  were 
unsuccessful,  producing  only  uncontrollably  high  critical  current 
densities.  Further  work  Is  planned. 

2.  The  use  of  thinner  films  (9)  should  Improve  cycling  durability 
and  Improve  edge  coverage. 

3.  A  niobium  anodization  process  using  photoresist  instead  of 
aluminum  would  speed  that  part  of  the  process. 
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Multiple-quantum  interference  ia  a  superconducting  interferometer  is  naad  tor  analog-to-digital 
conversion.  The  simple  felly  parallel  four-bit  converter  which  is  described  it  the  first  known  uee 
of  thi»  effect  m  a  digital  circuit.  Sampling  ram  of  2x  ltf  par  second  ware  achieved,  and  much 

a-;_a- - - «■*. 
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PACS  uaibm:  06.60JH,  74.3a  +  r,  06.70.Dn,  06.70.Ep 


This  latter  describes  the  Ant  known  use  in  a  digital  cir¬ 
cuit  of  the  multiple  lohes  of  the  threshold  curve  of  a  super¬ 
conducting  interferometer.  Although  digital  applications  of 
such  interferometers  have  been  known  for  several  years,14 
only  the  ftrst  lobe  has  been  used  previously.  An  analog -to- 
digital  (A/D)  converter  using  multiple  lobes  has  been  pro¬ 
posed.'-4  We  demonstrate  a  realisation  of  this  concept  in  a 
four -bit  A/D  converter  operating  at  a  rate  of  200  megasam- 
pies  per  second.  An  increase  to  rates  well  above  one  gigasam- 
pk  per  second  may  be  poaaibk. 

This  experiment  makes  use  of  four  current  compara¬ 
tors,  each  consisting  of  a  Joeephson  three-junction  interfer¬ 
ometer  which  has  the  junctions  connected  in  parallel  by  su¬ 
perconducting  loops.'  A  photograph  of  one  of  the 
interferometers  is  shown  in  Fig.  I  (a).  The  center  junction 
has  twice  the  critical  current  of  each  of  the  symmetrically 
arranged  outer  junctions.  Such  an  interferometer  has  a  cur¬ 
rent-voltage  characteristic  like  that  of  a  single  junction.  Its 
maximum  0-V  current  lm  (threshold  current)  is  modulated 
by  an  applied  magnetic  Held  in  a  way  determined  by  the 
junction  critical  currents  and  the  loop  inductances.  The 
junctions  are  sufficiently  small  (about  4pm  square)  that  the 
magnetic  fields  used  have  negligible  effect  on  them  directly. 
By  choosing  Li0  s  jd0 .  a  threshold  curve  like  that  shown  in 
Fig.  2(a)  is  obtained.  Here  £  is  the  inductance  of  one  of  the 
superconducting  loops,  i0  is  the  critical  current  of  one  of  the 
outer  junctions,  and  40  is  the  flux  quantum  h  /It. 

The  periodicity  of  the  threshold  curve  is  a  manifestation 
of  quantum  interference  occurring  in  the  interferometer. 
The  net  phaae  change  of  the  wave  Auction  describing  the 
superconducting  loops  in  the  interferometer  incresses  in 
proportion  to  the  amount  of  magnetic  flux  applied  through 
the  loops  of  the  interferometer.  A  net  phase  change  of  2ir 
corresponds  to  one  additional  flux  quantum  penetrating  the 
loop.  The  threshold  curve  repeats  whenever  two  additional 
flux  quanta  are  added.  In  a  two-junction  interferometer,  the 
curve  would  repeat  for  each  additional  flux  quantum. 

To  use  this  interferometer  as  a  current  comparator,  a 
magnetic  field  is  applied  equally  through  both  loops.  The 
magnetic  field  results  from  passing  a  current  /c  through  a 
control  line  above  the  interferometer.  Its  effect  is  illustrated 
by  the  horizontal  arrow  in  Fig.  2(a).  The  current  through  the 


interferometer/,  is  then  raised  from  zero  to  a  predetermined 

level/, indicated  by  the  vertical  arrow.  If/, is  above 

the  threshold  curve  the  interferometer  switches  to 
the  voltage  state;  otherwise  it  does  not 

A  single  interferometer  is  used  to  gene  rate  a  single  bit  of 

the  output  of  an  A/D  converter,  by  choosing  /, _ so  that 

ths  interferometer  switches,  or  does  not  switch,  for  equal 
intervab  of  control  current  By  — the  control  current 
the  analog  signal  to  be  digitised,  the  voltage  output  of  the 
gate  V,  is  one  bit  of  the  digital  value  of  the  signal  as  illustrat¬ 
ed  in  Fig.  2(b).  To  take  repeated  samples,  /,  is  cyclically 

varied  between  0  and  /, _ 

The  next  most  significant  bit  can  be  generated  by  ex¬ 
panding  the  bit  pattern  of  the  least  significant  bit  by  a  factor 
of  2.  The  expansion  is  produced  by  reducing  the  coupling  of  a 
second  interferometer  to  the  control  line  by  a  factor  of  2. 
Similarly,  the  third  and  fourth  bits  are  generated  by  farther 
expanding  the  bit  patterns  of  third  and  fourth  interferomet¬ 
ers.  Thus  only  four  interferometers  ate  required  for  A/D 
conversion  of  four-bit  accuracy.  Moreover,  since  these  inter¬ 
ferometers  operate  in  parallel,  maximum  speed  is  achieved 
with  all  four  bits  being  produced  essentially  simultaneously. 
In  general,  only  AT  interferometers  are  required  to  digitize  to 
Af-bit  accuracy,  with  little  decrease  in  speed  needed  for  addi¬ 
tional  bits. 


•’Work  supported  in  pan  by  the  Office  of  Nsvtl  Research  under  contract  FIO.  I.  («)  Photograph  of  one  of  the  four  interferometers  in  tin  A/D  cow- 
number  NOOOI4-77-F-4S.  verier,  (b)  Photograph  of  s  complete  four-bit  convener. 


7» 


Appf.  Phys.  Lett.  3«9).  1  November  1979  0003-4951  /79/210720-08900.90 


7J0 


FIO. 2. (a) Tbnahoid curve thenriaf naunam  zero voha§e cwnM lm  m 

a  turner  -a  of  control  currant  /, .  Idwl  manitnum  curvt  currant  /, it 

thorn  by  daahad  bar  (b)  Output  vnbfr  for  idol  /, _ «*  function  of 

control  current,  (c)  Actual  bit  patent  for  tbi  four-bit  A/D  ccuvartar.  (d) 
Value  of  digital  output  of  eoautrur  varaaa  analog  ogaal  (mead). 


Figure  1  (b)  show*  a  complete  four-bit  converter.  The 
observed  bit  patterns  for  this  converter  are  shown  in  Fig. 
2(c).  This  pattern  is  a  Gray  code  representation  of  the  analog 
signal.  Gray  code  is  frequently  used  in  A/D  conversion  since 
only  one  bit  changes  as  the  signal  crosses  from  one  digital 
level  to  the  next  Another  significant  feature  of  this  A/D 
convertor  can  be  realized  by  adding  a  second  control  line  to 
each  interferometer.  By  applying  a  dc  bias  through  ooe  of 
these  control  lines,  the  bit  pattern  from  that  gate  can  be  shift¬ 
ed  to  values  of  higher  or  lower  current  Using  this  feature,  a 
variety  of  unique  codes  can  be  produced. 

Superconducting  A/D  converters  of  the  type  described 
have  been  fabricated  with  two  converter*  on  each  6.4-mm- 
square  chip.  The  interferometers  are  much  like  those  de¬ 
scribed  in  Ref.  5,  although  they  are  scaled  up  in  linear  size  by 
a  factor  of  4  to  make  fabrication  simpler.  The  circuits  con¬ 
tain  resistors  for  damping  the  interferometer  resonances  and 
for  terminating  transmission  lines.  The  impedances  of  the 
striplines  on  the  chip  are  choecn  for  proper  matching  to  the 
interferometer  impedances.  The  coupling  of  the  signal  to  the 
interferometers  is  varied  by  adjusting  the  length  of  the  con¬ 
trol  line  above  the  interferometers  in  the  ratios  1:2:4: 8.  As 
can  be  seen  in  Fig.  1,  the  control  line  is  diverted  around  the 
portion  of  the  wide  couatenlectrode  passing  over  the  inter¬ 
ferometer.  In  this  region  the  counterelectrode  provides 
shielding  which  is  not  readily  calculable  Furthermore,  this 
control-line  arrangement  causes  the  coupling  to  each  half  of 


each  interferometer  to  be  insensitive  to  slight  —i— ligmna — 
during  fabrication. 

The  interferometers  have  been  tested  with  a  gate  cur¬ 
rent  having  a  frequency  (Le.,  sampling  rate)  of  200  MHz. 
Due  to  equipment  limitations  only  two  of  the  four  bits  were 
tested  at  this  speed.  They  were  found  to  produce  outputs 
which  repeat  synchronously  with  a  nominally  30-MHx  sig¬ 
nal  In  this  sense,  the  converter  has  been  daaMaatraiad  to 
operate  at  200  megaaamples  per  second.  Baaed  on  the  ob¬ 
served  switching  times,  we  believe  that  the  devices  fabricated 
will  operate  properly  at  rates  up  to  about  700  megasampliu 
per  second.  Improving  the  ground  plane  contacts  should  al¬ 
low  the  converter  to  operate  at  much  higher  rates.  Reported 
switching  times’  of  42  ps  for  these  devices  suggsst  rates  well 
above  1  gigasample  per  second. 

To  assess  the  accuracy  of  the  full  device,  conventional 
circuitry  was  used  to  convert  from  Gray  code  back  to  an 
analog  signal  The  resulting  output-vs-input  curve  [Fig. 
2(d)]  is  monotonic  with  a  maximum  error  in  step  width  of  0.7 
least  significant  bits.  These  measurements  were  performed 
at  about  10  megasamples  per  second  with  the  signal  swept  at 
a  few  hundred  hertz. 

The  coupling  of  the  signal  to  the  interferometers  was 
determined  to  be  in  the  ratio  13:4:8  within  about  1%,  indi¬ 
cating  that  this  design  may  be  suitable  for  converters  of  6  to  8 
bits. 

The  initial  tests  of  this  A/D  converter  design  suffer 
from  variations  in  interferometer  characteristics,  trapped 
flux,  and  the  general  difficulties  of  high-speed  testing.  In 
addition  the  ultimate  accuracy  of  the  device  may  be  limited 
by  dynamic  effects  including  subthreshold  swindling  and 
signal  line  transients  earned  by  switching.  Nevertheless,  the 
design  appears  to  offer  great  promise  for  significantly  higher 
speeds  in  the  ftiture. 

It  is  a  pleasure  to  acknowledge  hdpAil  conversations 
with  T.  Fulton,  L.  Geppcrt,  D.  HcrreU,  and  H.  Zappe.  H 
Peterson  assisted  with  the  understanding  of  the  threshold 
curves,  D.  McDonald  provided  constant  encouragement, 
and  L.  Mullen  assisted  with  the  febrication. 
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Threshold  curves  for  multqunction  superconducting  interferometers  have  been  calculated 
previously,  showing  general  agreement  with  observed  features,  especially  in  symmetric  caeet,  We 
here  add  some  more  details  to  the  analysis,  paying  particular  attention  to  the  effects  of 
asymmetries  in  coupling,  inductance,  or  critical  currents.  Feed-loop  inductance  and  flux 
quantization  in  the  feed  loop  can  be  important  A  changing  lobe  pattern  over  many  periods, 
asymmetries  within  a  period,  shifting  patterns  between  runs  a  warm-up,  and  sudden 

changes  in  pattern  bemuse  of  noise  in  the  environment  are  all  quantitatively  explainable  on  the 
basis  of  this  modeL  By  use  of  a  single  “calibration  curve”,  the  inductance  for  symmetric  two-  or 
three-junction  interferometers  can  be  obtained  immediately. 


PACS  numbers:  74.50.  4-  r,  85.25.  4-  k 
I.  INTRODUCTION 

A  superconducting  interferometer — a  device  contain¬ 
ing  two  or  more  Josephson  junctions  connected  in  supercon¬ 
ducting  loops — is  attractive  for  use  in  logic  circuits  because 
of  its  high  sensitivity,  high  speed,  good  current  transfer,  and 
small  heat  dissipation.  It  can  also  be  modeled  quite  accurate¬ 
ly,  as  we  hope  to  emphasize  in  this  paper.  One  characteriza¬ 
tion  of  such  a  device  is  its  threshold  curve,  /„  (ft)  or  /„  (/,), . 
where/,,  is  the  maximum  bias  current  the  device  can  accom¬ 
modate  before  developing  a  voltage,  and  Ie  is  a  control  cur¬ 
rent,  whose  magnetic  induction  B  threads  the  loops  of  the 
interferometer.  The  calculation  of  these  threshold  curves,  or 
certain  properties  of  them,  has  been  the  subject  of  many 
papers.  1*w 

The  purpose  of  the  present  paper  is  to  extend  the  analy¬ 
sis  of  the  threshold  curves  in  several  directions.  The  empiri¬ 
cal  curve  obtained  by  DeWaele  and  DeBruynOuboter*  for  a 
double-junction  interferometer,  which  provides  a  method 
for  the  rapid  determination  of  the  inductance  of  symmetric 
interferometers,  is  shown  to  be  valid  for  the  three-junction 
interferometer  as  well.  Approximate  analytical  expressions 
for  this  curve  are  developed.  We  show  the  importance  of 
including  parameters  of  the  feed  loop— both  its  inductance 
and  passible  trapped  flux — in  the  analysis.  The  effects  of 
various  asymmetries,  whether  in  junction  critical  currents, 
loop  inductances,  or  flux  coupling,  are  discussed  in  some 
detail.  We  also  summarize  the  salient  features  of  the  thresh¬ 
old  curves,  many  of  which  have  been  mentioned  before  but 
not  in  one  place.  Finally,  we  make  comparisons  with  some 
threshold  curves  measured  for  three-junction  interferomet¬ 
ers  with  slightly  asymmetric  coupling. 

IL  THEORY 

The  technique  of  constraint  maximization,  which  is  the 
calculations!  method  we  have  used,  has  been  described  ade¬ 
quately  before.”  Briefly,  one  wishes  to  calculate  the  maxi¬ 
mum  or  threshold  value  Im  of  the  gate  current 

l,  «  rind, ,  (1) 

which  can  flow  through  all  the  junctions  before  a  voltage 
develops.  In  Eq.  (1),  4,  is  the  phase  difference  across  the  /th 


junction  of  critical  current /„,,  and  the  summation  is  over  all 
junctions  in  the  interferometer.  This  relation  is  subject  to  the 
constraints  imposed  by  fiuxoid  quantization  within  each 
loop,  which  may  be  written,  for  loop  k, 

fk  +  £(£,/,  +  d,*o/2ir)  -  ijf»# o 

-0.  (2) 

In  Eq.  (2),  <txk  is  the  external  flux  in  loop  k,  whether  inten¬ 
tional,  as  from  a  control  line,  or  accidental,  as  from  flux 
trapped  nearby  in  a  ground  plane.  Also,  /,  is  the  current 
flowing  through  inductance  L,,  <P0  »  2.0685  X 10“ 15  Wb  is 
the  flux  quantum,  and  A,  is  an  integer.  The  summation  is 
over  all  inductances  and  junctions  in  loop  k,  and  may  also 
include  mutual  inductances  from  other  loops.  To  And  the 
extrema  of  /,,  one  uses  the  technique  of  Lagrange  multipli¬ 
ers,  forming 

Hd. ,  #, ,  ,  ...)=/,  +  £2*  /*  ,  (3) 

and  requiring  dr  /d> l,  *  0  and  dT  /d+,  —  0.  The  first  of 
these  conditions  simply  reproduces  Eq.  (2).  The  second  cre¬ 
ates  a  set  of  equations  which  overdetermine  the  Xk ,  thus 
imposing  relations  among  the  d,- The  sets  |d, ,  di.  — )  con¬ 
sistent  with  these  relations  are  the  allowed  solutions.  The 
threshold  current  Im  is  found  by  substituting  these  values 
into  Eq.  ( 1 ).  Generally  the  sum  of  all fk  —  0  equations  results 
in  a  convenient  equation  into  which  the  solution  set 
[di,di>  — 1  is  substituted,  thus  giving  the  value  of/c  or  total 
flux  d,  to  be  paired  with  the  calculated  Im. 

This  technique  yields  unstable  as  well  as  stable  solu¬ 
tions.  The  unstable  points  can  be  rejected  by  evaluating  the 
Hessian  of  the  potential  energy*-11  (second  derivatives  of  the 
Gibbs  free  energy  G  relative  to  principal  axes)  and  retaining 
only  the  positive  values.  Although  relatively  simple  for  two- 
junction  interferometers,  evaluation  of  the  Hessian  can  be 
quite  involved  for  three-junction  cases,  especially  when  feed 
loops  are  considered.  Landman*  used  an  alternative  approxi¬ 
mate  procedure.  We  have  found  that  use  of  the  necessary, 
but  not  sufficient,  condition 

1  +  (LI,*/# )  oosd,  >  0 ,  (4) 

will  reject  most  of  the  unwanted  points.  Here  and  subee- 
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quently  we  use  the  convenient  “phi-ber”  symbol  G—  <Pa 
flv.  The  inductance  L  in  Eq.  (4)  is  typically  the  inductance 
through  which  the  junction  current  flows.  Equation  (4)  is 
derivable  from,  but  not  generally  identical  to,  d*G  /d+  f  >  0 
(Le^  the  derivative  condition  before  transformation  to  prin- 
cipal  axes).  In  special  cases,  the  latter  is  simple  to  evaluate 
and  is  an  improvement  over  Eq.  (4).  However,  we  find  that 
the  simple  Eq.  (4)  is  usually  quite  satisfactory. 

III.  TWO-JUNCTION  INTERFEROMETER 

Various  facets  of  the  threshold  curves  for  two-junction 
interferometers  have  been  given  earlier.1**4*1*"’  Fulton  «t 
at."  and  Schulz-DuBois  and  Wolf'  have  given  rather  de¬ 
tailed  analyses,  and  we  shall  add  just  a  few  new  points  here, 
while  repeating  some  earlier  material  bearing  on  measurable 
quantities. 

For  the  twojunction  interferometer  as  depicted  in  Fig. 
1(a),  the  equations  are 

Cj  »  —  Cx/(Im/IQt  +0C,),  (5) 

/„  »/oiS|  +/<a-S j  ■  (6) 

+  &S, -4x  -0xSt  +2vN ,  (7) 

where 5,  m sin#,, C,  mcat4,.0*0i  +0,Iat/IOi, 

0,  »  LJv/t,  and  6  -  .  Equations  (5)— (7)  show  that 

Im  is  periodic  in  $  with  period  2w.  Commonly,  the  phases#, 
ate  restricted  to  the  range  —  r<4i  <«*,  and  the  integer  .V  of 
Eq.  (7)  is  changed  to  obtain  the  contour  of  a  different  fluxoid 
quantum  state.  However,  there  are  ranges  of  parameter  val¬ 
ues,  such  as  observation  (3)  below,  where  this  identification 
of  fluxoid  states  is  not  very  meaningful,  as  Schulz-DuBois 
and  Wolf  "  also  point  out  Clearly,  one  can  alternatively  set 
N  m  0  and  allow  the  4,  to  span  an  extended  range. 

We  now  make  a  few  observations  about  double-junc¬ 
tion  interferometers  which  relate  generally  to  measurable 
features,  and  from  which  the  junction  parameters  may  be 
deduced. 

(l)4i  m4tm  ±4**  are  always  solutions,  even  for 
asymmetric  feeds  or  critical  currants,  as  Eq.  (3)  shows.  For 
these  phases,  !m  -  ±  (/«,  4-  /«).  Thus  the  maximum  cur¬ 


rent  is  always  the  sum  of  the  junction  critical  currants,  which 
is  not  the  case  with  the  three-junction  interferometer.  Also, 

±<fc-*.)  +  2wtf,or#,«  ±  (£,/„-£,/„) 

+  N0O,  at  these  extrema.  Measurement  of  the  relative  dis¬ 
placement  of  the  minimum  and  maximum  of  Im  thus  pro¬ 
vides  one  measure  of  the  asymmetry.  The  three  curves  of 
Fig.  2  show  this  displacement.  These  points  have  been  made 

(2)  An  asymmetry  of  the  feed  injection  point  (£,  #£,) 
while  the  junction  critical  currants  are  kept  equal  will  leave 
the  cusp  ends  of  the  lobes  on  the  /„  —  0  axis,  ss  curve  (a)  of 
Fig.  2  shows.  Where  there  is  an  asymmetry  of  the  critical- 
currants,  these  ends  will  lift  off  from  the  Im  —  0  axis  to  near 
the  values  ±  (/0,  — 7W).  These  currants  are  actually  the 
points  of  zero  slope,  as  evaluation  of  dlm/d$  shows,  and 
occur  at,  e.g.,  di  »4»*withdj  “  ^  jo*,  giving 

^(r  +  0i  +  0,)  +  2vN.  Curve  (b)  of  Fig.  2  and  the  in¬ 
set  illustrate  this  effect  In  principal,  a  dynamic  experimental 
procedure  might  be  used  to  locate  these  extrema.  These  to¬ 
gether  with  the  measured  maximum  of  Im  completely  speci¬ 
fy  the  critical  currents.  Reference  1 1  supplies  additional  de¬ 
tails  about  the  nature  of  the  cusps. 

(3)  Both  from  experiments  and  graphical  analysis,  Ful¬ 
ton  tt  ai.'  noted  that  double-junction  interferometers  can 
exhibit  nearly  sinusoidal  Im  (Ic)  behavior.  This  result  is  easi¬ 
ly  deduced.  When 

/<h//oi  +0<l.  (8) 

Eq.  (5)  has  a  solution  for  di  only  for  a  limited  range  of  d, 
near  £*,  r,  etc.  Setting  4,  «  Js-  and  noting  that  Eq.  (8) 
requires /Sj  <1(0,  is  not  required  to  be  small  however),  we 
getdj~0  +0,  +  4**  from  Eq.  (7).  Hence  Im 
m/oi  +  /<n  a"(0 + 0, ),  showing  the  sinusoidal  behavior. 
Similarly,  with  4\  *  —  4»".  one  gets  Im 
a ■  — /0,  —lot  coHd—0,).  A  symmetric  interferometer 


FIO.  2.  Thmhold  curv  for  thrw  doubte-pntctwn  iatttftrottwi  Inving 
the  ium  tool  indocuact  and  total  critical  cumnl  (a)  L,  -0.9  pH. 

Li  —0.1  pH ,/«  —  An  —  0.5  nA.  (b)  L,  —  L,  —  0.5  pH,  An  —  0.7  mA . 
4i  —  0.3  nA.  (c)  i|  —  L,  —  0.5  pH,  A,,  —  0.9  nA.  An  —  0.1  mA. 
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FIG-  3.  Ttaohold  curves  for  two  symmetric  double- junettoo  interferomet¬ 
ers.  (■)£  *0.2  pH./,,  —  0.5  mA.fmng/?«  0.30.  (b)£  —  4.0pH,/„  —0.3 
mA.  fivin|£  —  6.1. 


can  never  exhibit  this  behavior  since  it  cannot  satisfy  Eq.  (8). 
Curve  (c)  of  Fig.  2  shows  that  the  threshold  curve  can  be 
quite  sinusoidal  even  when  Eq.  (8)  is  not  particularly  well 
satisfied  (/OT  //01  +0  =*0.41). 

(4)  When  the  interferometer  is  symmetric,  thane  is  al¬ 
ways  a  slope  discontinuity  midway  between  the  lobes  at 
±1*0.  etc.  The  interferometer  is  now  governed  by  a 
single  internal  parameter  A  which  has  the  value  of  U0/$ 
where Z.  =»  L,  ■+■  L ,  is  the  total  loop  inductance  and  /0  is  the 
common  junction  critical  current.  Figure  3  shows  two 
threshold  curves  for  symmetric  junctions  with/0<l  and 
/?>!.  The  threshold  current  at  the  midpoints  is  a  convenient 
measure  of  0.  Curves  of  this  current  versus  0  based  on  nu¬ 
merical  calculation  have  been  given.’ 11  Although  an  exact 
formula  for  this  curve  has  not  been  determined,  formulas  in 
the  two  limits  can  be  established  fairly  easily.  Figure4shows 
our  calculated  curve  together  with  some  analytical  results  in 
limiting  cases. 

For/?<l,  Eq.  (3)  shows  that  di  and  are  separated  by 
about  w  rad.  Linearizing  in  0,  one  finds 

d2  *  +0C  J  /S\  . 

Substitution  of  this  d2  into  Eq.  (7)  with  9  *  ir  gives  C  ? 
«S{.  Equation  (6)  then  results  in 

0<t-  (9) 

When 1,  Eq.  (3)  shows  that  C,  m  -  1/0  whence,  >1. 
Substitution  ofd2  *  i *+  1/0  +  —  into  Eq.  (7)  with  <?  —  rr 
then  gives 

*x  +(^)S,  •V-iir  +  O<0-'), 
which  gives  *  J*’  —  (4ir/0)l/i  +  0(0  ~ ').  Hence 

|/.m//— 00) 
This  is  equivalent  to  a  modulation  depth  of  <£<,/•£,,  a  result 
noted  earlier.1’ 


However,  a  much  better  result  can  be  obtained  at  fol¬ 
lows,  and  now  we  generalize  to  permit  asymmetric  feed,  al¬ 
though  retaining  70l  «/«.  At  large#,  the  lobes  degenerate 
toward  parallelograms.  Noting  where  the  corners  lie  from 
the  discussion  of  observations  (1)  and  (2),  one  readily  finds 
the  equations  for  the  straight-line  segments  of  these  parallel¬ 
ograms.  The  intersections  lying  highest  (there  can  be  more 
than  one  set  of  intersections  at  large/?,  as  Fig.  3  shows)  give 

| -£/(*■+*).  0>i  (ID 

and 

9^~(2n+l)ir+0(0l-0ly{v  +  0).  (12) 

Note  that  the  current  at  the  intersection  is  a  function  of  0 
only,  i.e.,  independent  of  the  asymmetry,  in  this  approxima¬ 
tion.  Fulton  et  al.'  and  Clarke  and  Paterson*  had  earlier  no¬ 
ticed  this  independence  from  their  numerical  work.  Figure  4 
shows  that  this  formula  gives  a  better  result  than  the  preced¬ 
ing,  and  in  fact  is  a  fair  approximation  to  the  correct  curve 
through  the  entire  range  of  0.  It  is  a  useful  analytical  result 
since  it  allows  a  rapid  determination  of  the  inductance  L 
from  the  threshold  curve. 

An  alternative  circuit  for  a  double-junction  interferom¬ 
eter  is  one  in  which  the  control  current  is  injected  directly 
into  the  interferometer,  as  shown  in  Fig.  1(b).  Here,  one 
must  know  what  fraction  of  the  total  inductance  is  coursed 
by  the  control  current,  as  illustrated  in  Fig.  1(b).  It  is  this 
fraction  which  determines  the  period  in  Ie .  The  of  the 
above  formulas  is  now  replaced  by  ( p ,  L,  ■¥p1L1)Ic.  The 
remaining  formulas  are  unaffected.  Similarly,  when  the  con¬ 
trol  line  overlies  the  plane  of  the  interferometer,  the  coupling 
is  not  generally  perfect.  When  the  widths  of  all  lines  are 
much  greater  than  the  distances  separating  them,  and  film 
thicknesses  are  much  greater  than  their  penetration  depths, 
the  control  current  can  be  treated  as  if  directly  injected  and 
passing  through  all  of  the  inductance  of  that  length  of  line 
overlain  by  the  control  line.  Since  these  conditions  may  not 
be  met  in  some  practical  cases,  the  degree  of  coupling  may  be 
somewhat  uncertain,  and  only  a  fraction  of  the  inductance  of 


FIG.  4.  Curve  of  the  threshold  current  midway  between  maxima  *s  a  hinc- 
tionof  0mUa/$  (solid  curve).  The  curves  an  valid  for  symmetric  two-  or 
three-junction  interferometer!  t»  defined  in  the  text. 
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FIG.  5.  Circuit  diatom  for  a  thicejuactioa  interferometer  with  ■  feed  loop. 
Magnetic  fluxes  and  direct  injection  aC  control  current  are  shown  together. 


the  overlain  line  is  involved.  In  these  cases  the  period  of  Ie  in 
the  lobe  pattern  is  not  a  reliable  measure  of  inductance,  but  is 
however  a  useful  measure  of  the  degree  of  coupling  between 
control  line  and  interferometer.  Measurement  of  Im  at  its 
lowest  value  between  the  lobes  is  superior  for  measurement 
of  inductance  in  symmetric  interferometers. 

IV.  THREE-JUNCTION  INTERFEROMETER 

An  example  of  a  three-junction  interferometer  circuit  is 
shown  in  Fig.  S,  in  which  a  control  current  is  injected  direct¬ 
ly  into  the  interferometer,  on  the  side  containing  a  feed  loop. 
Since  the  feed  loop  is  entirely  superconducting,  flux  quanti¬ 
zation  within  it  must  be  included  in  the  analysis.  This  may  be 
written 

L.I.  +L"I"  —L  ;i;  -  LbIb  +  4.f  -  .  (13) 

where  we  have  used  the  notation  indicated  in  Fig.  5.  iM/  is 
feed  loop  flux  from  sources  other  than  the  control  current, 
for  example  from  flux  trapped  in  the  ground  plane,  and  Nf  is 
the  number  of  flux  quanta  trapped  in  the  feed  loop.  One  may 
eliminate,  e.g., and  1 7  in  favor  of  I ",  Ic ,  I, ,  /, ,  and  /3 
by  using  the  Kirchoff  current  equations,  and  solve  Eq.  (13) 
for  /  7-  Substitution  into  the  two  fluxoid  quantization  condi¬ 
tions  for  the  two  loops  containing  the  junctions  then  gives 

f  =*^xi  +  «id*/  +  !  +«i(£.  +£*)& 

+  (^l  +£|  +'<*|f"a)fj  “(Jf+B|ft)/j 
-  a,  (L  ?  +  Lb)I}  +  *(d,  -  *3)  -  W  +  a,  N,)4>0 
-0.  (14) 

A  * 4xl  +  a2d*/  +  [LJ  +a2(I.a  +£*)]/, 

+  (M  +  a2Lb)I,  —{Lj  +L  J 
+  a:(L,  + 1  7)/j  +  *(d,  -  d2)  -  W  +  <*:*,)*„ 

-0,  (15) 

where  a,  «  L  "/Lf,  a2  »  L  2/Lf,  and  Lf~Lm+  Lb 
+  L  "  +  L  x- .  4*1  snd  4mi  have  the  same  meanings  as  de¬ 
scribed  for  4,f  above. 

Following  Tsang  and  Van  Duzer,'  we  have  included  a 
“mutual  inductance”  M,  representing  coupling  between  the 
two  loops.  One  sees  from  the  terms  a,  Lb  and  a2L,  added  to 
Af  that  the  feed  loop  also  facilitates  coupling  between  the 
loops.  The  terms  containing  Nf  show  that  flux  trapped  in  the 
feed  loop  adds  fractionally  to  the  flux  in  the  other  loops. 
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The  “Lagrangian  current”  is  now  formed  according  to 
Eq.  (3).  The  equations  dr  /34,  “  0,  with  /  m  1, 2, 3  give 
three  equations  containing  the  two  Lagrange  multipliers. 
Two  of  these  equations  are  used  to  determine  A,  and/lj. 
When  these  ate  substituted  into  the  third,  a  condition  is  im¬ 
posed  among  the  d, .  which  may  be  written 

C  =  +f<nCi  +  A,2C,Ci) 

/w  +  A,2C,  +AuCi  +  At23C,C2 

where 

A 13  “/q i0\  +foiA  • 

Au  ™  As0i  +  I<n(J  3  , 

A n  *  /ot  (02  +  Vi)  +  +  7i )  > 

Am  -/o.0?iA  +*li0i) 

+  I<n(V \0i  +  0\0i)  +/<»(0!&  ~ViVz)< 

0,  —  (f-i  +•  L  1  +  otL,)I0l /<P , 

02  *  (I2  +  L  j  +  , 

0y  “  ai  (L  j  +  L>Vo,  > 

0i  m  <*!((■  7  +  LMn/*  1 
V\  »(M  +  a1L.)I0,/<P, 

Vi  =*(M  +  a,Lb)I01/<P. 

Equations  (14)  and  (IS)  impose  another  condition 
among  the  4, ,  which  is  found  by  eliminating  Ic  between 
them.  At  this  point,  we  shall  reduce  the  great  generality  of 
the  above  by  imposing  the  restriction 

Lm/Lk»L:/L1, 

which  is  a  common  design  goal  and  which  makes  f$\  —0,, 
and  a\Lm  =*  a,Lb.  Combining  Eqs.  (14)  and  (15)  then  gives 

(1  +  F)B,  ~i?2  +  {N2  <P0  -  drt)/0  -  E, 

+  F[F,  -(N,<P0 -+x,)/<P-E2] 

+  (a2  -  Fa, )  (Nf<Pa  -  *„)/* ,  (17) 

where 

F-  [L  J  +  a2{L,  +  Lb)\/[L  J  +  a, (i.  +  /.»)] , 

B 1  “d,  +0,S, , 

E,  =  V,S<  ■  (18) 

The  coupling  factor  F  is  the  ratio  of  fluxes  from  the  control 
current  entering  the  two  loops  of  the  interferometer. 

Our  procedure,  which  is  more  direct  and  apparently 
simpler  than  those  described  by  Landman*  and  Schulz-Du- 
Bois  and  Wolf,"  is  to  choose  d,  and  d2 ,  and  determine  C3 
from  Eq.  (16)  and  B}  from  Eq.(  17).  The  values  of  C,  and  21 3 
each  specify  a  set  of d2.  Only  certain  pain  f  d, ,  d2 }  allow 
equality  between  mem  ben  of  these  two  sets.  When  the  al¬ 
lowed  di»dj.  and  d3  are  found,  they  are  substituted  into  Eq. 
( 1 )  and  into  Eq.  ( 1 4),  or  Eq.  ( 1 5),  or  the  sum  of  the  latter  two, 
which  is  simpler 

9  “  B2  +  E2  —  B,  —  E,  +  2tr[2V|  +  V2  +  (<X|  +  a2)Afr] , 

(19) 

where 

9=[L\  +Li  +(a,  +«,)(/,.  +L*)]/,/$ 

+  [d*i  +du  +(a,  +02^,/]/^.  (20) 
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Note  that  the  coefficient  of  /,  /#  in  Eq.  (20)  is  the  total 
inductance  through  which  the  control  current  passes,  in¬ 
cluding  that  of  the  feed  loop.  If  the  control  current  were 
injected  in  the  opposite  side  of  the  interferometer,  as  indicat¬ 
ed  by  the  dashed  lines  in  Fig.  S,  the  coefficient  of  Ic  /0  would 
be£,  +  L2 .  None  of  the  other  formulas  are  affected.  In  the 
explicit  flux  term  of  Eq.  (20),  no  assumptions  concerning  the 
apportionment  of  the  flux  from  other  sources  among  the 
separate  loops  is  made;  the  coefficient  of  d,/  results  from  the 
flux  quantization  condition  in  the  feed  loop. 

When  the  control  current  is  not  injected  into  the  inter¬ 
ferometer,  but  lies  above  it,  as  is  usually  the  case,  one  could 
drop  the  Ie  term  in  Eq.  (20)  and  treat  the  d*.  as  due  to  the 
control  current,  and  stray  flux  if  necessary.  This  requires  a 
change  in  some  of  the  terms  in  Eq.  (17).  The  coupling  factor 
Fin  Eq.  (18)  is  the  ratio  of  the  coefficients  of  Ic  in  Eqs.  (14) 
and  (IS).  Dropping  the  Ic  terms  there  requires  taking  the 
control  current  contributions  d  «  to  the  fluxes  as  propor¬ 
tional  to  each  other.  Thus,  setting 

dii  ““ndii  • 
casts  F  into  the  form 

F*  (a2l  +  as x,)/(l  +  as x, ) ,  (21) 

of  which  Eq.  (18)  is  a  special  case.  The  form  of  Eq.  (17) 
remains  intact  if  we  take  the  d«,  appearing  there  as  possible 
stray  fluxes.  a2l  and  ay  are  determined  from  the  geometrical 
arrangement  of  control  lines  over  the  interferometer  and 
feed  loop. 

The  integers  IV,  and  N2  are  immaterial  in  determining 
the  threshold  curves  if  di  and  di  are  allowed  to  spaa  an 
arbitrary  range.  That  is,  in  both  Eqs.  (17)  and  (19),  AT,  and 
appear  only  in  the  combinations  d,  ±  2jhV,  .  The  same  is 
not  true  for  Sf,  the  number  of  flux  quanta  trapped  in  the  feed 
loop.  A  changed  value  of  Nr  will  not  only  translate  the 
threshold  curve  along  the  Ic  axis  in  an  amount  not  generally 
equal  to  a  multiple  of  2tr  [in  the  units  of  Eq.  (19)],  but  will 
also  change  the  shape  of  the  curve  if  a2  /fa, ,  as  Eq.  (17) 
shows.  Stray  flux  will  likewise  both  translate  and  change  the 
shape  of  the  curves. 

We  now  make  some  specific  observations  about  the 
threshold  curves  for  three-junction  interferometers. 

(1)  Although  the  set  di  *di  »dj  *  ±  iw  satisfies 

Eq.  (16),  it  is  generally  inconsistent  with  Eq.  (17).  Thus  the 
maximum  Im  is  not  the  sum  of  the  separate  junction  critical 
currents,  except  for  special  cases.  Even  for  a  symmetric  in¬ 
terferometer  (/„,  ■/„  *I0,Lm  *L2,L\  *L\, 

and  L  '  *  L  J),  the  sum  of  the  individual  critical  currents 
will  be  realized  only  if 

I0}L;-2I0(L,  +L[  -Af),  (22) 

as  deduced  from  Eq.  ( 17),  and  provided  there  is  no  stray  flux. 

(2)  Consider  the  symmetric  interferometer  with  Ios 
—  2/0,  having  negligible  coupling  M  between  loops,  feed 
taps  positioned  such  that  each  current  path  contains  the 
same  amount  of  inductance  (L  "  *  L,  +L\),  and  with  no 
stray  flux.  This  interferometer  will  have  a  maximum  Lot 
4/0  when  Af2  —  .V, .  For  this  case,  one  has/9,  - 


-  Vi  * 0 j  “  where0  =  Lfo/$  and 

L*L{  +  L ;  +L  7  is  the  total  inductance  ofone  loop.  This 
0  thus  has  the  same  meaning  as  that  of  the  symmetric  two- 
junction  interferometer.  Equations  (17)  and  (19)  now  are 

2(dj  +ASj)*d)  +di 

+0>(02  +S,)  +  2v(N2  -AT,), (23) 

0  >d]  —  di  +0'(0i  ~^i) 

+  hrffft  +N2  +  2a,  Afy) ,  (24) 

where#' w>£,  + 1) , .  These  equations  show  that  flux  trapped 
in  the  feed  loop  translates  the  threshold  curve  along  the  Ic 
axis  without  distortion  in  this  case.  The  maxima  are  separat¬ 
ed  by  4ir  in  0.  For  Nf  *  0, 9  is  2tr  midway  between  the  first 
two  maxima,  which  can  happen  only  when  di  =di  4-  2tr 
according  to  Eq.  (24).  Then  Eq.  (16)  becomes 

-C,/(1+#C,),  (25) 

which  is  just  the  relation  for  a  current-symmetric  two-junc¬ 
tion  interferometer.  This  must  be  made  consistent  with  Eq. 
(23)  at  d2  *d.  +  2tr. 

d)  +  “ir  +  di  +  •  (26) 

What  is  important  to  note  about  these  two  equations,  valid  at 
the  midpoints  between  major  lobes,  is  that  the  feed  loop  in¬ 
ductance  has  dropped  out.  The  threshold  current  at  the  mid- 
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FIG.  6.  Threshold  curves  for  t  symmetric  three. junction  interferometer 
with  (ate  current  injected  close  to  the  outside  junctions  (see  Fig.  5).  L  ; 
-Zi  -0.01  pH.Zr-Zr-0.79pH.Z,  -Z,  -0 -/*  -J/0, 

—  O.S  mA.  No  stray  flu*.  snd/ty-0.  Lower  graph:  Z.  «  Z,  -  0.79  pH. 
Upper  graph  (with  ravened  threshold  current  not  shown):  Z.  —  Z,  —O' 
pH. 
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FIG.  7.  Threshold  curve  for  a  three-junction  interferometer  having  sym¬ 
metric  critical  currents,  hut  asymmetric  inductances  and  fluxes,  showing 
removal  of  lobe  cusp  ends  from  the  axis,  and  a  changing  pattern.  /„,  »/„, 
-J/n,  m  0.3  mA.  L ;  -IT-t,  «  0.4  pH;  L  i  mL~  ~Lh  -0.3  pH; 
L,  mL,  «  0;  .V,  —  Os  no  stray  flux. 


point  thus  provides  a  convenient  measure  of  the  loop  induc¬ 
tance  alone.  The  current  ratio  I^/I^  in  fact  is  precisely 
that  of  the  symmetric  two-junction  interferometer:  Not  only 
is  Eq.  (25)  the  same  as  the  symmetric  form  of  Eq.  (5),  but  Eq. 
(26)  is  just  the  same  as  Eq.  (7)  with  8  **  ir,  the  midvalue  for 
that  case.  Thus,  Fig.  4  applies  also  to  the  symmetric  three- 
junction  interferometer  as  defined  in  this  paragraph,  and,  we 
may  speculate,  would  apply  generally  to  the  ^-junction  sym¬ 
metric  interferometer  whose  internal  critical  currents  are 
twice  those  of  the  ends. 

(3)  If  the  feed  arrangement  for  supplying  gate  current  to 
the  interferometer  forms  a  superconducting  loop,  it  is  impor¬ 
tant  to  take  the  feed  loop  inductance  into  account  It  is  not 
just  the  control  current  scale  which  is  affected,  but  also  the 
shape  of  the  threshold  curve,  since  the/S's  are  affected.  Fig¬ 
ure  6  illustrates  a  hypothetical  case  in  which  the  gate  current 
is  injected  close  to  the  two  outside  junctions  by  means  of  a 
superconducting  feed  loop.  For  one  curve  the  feed-loop  in¬ 
ductance  is  made  rather  small,  and  for  the  other  it  is  taken 
equal  to  the  interferometer  inductance.  Not  only  is  the  con¬ 
trol  current  scale  much  different  for  the  two  cases,  but  the 
connectedness  of  the  threshold  curves  is  also  altered. 

(4)  In  contrast  to  the  two-junction  case,  the  threshold 
pattern  for  the  three-junction  interferometer  will  have  the 
lobe  cusp  ends  removed  from  the  Im  0  axis  when  induc¬ 
tances  are  asymmetric.  This  can  be  seen  in  Fig.  7,  in  which 
asymmetric  inductances  but  symmetric  critical  currents  are 
chosen.  Of  greater  practical  importance,  however,  is  the  fact 
that  there  is  now  a  “long-range”  change  in  the  lobe  pattern. 
This  would  ordinarily  be  undesirable  in  device  operation. 
For  example,  a  redistribution  of  trapped  flux  in  the  vicinity 
of  the  interferometer  between  runs  could  then  drastically 
alter  the  shape  of  the  threshold  current  near  zero  control 
current.  Note  that  the  maximum  threshold  current  occurs  at 
I'  -  0,  and  that  there  is  point  inversion  symmetry  about  the 
origin. 

(5)  Figure  7  does  not  reveal  an  important  point  which 
has  practical  consequences,  which  is  the  fact  that  it  is  not  the 
asymmetry  in  inductance,  per  se,  which  creates  the  long- 


range  changing  pattern,  but  rather  the  unequal  amount  of 
flux  penetrating  the  two  loops.  This  is  seen  in  simulations, 
not  shown  here,  in  which  the  threshold  curve  is  calculated 
for  unequal  inductances  but  equal  fluxes  through  the  two 
loops.  (It  would  be  unusual  of  course  to  realize  such  a  situa¬ 
tion  experimentally.)  The  result  is  an  identically  repeating 
pattern  of  the  same  period  as  in  symmetric  cases.  Thus,  care¬ 
ful  arrangement  of  control  lines  over  the  two  loops  of  the 
inu^ferometer  is  important  in  avoiding  a  lobe  pattern  which 
changes  over  many  periods. 

(6)  Asymmetry  of  critical  currents  creates  an  asymme¬ 
try  of  the  lobe  pattern  within  a  basic  period,  just  as  for  the 
two-junction  case. 


The  usefulness  of  the  foregoing  observations  is  illustrat¬ 
ed  by  considering  some  experimental  data  obtained  from  a 
three-junction  interferometer  fabricated  in  our  laboratory, 19 
and  shown  in  Fig.  8.  This  device  was  fabricated  with  a  nine- 
level  vapor  deposition  process  on  an  oxidized  silicon  sub¬ 
strate.  The  nine  levels,  in  order  of  deposition  are  as  follows: 
(1)  4000  -A  niobium  ground  plane;  (2)  1200  -A  Nb,  O, 
formed  by  anodization;  (3)  2000  -A  SiO  dielectric;  (4)  435-A 
Auln, ,  1 12  /□,  resistive  intermetallic;  (5)  3200  -A  lead-indi- 
um-gold  base  electrode  alloy;  (6)  3660  -A  SiO  to  define  the 
junction  windows;  (7)  4000  -A  lead-gold  alloy  counterelec¬ 
trode;  (8)  7200  -A  SiO  control  line  insulation;  (9)  8000  -A 
lead-indium-gold  alloy  control  lines.  Patterning  was  by  etch¬ 
ing  for  the  ground  plane  and  lift-off  for  the  remaining  eight 
levels.  The  interferometer  occupies  an  ares  of  220 X  70 ftm 
and  the  minimum  line  width  is  S  fim.  The  two  interferometer 
loops  are  folded  back  on  themselves  so  that  the  junctions  are 
clustered  in  the  center  region. 

The  inductance  of  the  various  strips  are  given  approxi¬ 
mately  by 

!»*,/(</ -M,  +  ^)/  W,  (27) 

where  /  is  the  strip  length,  d  is  the  insulation  thickness,  and 
IF  is  the  strip  width.  Equation  (27)  is  valid  for  ff\d  and  film 
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FIG.  S.  Photograph  of  the  type  o f  three-junction  interferometer  whose 
threshold  curves  are  depicted  in  Ftp  9  and  10.  Note  the  very  slight  unequal 
coupling  to  the  two  halves  of  the  interferometer  because  of  the  arrangement 
of  the  control  line 


V.  COMPARISONS  WITH  EXPERIMENT 
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thicknesses  much  greater  than  their  respective  penetration 
depths.!,  and  .  These  conditions  are  only  marginally  sat- 
isfied  for  our  interferometer.  Two  control  lines  overlie  the 
interferometer  loops.  Figure  9(a)  is  a  threshold  curve  ob¬ 
tained  using  the  outer  control  line  of  this  device.  The  essen¬ 
tial  symmetry  within  a  period  shows  that  the  critical  cur¬ 
rents  are  quite  symmetric  and  that  the  designed  1:2:1  ratio 
for  critical  currents  has  probably  been  achieved.  The  chang¬ 
ing  pattern  on  a  long-range  scale  indicates  that  unequal 
fluxes  are  penetrating  the  two  interferometer  loops.  This  is 
probably  a  result  of  the  asymmetry  in  the  path  of  the  control 
lines.  The  long-range  pattern  has  a  period  of  about  14  lobes, 
indicating  an  asymmetry  on  the  order  of  7%.  (The  long- 
range  pattern  repeats  each  time  the  difference  in  flux  pene¬ 
trating  the  two  loops  changes  by  0O-) 

Since  the  feed  loop  connects  to  the  interferometer  at  the 
midpoint  of  each  loop,  we  may  assume  that  L  "mLK  +  L ; 
as L  +  L  [,  and  Mm 0.  Thus  Eq.  (22)  is  satisfied  and 
therefore or  about  63/4  A  in  the  case  of  Fig.  9(a). 
Even  though  this  interferometer  is  slightly  asymmetric,  the 
value  of  can  be  estimated  by  measuring  the  ratio /aM //-u . 
1^  is  the  current  at  the  first  midpoint  on  either  side  of  /*»,. 
The  data  of  Fig.  9(a)  together  with  Fig.  4  yield  £  —  1.7,  and 
thus  a  single-loop  inductance  of 


L~LX  +L|  +Lr«ff$//0«S.9  pH. 

The  control  current  couples  to  L\  +L  J  and  the  parallel 
combination  of  L  7  +  L  J  and  L ,  +  £»,  as  seen  in  Eq.  (20). 
The  product  of  this  coupling  inductance  with  the  basic  peri¬ 
od  AIC  for  a  symmetric  interferometer  is  2d>„.  Taking  this 
period  from  Fig.  9(a),  one  finds 

L\  +L;L./(L7+L.)~*0/AIt-6.3  pH. 

One  fiHirin«»ai  assumption  is  necessary  to  uniquely  deter¬ 
mine  all  of  the  L  values.  Since  the  inductance  of  a  stripline  is 
proportional  to  its  length/width  ratio  [see  Eq.  (27)],  we  can 
estimate  from  the  geometry  of  Fig.  8  that  I.,  +  L\  +  L“ 
■£..  Combining  these  results,  we  obtain  the  inductances 

Lx  -  1.1,  LJ  -3.4,  L7-4.S,  L.  -8.9  pH. 

Our  control  line  geometry  suggests  that  the  7%  asym¬ 
metry  is  a  result  of  asymmetric  coupling  to  the  interferomet¬ 
er,  rather  than  unequal  inductances.  Simulation  experience 
has  shown,  however,  that  for  asymmetries  as  small  as  this,  it 
makes  very  little  difference  whether  one  chooses  in  the  mod¬ 
el  symmetric  inductances  and  unequal  coupling  [using  Eq. 
(21)],  or  simply  makes  the  inductances  asymmetric.  For  the 
theoretical  curve  of  Fig.  9(b),  we  chose  to  use  the  control 
current  form  of  Eqs.  (17H20).  and  achieved  the  coupling 

asymmetry  by  increasing  and  decreasing  the  inductances  by 
3.5%  on  left  and  right,  respectively.  The  resulting  curve  is 
quite  a  good  representation  of  the  data,  both  in  shape  and 
absolute  scale.  The  calculated  inductance  values  based  on 
Eq.  (27)  are  about  30%  higher  than  the  values  given  above. 
A  fringing-field  correction  lowers  the  calculated  value  by 
about  10%.  This  is  reasonable  agreement  in  view  of  the  un¬ 
certainty  of  the  inductance  calculation  at  the  corners,  and  in 
particular,  uncertainties  in  the  values  of  X ,,  ,  and  the  insu¬ 

lation  thickness  d  (—25%). 

Often  the  maximum  of  threshold  current  does  not  coin¬ 
cide  with  the  zero  of  control  current,  thus  indicating  the 
presence  of  flux  trapped  nearby.  In  another  experimental 
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ran,  shown  in  Ft*.  10,  the  local  maximum  increased  with 
increesing/,  near/,  —  0,  suggesting  that  trapped  flux  was 
opposing  that  from  the  control  current.  Postulating  values 
for  the  fluxes  d,i .  and  ds/  resulting  from  an  uncertain 
distribution  of  trapped  flux  is  not  very  practical  because  of 
the  great  variation  possible.  The  feed  loop,  however,  offers  a 
convenient  analytical  tool  for  including  the  effects  of 
trapped  flux.  We  found  that  the  nearby  trapped  flux  was  the 
equivalent  of  about  30  flux  quanta  in  the  feed  loop  in  this 
case.  This  number  is  interesting  because  it  is  about  the  num¬ 
ber  which  would  be  present  in  an  area  the  size  of  the  feed 
loop,  from  the  ambient  flux  pollution  field  in  the  Dewar. 
Using  a  procedure  similar  to  that  described  in  the  preceding 
paragraphs,  we  deduced  critical  currents  and  inductances 
and  achieved  a  very  satisfactory  match  to  the  experimental 
curve,  both  in  shape  and  scale. 

In  this  particular  ran,  the  curve  shifted  abruptly  from 
the  solid  curve  to  the  dashed  curve  in  Fig.  10,  coincident 
with  an  electrical  distrabtance  in  the  room.  The  reduced 

and  the  asymmetry  within  a  basic  period  show  that 
at  least  one  of  the  junction  critical  currents  has  been  reduced. 
The  decreasing  local  maxima  with  increasing  /,  near/,  »  0 
suggests  that  now  the  trapped  flux  is  penetrating  the  interfer¬ 
ometer  from  the  opposite  direction.  We  found  that  changing 
the  junction  la’i  from  the  previous  0. 10, 0.20,  and  0. 10  mA 
to  0.06, 0.1 1,  and  0. 10  mA,  and  taking  Nf  —  —  25,  along 
with  the  previous  inductances,  satisfactorily  reproduced  the 
new  curves.  Apparently  the  electrical  disturbance  had  rear¬ 
ranged  the  nearby  trapped  flux  and  had  even  caused  some  to 
lodge  in  the  junctions  themselves  even  though  the  junctions 
were  “small”  (3/jm  X  3/zm).  Thus  even  an  unusual  situation 
has  been  found  amenable  to  accurate  simulation. 

Trapped  flux  in  these  devices  can  thus  be  a  considerable 
problem.  Unless  care  is  taken  to  cool  the  device  in  a  very 
small  field,  the  shift  caused  by  trapped  flux  can  result  in 
radical  variations  in  the  threshold  curve  near  the  origin.  To 
obtain  a  number  of  trapped  flux  quanta  less  than  unity  for 
devices  such  as  that  in  Fig.  8,  one  requires  a  field  less  than 
about  10-7T(1  mG). 

It  is  sometimes  desirable  to  control  the  lobe  spacing  of 
the  threshold  curve  by  varying  the  path  of  the  control  line 
which  passes  over  the  interferometer  loops.  '•  The  lobe  spac¬ 
ing  may  be  made  arbitrarily  large  simply  by  reducing  the 
length  of  control  line  passing  over  the  interferometer  line. 
When  doing  this  it  is  important  to  note  that,  because  of  the 
shunting  effect  of  the  feed  loop,  an  interferometer  such  as 


that  in  Fig.  8  is  less  sensitive  to  control  current  across  the 
upper  half  of  its  loops  than  to  control  current  across  the 
lower  half. 

VI.  SUMMARY 

Threshold  curves  for  superconducting  interferometers 
can  be  calculated  on  the  basis  of  a  fairly  simple  theory  since 
only  the  /0  sind  term  is  needed  for  characterizing  the  Jo- 
sephson  junctions;  details  of  the  quasiparticle  curve,  resis¬ 
tances,  and  capacitance  are  immaterial  since  a  voltage  is  not 
developed.  Feed-loop  parameters  and  flux  quantization 
within  the  feed  loop  are  important  and  have  been  included  in 
■  the  present  analysis.  Experimental  curves  can  be  reproduced 
in  detail,  even  for  highly  asymmetric  cases  as  well  as  for  very 
slightly  asymmetric  cases.  The  effects  of  flux  trapping  can  be 
calculated.  A  “calibration  curve”  applicable  to  both  two- 
and  three-junction  symmetric  interferometers  provides  a 
rapid  means  for  deducing  interferometer  inductance.  An  ap¬ 
proximate  simple  analytical  expression  for  this  curve  is 
given. 

ACKNOWLEDGMENTS 

We  wish  to  thank  R.E.  Harris  for  helpful  discussions 
throughout  the  course  of  this  wor*i,  and  F.  Lloyd  and  L.O. 
Mullen  for  fabrication  of  the  experimental  devices. 


J.£  Zimmerman  ind  AH.  Silver.  Phys.  Rev.  141. 36?  (1966). 

•J.  Clarke.  Philo*.  Mag  13.  1 13  (1966V 

>A  Th  A.M.  De  Wade  and  R.  De  Broyn  Ouboter.  PhyiK*  41.  223  (1969); 
42.  626 ( 1 969);  Pivgraa  in  Law  Ttmptnmn  Physics,  edited  by  CJ.  Gorier 
(North-Holland.  Amsterdam.  1970).  Vol.  6.  p.  243. 

*T.A.  Fulton,  Solid  State  Commun.  S.  1333  (1970). 

T.A  Fulton,  L.N.  Dunkleberjer,  and  R.C.  Dyne*.  Phys.  Rev.  B  6.  833 
(1972). 

•3.  Clarke  and  3.P.  Paterson.  Appi.  Phy*.  Lett.  19, 469  (1971). 

’D.L.  Stuchm  and  C.W.  Wilmscn.  AppL  Phys.  Lett.  20. 436  (1972). 
•Won-Tien  Tseng  and  T.  Van  Duzer.  3.  Appl.  Phys.  46. 4373  (1973). 

'B.S.  Landman.  IEEE  Trans.  Magn.  MAG-13. 871  (1977). 

•H.  Beha.  Electron.  Lett- 13. 218  (1977). 

"E.O.  Schutz-DuBots  and  P.  Wolf.  Appl-  Phy*  16.  317  (1978V 
irW.H.  Henkel*.  Appl.  Phys.  Lett.  32. 829  (1978). 

"H.H.  Zappe  and  B.S.  Landman.  3.  Appl.  Phys.  49.  4149  (1978). 

'•L.M.  Oeppert.  3.H.  Greiner,  D.3.  HerrelL  and  S.  KJepner.  IEEE  Trans. 
Magn.  MAG-15. 412  (1979). 

"This  device  is  similar  to  one  described  in  Ref.  14. 

"R.E.  Harris.  C.  A  Hamilton,  and  F.L.  Lloyd.  Appi.  Phy*.  Lett.  35.  720 
(1979). 


6142 


J.  Appl.  Phys,  Vol.  50.  No.  12.  December  1*73 


Peterson  end  Hsmston 


•142 


I 

I 

I 

i 


APPENDIX  C 
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A  superconducting  loop  containing  a  Jocephaon  junction  develops  voltage  pulses  when  it 
admits  and  expels  magnetic  flux  quanta,  and  thus  may  be  used  as  an  A/D  converter. 
We  here  develop  and  discuss  several  conditions  which  must  be  satisfied  for  the 
generation  by  the  SQUID  of  an  unambiguously  countable  pulse  train,  from  which  the 
analog  signal  can  be  faithfully  reconstructed.  These  conditions  can  be  — over  « 
broad  range  of  realizable  values  of  inductance  and  resistance.  The  ice, 

however,  must  be  carefully  controlled.  The  results  of  simulations  are  also  presented, 
illustrating  the  various  ways  in  which  the  pulse  trains  are  affected. 


PACS  numbers:  74.50.  +  r 

I.  INTRODUCTION 

A  superconducting  loop  containing  a  Josephson  junc¬ 
tion — a  single-junction  SQUID  (SJ-SQU1D) — will  admit  or 
expel  magnetic  flux  in  quantized  units  with  concurrent  de¬ 
velopment  of  voltage  pulses.  These  pulses  may  be  of  picose¬ 
cond  width  at  the  junction.  Figure  1(a)  shows  one  possible 
arrangement  of  such  a  SQUID  circuit  That  this  behavior 
might  be  used  for  the  conversion  of  a  fast  analog  signal  to  a 
digital  one  has  doubtless  occurred  to  many  workers  in  the 
field.  A  patent  based  on  the  concept  has  recently  been  grant¬ 
ed,1  following  a  technical  disclosure.1  Hurrell  and  Silver1 
have  considered  a  scheme  in  which  the  pulses  from  a  SI- 
SQUID  are  counted  by  means  of  a  linear  array  of  coupled 
double-junction  SQUID’S. 

However,  there  are  many  conditions  which  must  be  met 
before  this  idea  can  be  realized.  The  purpose  of  this  paper  is 
to  examine  in  detail  those  conditions  necessary  for  the  pro¬ 
duction  of  a  train  of  well-resolved  pubes,  free  from  the  dis¬ 
tortions  which  prevent  faithful  reconstruction  of  the  analog 
signal.  We  also  present  the  results  of  some  simulations  illus¬ 
trating  the  various  effects.  Our  conclusion  is  that  such  a  de¬ 
vice  can  be  made  to  work,  in  the  sense  that  the  SQUID  pa¬ 
rameters  (RJL,C)  can  be  chosen  to  meet  all  the  conditions 
simultaneously.  A  wide  range  for  R  and  L  are  possible,  but  C 
must  be  carefully  controlled. 

Other  A/D  converter  designs  utilizing  the  Josephson 
effect  could  be  somewhat  faster.  One  can  readily  imagine  a 
network  of  sampling  circuits  operating  in  parallel,  each  cir¬ 
cuit  containing  a  Josephson  junction  sensitive  at  a  different 
threshold.  The  SQUID  A/D  converter  discussed  here  pro¬ 
duced  pulses  sequentially,  which  then  must  be  processed  se¬ 
rially.  However,  the  inherent  simplicity  and  great  linearity 
of  a  SQUID  A/D  converter  makes  it  appealing  for  use  in 
cases  where  ultimate  speed  is  not  required.  Even  so,  it  should 
be  faster  than  any  existing  A/D  converter. 

The  SQUID  model  we  shall  analyze  is  shown  in  Fig. 
1(b),  in  which  the  analog  signal  is  an  external  current  IM(t ) 
applied  directly  to  the  loop.  The  three  elements  on  the  right 
comprise  the  Josephson  junction  with  shunt  capacitance  C 


and  phase-dependent  conductance  G  (d  ).  The  superconduct¬ 
ing  phase  difference  4  across  the  junction  is  related  to  the 
junction  voltage  K(r )  by 


\2r)  dt' 


(1) 


where  «  2.068  X 10"’  V  s  is  the  flux  quantum.  The  cur¬ 
rent  I(t )  where  indicated  in  Fig.  1(b)  was  first  derived  by 
Josephson*  as 

/-/.sind  +  U(d)  V,  (2) 

where 


<r  (d  )  —  <7*  +  G|  coed-  (3) 

Here  /.  is  the  junction  critical  current,  <?.= l/R  is  the  nor¬ 
mal-state  junction  conductance,  and  G,  V  coed  is  often  re¬ 
ferred  to  as  the  coad  term,  or  the  interference  term  between 
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no.  2.  Graphics!  solution  of  Eq.  (I)  for  0  - 1  showing  hyturau  and  a  Equation  (8)  is  readily  solved  graphically,  as  shown  in  Fig,  2 

windoar  uncertainty  factor  i.  for#- 3.  Let  da  be  the  principal  value  of  cos*‘(  -  1/0),  lying 

between  Jr  and  v.  As  the  current/,  is  increased,  ^increases 
and  at  values 


the  supercurrent  and  the  quasi  particle  current  G,V.  It  is  con¬ 
ventional  in  analyzing  SQUID  behavior  to  use  Eq.  (2)  but  to 
ignore  the  cosd  term,  and  to  treat  /,  and  G,  as  constants. 

We  begin  by  considering  only  the  loop  inductance  L 
together  with  the  junction  supercurrent  since  the  basic  ideas 
of  pulse  production  and  spacing  are  contained  therein.  Re¬ 
sistance  and  capacitance  are  then  inserted  and  their  effects 
analyzed.  Finally,  we  discuss  the  effects  of  the  coad  term  and 
the  voltage  dependence  of  the  coefficients.  In  Appendix  A 
we  provide  some  of  the  details  for  estimating  the  pulse 
widths  when  inductance  is  important,  and  in  Appendix  B  we 
discuss  the  factors  affecting  the  linearity  of  a  SQUID  A/D 
converter. 

Once  a  satisfactory  sequence  of  pulses  is  produced,  the 
next  step  in  building  a  converter  is  to  use  the  pulses  to  acti¬ 
vate  a  counter  or  shift  register,  which  is  then  sampled.1  It  is 
beyond  the  scope  of  this  paper  to  discuss  these  devices, 
which  would  themselves  be  circuits  containing  Joaephson 
junctions.  Likewise,  the  various  techniques  and  algorithms 
which  might  be  used  for  reconstructing  the  signal  will  not  be 
discussrd  in  much  detail  in  this  paper,  although  we  do  give 
two  simple  examples. 

II.  CONDITIONS  FOR  ADEQUATE  PULSES 

The  differential  equation  appropriate  to  Fig.  1(b),  but 
without  the  coed  term,  is 

(4) 

where 

d-d  +0ua4.  (5) 

4,(t)-2rU,(t)/*»  («) 

0  —  lxUt/Pf  (7) 

The  essential  aspects  of  pulse  production  may  be  de¬ 
duced  by  considering  just  the  inductive,  supercurrent,  and 
driving  terms  of  Eq.  (4),  Le., 


*21«#.  +  2tn..  «-Q,l,2,~,  (9) 

jumps  discontinuously  to  a  higher  value.  When  R  or  Can 
included,  this  jump  is  not  discontinuous.  The  jumps  in  d 
correspond  to  voltage  pulses,  and  the  entry  of  flux  quanta, 
according  to  Eq.  ( 1).  When  the  signal  decreases,  the  jumps  in 
d  initiate  at 

*  -  2»-  -  d.  +  2tnr,  n  «  0, 1,2,-.  (10) 

The  voltage  pulses  are  of  reversed  polarity.  The  hysteretic 
effect  seen  in  Fig.  2  causes  an  abnormal  sparing  between 
pulses  of  opposite  polarity. 

When  0  <  1 ,  there  are  no  extrema  in  the  I,  -vs-d  plot. 
Nevertheless,  di  /dt  is  highest  when  the  slope  is  smallest  (at 
d  —  v,  3tr,— ),  and  the  voltage  will  tend  to  develop  pulses. 
These  pulses  in  fact  can  be  quite  respectable,  as  Fig.  3  shows 


TIMLm 

FIG.  4.  Top:  Paha  from  a  wry  hysttretic  SQUID,  driven  by  a  sinusoidal 
current  of  amplitude  /, »  0.2  mA  and  frequency  I  GHz.  /,  -  0.03  mA. 

Jt  ■  5  (i,  and  L  m  100  pH,  giving  0  m  13.1  Bottom:  The  drivint  current 
together  with  two  reconstructed  ugnais  baaed  on  the  counting  of  the 
SQUID  poises.  The  dashed  “curve"  it  contracted  by  addins  or  subtracting 
a  current  OJL  at  each  pulsn,  depending  on  its  polarity .  The  dotted  “curve" 
includes  hystetetic  corrections  as  discussed  in  the  ten.  The  dotted  and 
dashed  “curves”  are  drawn  so  ea  to  connect  the  leading  edges  of  each  step 
corresponding  to  a  pubs  count  The  steps  are  not  shown,  to  avoid  cluttering. 
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FIG.  5.  Top  Palm  Aram  •  SQUID  with  vvy  mall  hyiMirw.  drivon  by  i 
mmoomW  eomnt  of  MBplitado /,  ••  2.0  mA  aad  fracM.wqr  10  GHz. 

/.  a  0.1  nA.  A  —  20/2,  aad  L  «  4  pH.  L2I.  Bottom:  Tbcdnv- 

intcmrwuadwwwoaooBQCMd  agaalawithtaowo— atassaamFit. 
4. 


for  0  —  0.91.  However,  the  pulse  heights  and  widths  an 
rather  sensitive  to  the  rate  of  change  of  IM(t\ which  would 
likely  be  undesirable.  On  the  other  hand,  0>1  introduces 
the  hysteretk  effect  mentioned  above.  A  rather  extreme  case 
is  shown  in  the  top  half  of  Fig.  4,  together  with  the  badly 
distorted  reconstructed  signal  in  the  bottom  half  of  the  fig¬ 
ure.  Clearly,  large^  should  be  avoided. 

Figures  3  and  4  also  show  what  is  generally  true:  then  is 
no  cumulative  gain  or  loss  of  pulses  during  successive  sign 
changes  of  dl„  /dt. 

Pulse  spacing,  signal  reversal  effects,  and  “baseline’' 
variation  can  also  be  obtained  from  Eq.  (8),  as  follows. 

For  a  linearly  rising  singal,  /,(r)  —  I#it,  the  times  of 
occurrence  of  the  voltage  pulses  are,  from  Eqs.  (8)  and  (9), 

f.  -  +  ir-0,1,2,—,  (11) 

IvLIia 

when 

*,»*  +  <0»_l),/».  (12) 

The  separation  between  pulse  centers  is  thus 

At  -  (13) 

To  maintain  a  given  separation,  fast  signals  thus  require 
small  SQUID  ana  or  small  or  attenuated  signal  amplitude 
/,.  Pulse  separation  is  compared  to  pulse  width  below,  when 
resistance  is  introduced. 

When  the  continuous  signal  reverses,  Le.,  when  dlx  /dt 
changes  sign,  the  spacing  between  pulses  is  "abnormal”  both 
because  of  the  hysteretic  effect  mentioned  above  and  the 
“window  quantization”  effect,  common  to  all  digitization 
schemes,  in  whkh  it  is  uncertain  just  when  within  a  window 
the  turnaround  has  taken  place.  (Of  course,  since  no  signal 
turnaround  is  instantaneous,  the  pulse  spacing  will  also  in¬ 
crease  as  dl,  fit  decreases).  It  is  of  interest  to  compan  win¬ 
dow  quantization  and  SQUID  hysteresis  in  their  effects 
upon  the  spacing  between  reversed  pulses.  This  is  easily  oh* 
tamed  in  the  present  case  where  we  an  not  yet  considering  R 


and  C.  For  a  signal  rising  linearly,  stopping  suddenly,  and 
then  falling  at  the  same  rate,  one  gets  using  Eqs.  (9)  usd  ( 10) 

(14) 

when  dr  is  the  normal  pulse  spacing  given  by  Eq.  (13),  and 

d,  is  given  by  Eq.  ( 12).  The  term  6  /*  is  the  window  quantiza¬ 
tion  factor,  aad  the  remainder  is  due  to  hysteresis.  The  hys- 
tentk  factor  is  zero  at  0  —  1,  and  about  equal  toff  for 0V1. 
The  value  of  £  is  uncertain,  by  definition,  within  the  range  0 
to  2sr.  Thus  the  hysteretk  effect  dominates  the  window 
quantization  only  for  0  larter  than.  sav.  4. 

The  hysteretk  effect  on  pulse  spacing  can  be  compen¬ 
sated  in  part  during  signal  reconstruction  by  appropriately 
weighting  the  first  reversed  pulse  in  each  signal  reversal.  Fig¬ 
ure  4  shows  a  calculated  response  to  a  sinusoidal  signal,  in 
which  there  is  large  hysteresis  (0s  13).  The  lower  half  of  Fig. 
4  shows  the  associated  input  signal  together  with  two  meth¬ 
ods  for  restructuring  the  signal  from  the  pulses.  The  “naive” 
method  simply  adds  or  subtracts  a  constant  current  incre¬ 
ment  6t/L  tor  each  detected  pulse.  The  other  method  notes 
that  the  current  at  the  first  pulse  is(0»/L  XP«/2r),  according 
to  Eq.  (8),  and  that  the  current  drop  between  reversed  pulses 
i s(0«/£X&/*)—  1- Figure  3  shows  another  pulse  train  and 
reconstructed  signals  at  0—  1.  Figures  4  and  3  illustrate  the 
fidelity  with  whkh  signals  can  be  digitized  and  reproduced, 
assuming  perfect  counting.  Note  that  in  Fig.  3,  the  pulse 
number  can  be  increased  considerably  without  detriment, 

e. g.,  by  increasing/,. 

Even  for  well-separated  pubes,  the  voltage  does  not  re¬ 
turn  to  zero  between  pubes,  as  Figs.  3-3  show.  The  reason  is 
discovered  from  Eq.  (8)  by  setting  d  —  2vn  +  dp,  and  letting 
\8+ 1  <  1  in  order  to  examine  the  situation  between  pubes  (re¬ 
fer  to  Fig.  2).  By  taking  the  time  derivative  of  Eq.  (8)  and 
using  the  voltage  relation  (1),  one  gets  lor  the  voltage  be¬ 
tween  pubes 

(15) 

where 

L,«L/(l+0).  (16) 

Ly  b  the  loop  inductance  reduced  by  the  junction  parametric 
inductance  factor  1  +0,  and  appears  frequently  in  the  sub¬ 
sequent  analysis.  From  the  equations  incorporating  resis¬ 
tance  and  capacitance,  one  can  show  that  the  presence  of 
resistance  subtracts  approximately  the  term 

L]  d'IM 
R  dt  *’ 

from  the  right-hand  side  of  Eq.  (13).  Capacitance  subtracts 
another 


In  the  linearly  rising  or  falling  portions  of  the  signal,  the 
latter  two  make  no  contribution,  whereas  in  the  turn-around 
portions,  they  are  the  most  important 

The  reason  for  concern  about  baseline  variation  b  that 
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if  a  pulse  discnmiaator  has  a  constant  threshold  level  dut 
level  must  be  set  above  the  maximum  basehne  voltage;  if  this 
mtrimnn  is  high,  th— »  puke*  occurring  wheo  the  baseline 
voltage  is  low  might  be  rr:*r*‘4  However,  then  are  two  rear 
sons  for  expecting  that  baseline  varietioa  should  not  be  a 

substantial  problem.  First,  an  expression  for  K..,,  the  maxi¬ 
mum  voltage  of  a  pulse,  may  be  obtained  from  VamtTf  ss$0, 
when  /is  a  measure  of  the  pulse  width.  For  wett  resolved 
pulses,  r'SAt  /M,  when  M  is,  say,  4  or  5-  Noting  the  rela¬ 
tionship  between  4/ and  Kfl  given  by  Eqs.  (13)  and  (15),  one 
thus  gets 

+ff)F,  (17) 

Thus  for  well-separated  pubes,  the  pulse  height  will  be  sever¬ 
al  tiiM  greater  than  the  baseline  voltage  during  the  linear 
portions  of  the  signaL  Similar  arguments  can  be  used  for  the 
changing  parts  of  die  signal,  with  similar  conclusions.  This, 
of  course,  does  not  guarantee  that  a  given  pulse  height  will  be 
greater  than  the  maximum  baseline  voltage.  Of  mom  impor¬ 
tance  is  the  fact  that  it  should  be  possible  to  design  the  dis¬ 
criminator  level  to  vary  according  to  Fa(r),  or  to  filter  out 
that  background  on  the  basis  of  frequency  in  the  circuitry 
following  the  SQUID. 

A  junction  with  finite  resistance  R  is  known1  to  switch 
in  and  out  of  its  voltage  state,  when  driven  with  a  current 
larger  than  /« in  the  time 

r,  -  <t>JirRK  <l») 

This  is  the  voltage  pulse  width  at  half-height  We  wish  to 
determine  now  by  how  much  this  pulse  width  is  affected  by 
the  presence  of  the  SQUID  inductance  L.  By  consulting  Fig. 
2  and  Eq.  (4)  without  the  capacitance  term,  we  see  that  the 
state  of  the  system,  represented  by  a  “particle"  riding  on  the 
9  curve,  does  now  not  simply  follow  /,(»),  but  lags  by  an 
amount  (L  /R  YMd  /dt ).  Thus,  discontinuous  jumps  in  d  dc 
not  occur;  rather  the  particle  makes  a  fast  ride  through  the 
local  valleys.  This  “inertial"  effect  can  also  be  manifested  in 
the  appearance  of  a  final  positive  pulse  just  after  an  increas¬ 
ing  signal  has  started  back  down. 

Good  pulse  resolution  requires  that  the  pulse  width  r  be 
much  less  than  pulse  spacing  At.  The  finite  resistance  now 
being  considered  still  not  alter  pulse  spacing;  i.e.,  when  the 
transients  are  finished,  the  “steady-state"  response  to  a  lin¬ 
early  increasing  signal  must  be  a  series  of  pulses  with  the 
same  periodicity  as  when  the  resistance  is  infinite.  Thus  Eq. 
(13)  for  the  spacing  At  is  still  valid. 

The  effect  of  loop  inductance  on  pulse  width  is  not  sim¬ 
ple  to  analyze.  An  approximation  is  developed  in  Appendix 
A.  The  expressions  obtained  there  do  not  yield  a  simple  for¬ 
mula  for  r  except  at ^>1.  However,  whenff  1.3,  r  and  r, 
are  not  so  different  that  we  cannot  replace  r  by  r,  in  the  puke 
resolution  condition.  The  simulations  show  that  At  —  4r,  is 
about  the  minimum  condition  for  adequate  resolution.  Thus 
from  Eqs.  (13)  and  (18),  we  arrive  at 

L/ffS70/V7,  (19) 

as  a  reasonable  condition  for  satisfactory  pulse  resolution, 
where/—  a/1  w. 
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Figures  3-3  arc  representative  of  simalaboasoffrQUID 
behavior  whaa  junction  (or  sham)  ramataaos iffiamdudad, 

developed  pulaaa  from  a  aoohystaretk  SQUID  (ffasO.9). 

tioa  (19)  is  well  satisfisd,  aad  all  puleas  are  wall  above  the 
maximum  value  of  Vt .  Dm  measured  1/e  width  ofMO- 

0.85  ps  for  the  first  pulse  (Le.,  in  the  region  of  a  linearly  rising 
signal)  is  more  than  a  factor  of  2  greater  than  the  siaapie 

junction  r„  but  doae  to  the  r  as  calculated  from  Eqs.  (Al)~ 

(A4).  The  third  aad  fourth  pulses,  occurring  when  the  signal 

is  changing  more  slowly,  arc  noticeably  broader,  as  expected 
for  a  aonhysteretk  SQUID. 

In  the  very  hysteretic  SQUID  13)  of  Fig.  4.  the 
pulses  are  well  developed  and  Eq.  ( 19)  is  just  satisfied  The 
long  delay  between  pulses  of  opposite  polarity  shows  princi¬ 
pally  the  effects  of  the  great  hysteresis.  Pulse  heights  here  are 
far  above  Vt ,  and  pulse  broadening  as  the  signal  stows  is  not 
discernible.  Now  r,  and  r  are  about  equal,  and  close  to  the 
measured  width. 

Figure  3  shows  the  response  of  a  slightly  hysteretic 
SQUID  (0s  1.2).  Only  the  inductance  is  dumped  from  Fig. 
3.  The  inductance  increase  by  one-third  results  in  a  pulse 
width  decrease  by  one-third.  A  slight  broadening  of  the 
pulses  in  the  signal  slow-down  regions  can  be  seen. 

Good  pulse  trains  of  course  also  occur  in  quite  different 
regimes  of  operation.  For  example,  simulations  for  a  rather 
typical  shunted  point-contact  SQUID  (L  —  500  pH,  R 
m  100/4/2, 1,  m  1  fiA)  driven  at  1  kHz  show  good  pulses. 
However,  the  very  low  peak  voltages  ( &  10-**  V)  for  such  a 
case  would  make  pulse  detection  perhaps  impossibly  diffi¬ 
cult  because  of  noise. 

Capacitance  is  the  bad  actor  in  pulse  production  from 
SQUID’S  and  must  be  kept  as  small  as  possible,  consistent 
with  an  adequate  current  /«.  The  main  effects  of  capacitance 
can  be  understood  qualitatively  by  the  device  of  rewriting 
Eq.  (4)  in  the  form  of  an  energy  equation,  in  which  a  particle 
can  be  pictured  as  sliding  along  a  time-dependent  potential 
well,  and  restricted  by  viscous  damping.  Thus.  Eq.  (4)  is 
equivalent  to 

«  -/*’-(*-*>„  (20) 

where  the  dot  over  a  symbol  refers  to  differentiation  with 

respect  to  t,=r /(LC)xn.  Here D  -  (LC)'n/RC, 

and 

rmW-4My-0  coed-  (21) 

The  “potential”  V  is  a  parabola  centered  at  4M,  with 
local  minima  which  become  progressively  shallower,  and  fi¬ 
nally  disappear,  as  one  moves  farther  up  the  sides  of  the 
parabola.1  As/,  is  increased,  a  particle  once  displaced  from 
its  vanished  minimum  drops  into  the  adjacent  well,  and  if  the 
damping  constant  D  is  large  enough,  becomes  trapped  in 
that  well,  and  oscillates  about  the  minimum  before  coming 
to  rest  Each  new  local  minimum  entered  corresponds  to  the 
generation  of  a  voltage  pulse,  and  the  oscillations  within  a 
local  minimum  are  the  plasma  oscillations.  In  the  present 
case,  /,  does  not  generally  change  slowly,  and  the  pulse 
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no.  7.  A 4-pH  SQUID  Smaa  at  10 OH*,  with/. « ai  nA. /,  •  2.0mA. 
tad*  >20/2. («)C>0t(b) C> 0.03 pF  Tht  plwau  aaaUsaoai  a  (fee 
tanar  can  prodaeas  aa  uaaccaptaMt  paka  train,  aad  Eq.  (23)  a  not 


no.  A  Aa  10-fH  SQUID  dri*M«2.S  OH*,  vith  /. »  0. 1  oA.  /,  >  0.3 
■A.  aad  K  m  3  a.  (*)  C  -  ft  (b)  C  -  02  pF;  (c)  C  -  OS  pF.  Tbe  «pw 

iUaatratan  tbaappaanaeaafaa  *m  paiaa(c)  whaa  the  capacKaaoa  is  large 
aaoagh  that  tin  daaipiag  it  I—  da*  thrahold  damping  tar  a*gk  Bujund 
aatry.  Eqaanoa  (23)  it  aot  mMid  a  (ha  saddle  tad  bowon  cuna. 


characteristics  an  determined  also  by  the  rate  of  change  of 
I,  as  well  as  by  the  value*  of  D  and  #,  as  Eq.  (20)  shows. 

If  the  damping  D  a  too  small  (capacitance  too  large),  at 
least  three  undesirable  things  can  happen.  One,  the  plasma 
oscillations  may  be  of  such  large  amplitude  that  they  could 
be  counted  as  pulses.  Second,  if  the  ringing  does  not  damp 
out  soon  enough,  it  can  interfere  with  the  effect  of  the  chang¬ 
ing  potential,  causing  the  particle  to  move  too  soon  or  too 
late  into  the  next  local  minimum.  Third,  even  for  a  very 
slowly  changing  signal,  if  D  is  smaller  than  a  threshold  value 
A,  the  particle  might  move  rapidly  through  more  than  one 
minimum,  producing  pulses  in  groups.  If  the  signal  at  a  given 
time  is  fast  enough  to  prevent  this,  extra  pulses  can  still  ap¬ 
pear  as  the  signal  slows.  All  these  effects  have  been  seen  in 
the  simulations.  Obviously,  any  of  them  could  produce  a 
digitized  signal  so  distorted  that  the  result  may  be  worthless. 

The  damping  threshold  A  mentioned  above  is  the  con¬ 
dition  for  single-pulse  (as  opposed  to  multiple)  generation  as 
the  signal  is  very  slowly  increased  through  a  value  which 
gives  a  pulse.  The  formula 

A- 1-1930 1/1 -2.530  (22) 

has  been  shown*  to  be  quite  accurate  for#  above  about  4.6 
and  is  certainly  accurate  enough  for  SQUID  A/D  conver¬ 
sion  considerations.  When  the  signal  changes  rapidly,  the 
threshold  between  single  and  multiple  quantum  transitions 
becomes  blurred,  however. 

Figures  6(*)-6(c)  illustrate  an  effect  when  D  <  A-  Here 
an  SOpH  SQUID  is  driven  at  2.5  GHz.  Figure  6(a)  for  C  ■  0 
shows  an  acceptable  pulse  train  even  though  pulse  spacing  is 


perhaps  not  large  enough.  At  C  —  0.2  pF  in  Fig.  6(b),  the 
pulses  are  still  regularly  spaced,  although  resolution  may 
now  be  s  problem.  Note  that  the  sixteenth  pulse  in  Figs.  6(a) 
and  6(b)  actually  occurs  after  the  signal  has  started  to  de¬ 
crease.  This  is  common  in  the  simulations  and  illustrates  the 
“inertial’’  effect  mentioned  earlier.  The  damping  D  is  4.0  in 
Fig.  6(b),  larger  than  the  threshold  A  of  3.3  Figure  6(c)  at 
C  ■  0.5  pF  shows  a  still  worse  train,  but  of  interest  here  is 
the  fact  that  there  are  now  17  pulses  (and  the  seventeenth 
occurs  before  signal  reversal).  This  is  because  now  D  «■  2.5, 
less  than  A- 

Figures  7(a)  and  7(b)  illustrate  plasma  oscillation  ef¬ 
fects  in  a  4-pH  SQUID  driven  at  10  GHz.  Figure  7(a)  shows 
a  clean  train  of  pulses  at  C  «  0.  However,  when  C  is  in¬ 
creased  to  just  0.03  pF  in  Fig.  7(b).  large-amplitude  plasma 
oscillations  appear.  Note  that  the  pulse  spacing  is  relatively 
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unchanged,  but  that  there  would  now  be  a  problem  of  bow  to 
avoid  counting  the  larger  of  the  piaama  oscillations  as  pulses. 

Figures  8(a)  and  8(b)  illustrate  the  “ringing  interfer¬ 
ence"  mentioned  above,  with  a  10- pF  SQUID  responding  to 
a  10-GHz  signal  At  C  «  0,  Fig.  8(a),  a  good  pulse  train 
appears,  but  at  C  »  0.0S  pF,  Fig.  8(b),  an  intolerable  situa¬ 
tion  has  already  developed.  The  pulses  are  now  irregularly 
spaced  because  of  the  ringing  interference  between  the  plas¬ 
ma  oscillations  and  the  signal  motion. 

Figures  9(a)-9(c)  show  a  case  in  which  a  5-pH  SQUID 
driven  at  10  GHz  haa  a  lair  pulse  train  at  C  •  0,  Fig.  9(a), 
which  is  not  materially  affected  when  a  capacitance  of  0. 1  pF 
is  introduced.  Fig.  9(b).  Figure  9(c)  is  addressed  later  in  con¬ 
nection  with  the  coed  term. 


The  conditions  for  maintaining  a  good  pulse  train  in  the 
presence  of  the  inevitable  capacitance  are  deduced  as  fol¬ 
lows.  Of  principal  interest  are  the  plasma  oscillation  period 
and  damping  time.  Thcae  may  be  analyzed  by  examining  Eq. 
(4)  near  the  midpoints  between  the  pulses,  at  d  -  2en  -+  54, 
where  | M  |  <  1-  The  homogeneous  solution  to  Eq.  (4)  linear- 
ixad  in  64  has  the  form  etpfpt),  where 


1/1 


and  the  redacad  inductance  L,  is 


(23) 

in  Eq.  (Id).  Whan 


having  the  period 

Tp  m  WL,C (4£,C  —  L  */R  r)  '"  (24) 

and  damping  time  2AC. 

The  various  sihets  seen  ia  Figs.  7-9  can  be  pfidiettd 
from  knowledge  of  the  SQUID  and  junction  parameters. 
First,  let  as  dispose  of  two  insufficient  conditions  for  good 
pulse  trains.  If  the  pulse  spacing  dr  satisfied  dt»2AC,  there 


would  be  no  ringing  interference,  and  the  pulses  would 
maintain  their  correct  spacing.  However,  this  condition  is 
too  weak  since  it  still  allows  the  possibility  that  initial  large- 
amplitude  plasma  osdUatioas  may  be  counted  as  pulses. 
This  could  be  avoided  by  requiring  the  plasma  period  to  be 
long  compared  with  At,  so  that  succeeding  pulses  occur  be¬ 
fore  the  plasma  oscillations  have  a  chance  to  develop.  But  it 
is  clear  that  this  condition  is  also  not  sufficient;  the  signal 
must  eventually  slow  and  reverse,  and  the  resulting  in¬ 
creased  pulse  spacing  would  allow  the  plasma  oscillation  ef¬ 
fects  to  occur  again. 

One  way  to  avoid  these  problems  is  to  require  the  “over- 
damped"  condition,  in  which  the  radicand  of  Eq.  (23)  is  posi¬ 
tive:  L,/R  >  ARC.  Small  R  and  Care  thus  desirable.  Howev¬ 
er,  modest  underdamping  can  certainly  be  tolerated,  as  our 
simulations  generally  show.  This  may  be  quantified  as  fol¬ 
lows.  Noting  that  the  failing  side  of  a  pulse  is  actually  the 
first  psrt  of  the  plasma  oscillation,  we  should  require  that 
{T, — the  approximate  time  at  which  the  voltage  has  its  least 
value— should  be  somewhat  greater  than  the  damping  time 
IRC,  in  order  that  the  oscillation  amplitudes  not  be  so  large 
as  to  be  possibly  counted  as  pulses.  Using  Tf  ia(2RC) 
where  a  is  somewhat  larger  than  2  for  a  minimal  condition, 
we  have  with  Eq.  (24) 


R  >  4a*  +a* 


-AC. 


(25) 


Simulations  show  that  a  ->  3  is  about  right.  The  coefficient 
of  AC  in  Eq.  (25)  is  then  about  J. 


It  may  be  necessary  to  use  a  shunt  resistance  to  satisfy 
this  condition.  Let  us  consider  an  example  which  would  ap¬ 
pear  to  be  fairly  typical  Take  C  -  0.3  pF,  an  attainable  val¬ 
ue.  For/,  —  0.2  mAandd  *  2  mV,  where  d  is  the  half-gap 
voltage,  the  normal-state  junction  resistance  A  j  would  be 
about  16  A,  according  to  Rjl,  —  Jwd.  If  the  pulse  voltage 
does  not  rise  above  about  2  mV,  the  actual  junction  resis¬ 
tance  will  be  considerably  greater  than  16  (1.  Thus  RCt  5 
pa,  requiring  L,  >  60  pH  to  satisfy  Eq.  (25).  But  taking  note 
of  the  definition  of  0,  we  see  that  L,  >  60  pH  cannot  be  satis¬ 
fied  for  positive  L.  In  feet,  Eq.  (25)  demands  a  maximum 
value  for  A  of  about 


A^  -  Q0SirUCV  (26) 

Figures  6-8  illustrate  cases  in  which  Eq.  (25)  is  not  sat¬ 
isfied,  whereas  in  Fig.  9  it  is  satisfied. 

To  sum  up,  Eqs.  ( 19)  and  (25)  should  be  simultaneously 
satisfied  for  good  pulse  production  from  a  SJ-SQUID. 


III.  OTHER  EFFECTS 

Equation  (2)  with  conductance  C(d)  given  by  Eq.  (3)  is 
derived*  for  constant  voltage  acroas  the  junction.  So  fer  we 
have  ignored  the  coed  term  and  have  treated  G,  —  1/A  and/, 
as  constants,  even  though  the  derivation  shows  them  to  be 
voltage  dependent  When  one  confronu  a  situation,  as  here, 
in  which  picosecond  voltage  pulses  occur,  whose  peak  values 
can  reach  or  exceed  the  gap  voltage,  several  questions  arise. 
Will  a  nonequilibrium  description  of  the  superconductors  be 
mstntial  for  a  quantitative  analysis?  Should  a  tune-depen- 
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dent  microscopic  formulation’4  (but  still  within  an  equilibri- 
um  theory)  be  used?  Should  the  voltage  dependence  of  the 
parameters  /«,  Gu  and  G,  be  used  in  the  constant- vohage 
formulation?  How  important  ia  the  coed  term?  We  thall  ad- 
dreaa  questions  in  the  order  asked. 

The  nonequilibrium  questioo  ia  difficult  to  aoawer  rig¬ 
orously.  As  long  as  the  superconductors  are  not  driven  far 
out  of  equilibrium  (say  with  voltages  rirmring  the  gap  volt¬ 
age),  or  if  the  equilibration  tunes  are  shorter  than  the  pulse 
widths  when  the  voltage  peaks  an  substantial,  one  would 
not  expect  a  grow  deviation  from  equilibrium-based  compu¬ 
tations.  Since  superconductor  relaxation  times  are  substan¬ 
tially  lew  than  picosecond  rimes,  and  since  the  pulse  peaks 
only  rarely  exceeded  gap  values  in  our  simulations,  we  ex¬ 
pect  that  nonequilibrium  effects  should  not  be  serious  in 
SQUID  A/D  conversion. 

As  for  going  to  a  time-dependent  formulation,  this  is  in 
fact  what  one  should  do  for  the  best  analysis  within  an  equi¬ 
librium  theory.  Such  an  analysis,  is.  of  course,  much  more 
expensive  than  the  present  one.  We  have  already  made  com¬ 
parisons  between  the  microscopic’  and  phenomenological’ 
theories  in  the  context  of  examining  the  possibility  of  using 
tunnel  junctions  as  picosecond  pubers  and  found  no  sub¬ 
stantial  differences  in  pulse  configurations.  Further,  it  has 
been  observed  for  many  years  that  the  phenomenological,  or 
resmivdy-shunted-junctioa.  model  is  respectable  in  a  wide 
variety  of  applications.  The  reason  is  probably  somewhat  w 
follows,  and  this  also  addresses  the  remaining  questions 
asked  above.  The  critical  current  /.  vanes  only  slightly  with 
V  for  K< d  and  thus  can  be  treated  w  a  constant  in  A/D 
converter  applications.  A  shunt  conductance,  which  would 
probably  be  mad.  will  usually  dominate  G+  Since  G,  and  G, 
are  comparable  in  magnitude,  the  voltage  dependence  of 
both  G,  and  (7,  is  not  then  significant.  The  time-dependent 
formulation  thus  probably  introduces  changes  only  in  detail 

Some  effects  of  the  coed  term,  and  a  new  way  for  mea¬ 
suring  G,  have  been  discussed  in  Ref.  6.  It  is,  however,  inter¬ 
esting  to  note  the  effect  of  the  coed  term  in  the  praecat  con¬ 
text  One  sew  that  its  contribution  to  Eq.  (20)  is  to  multiply 
the  damping  coefficient  D  by  the  factor  1  +  r  coed,  where 
*  —  GJGf  BCS  theory  applied  to  tunnel  junctions  states 
thatrs  +  1  for  V<A.  Thus  at  d  —  2m,  *  -  0,1,2-,  which 
are  the  positions  of  the  local  minima  in  the  ‘’potential’*  f, 
the  effective  damping  is  doubled;  the  bottoms  of  the  weds  are 
“stickier"  and  plasma  oscillations  would  be  lew  of  a  prob¬ 
lem;  that  is,  larger  capacitance  could  be  tolerated.  However, 
ifrs  -  1,  wit  is  in  the  theory  of  microbridges,  and  w  sever¬ 
al  eaperiwanlt  on  microbridges,  point  contacts,  and  tunnel 
junctions  have  suggested  (sec  Ref.  6  for  a  fuller  diacuasion 
and  references),  then  the  damping  in  the  local  walla  would  be 
very  much  smaller,  and  plasma  oactllarions  would  pom  a 
more  seven  problem.  We  show  an  example  of  this  ia  Fig. 
9(c),  in  which  the  parameters  are  the  same  w  in  the  middle 
curve,  except  that  the  coed  term  with  t  ■  -  1  is  added.  Al¬ 
though  Eq.  (23)  is  satisfied,  this  pubs  train  would  probably 
not  be  acceptable.  We  repeat  the  caveat  given  above,  that  in 
the  realm  of  rapid  and  hugs  voitags  changes,  the  oonstant- 
vohage  formulation  may  not  be  quantitatively  rahabla, 


However,  sve  would  abo  expect  it  to  be  not  completely 
wrong.  Thus,  if  it  turns  out  that  ts  —  1  for  tunnel  junctions 
at  V<A,  the  avoidance  of  disruptive  plasma  oscillations  may 
be  very  difficult  unlew  a  large  shunt  conductance  b 
provided. 
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APPENDIX  A 

In  this  Appendix,  we  analyze  the  voltage  pulse  width 
for  a  junction  with  resistance  but  no  capacitance,  taking  into 
account  the  loop  inductance. 

We  consider  a  linearly  rising  signal  Ix  m  before. 
From  Fig.  2  one  can  see  that  the  pulses  will  be  asymmetric. 
The  voltage  pulse  “begins”  when  8,  defined  in  Eq.  (3),  ceases 
to  rise  linearly,  and  rises  most  rapidly  when  8  decreases  moat 
rapidly,  at  about  8  m  ir,3r,~.  The  pulse  peaks  near  the  mini¬ 
mum  of  8  then  generally  falb  more  rapidly  than  it  rose.  We 
can  estimate  the  pube  width  by  replacing  8  by  a  series  of 
straight  lines  on  a  0-vs-d  plot,  with  the  lines  passing  through 
the  points  of  greatest  slope,  and  having  the  correct  slopes 
(I  ±  /9 )  at  those  points.  The  solutions  to  the  approximated 
Eq.  (4)  with  C  —  0  axe  then  increasing  or  decreasing  expo¬ 
nential  functions.  The  decreasing  part  of  the  pube  is  expo¬ 
nential  with  rime  constant 

r,  •  UR,  (Al) 

where  L,  a  the  reduced  inductance  of  Eq.  (16).  The  increas¬ 
ing  exponential  has  a  time  constant  (L/R  X/?-  1),  but  this 
dow  not  measure  the  “left”  half  pube  width.  If  we  equate  the 
voltage  on  the  rising  portion  to  the  l/<  point  on  the  falling 
side,  we  find 


L  1 
R  0-1 
where  X  and  Y  are  solutions  of 


(A2) 

(A3) 


,rm  1 +*"(** -I).  (A4) 

Many  simulations  have  shown  that  the  resulting  com¬ 
puted  pube  width  r  m  r,  +  r,  b  always  a  littb  smaller  than 
the  measured  width  at  the  l/e  level,  averaging  about  10%  of 
the  measured  value.  Ia  fact,  the  maple  expression  for  r,  ia 
Eq.  (lff)b  found  to  be  good  to  withm  about  30%  for  fft  1.3. 
Forff  S  1.3,  we  have  found  tome  caarnin  which  r,b  ia  error 
by  w  much  as  a  factor  of  3.  However,  r  remamt  “good"  at  all 
0 

The  pube  width  r  according  to  the  above  formulas  is 
thus  dependent  upon  the  signal  rate,  especially  nearff  •  1. 
as  oat  expects.  Whanffbl,  X  -  Fsl.  and 

t~1L  /0R  -  v  (A3) 
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FIG.  10.  Poroom  of  (wo  pouibfe  junction  amnfemcnts  in  SQUID's  con¬ 
structed  by  mkrofabncsikm  techniques. 


Thus,  the  very  hysteretic  SQUID  develops  voltage  pulses 
comparable  to  those  produced  by  the  simple  resistive  junc¬ 
tion.’  whereas  the  nonhystemic  or  only  slightly  hysteretic 
SQUID  gives  rate-dependent  pulses. 

APPENDIX  B 

Here  we  shall  examine  some  factors  affecting  the  linear¬ 
ity  of  the  SQUID.  As  was  seen  in  Sec.  II.  the  two  most  im¬ 
portant  SQUID  parameters  are  the  inductance  L  of  the  loop 
and  the  critical  current  /,  of  the  junction.  Both  of  these  pa¬ 
rameters  can  be  changed  by  the  signal  being  measured,  pro¬ 
ducing  a  nonlinearity. 

Typically,  the  magnetic  field  from  the  current  in  a  near¬ 
by  control  line  will  be  used  to  bring  the  critical  current  to 
some  convenient  operating  point.  The  addition  of  a  signal, 
whether  as  a  magnetic  field  or  as  a  current,  will  change  the 
number  of  (luxoids  in  the  loop,  which  in  turn  will  change  /« if 
the  associated  magnetic  induction  B  crosses  the  junction  (see 
Fig.  10).  We  need  to  know  by  how  much  I,  changes  for  a 
given  number  of  fluxoids  entering  the  loop  due  to  an  increas¬ 
ing  signal 

Consider  a  square  tunnel  junction  of  side  W  and  thick- 
ness d-  If  WSXJt  where  X,  is  the  Joaephaon  penetration 
depth,  the  critical  current  follows  the  Fraunhofer  diffraction 
formula1* 


/,  -  /„  |(suu)/x|. 

(Bl) 

when 

x~BW(2X  +dytr/+. 

(B2) 

and  X  is  the  Loodon  penetration  depth  of  the  magnetic  field 
into  a  superconductor.  A  skewed  /,-ve- B  curve,  rather  than 
Eq.  (Bl),  results  for  configurations  utilizing  a  ground  plane, 
but  Eq.  (Bl)  is  adequate  for  our  purpose.  Expanding  suu 
about  a  representative  operating  point  x%  -  2  rad,  and  re¬ 
placing  B  by  N+Jtt*,  where  r  is  an  effective  radius  of  the 


loop  and  ATis  the  number  of  fiux  quanta  enclosed  by  the  loop, 
one  obtains 


N-N, op- 


(B3) 


Here  we  have  dropped  4530  A  compared  to  A  £300  A,  and 
set  Afop  equal  to  the  number  of  enclosed  flux  quanta  at  the 
operating  point  Clearly,  the  number  of  enclosed  fiux  quanta 
is  quite  sensitive  to  the  loop  and  junction  dimensions.  For 
X  a  860  A  (niobium),  a  junction  with  IF*  0.1  mm  and  a 
loop  of  r  *  1  mm— dimensions  which  are  large  by  typical 
miciofabrication  standards — the  relative  suppression  of  /« is 
1.7x  10?*  for  each  new  fiux  quantum,  or  about  3X 10"*  for  a 
signal  producing  16  pulses  (4-bit  resolution).  That  is,  each 
pulse  will  be  mispositioned  relative  to  its  neighbor  by  a  rela¬ 
tive  amount  1.7  X 10"’,  and  a  train  of  16  pulses  of  the  same 
polarity,  from  a  linearly  increasing  positive  signal,  will  be 
compressed  by  3  parts  in  10*.  Stated  another  way,  if  a  nomi¬ 
nal  16-pulse  signal  were  within  3  parts  in  10*  of  its  17-pulse 
value,  it  would  generate  that  seventeenth  pulse  because  of 
the  decrease  of 


The  error  increases  as  the  r  and  W  dimensions  are  re¬ 
duced  to  the  microfabrication  scale.  However,  the  nonlin¬ 
earity  is  considerably  reduced  if  one  would  operate  near  the 
top  of  the  7,-vs-ff  curve,  where  the  slope  is  small.  For  exam¬ 
ple,  at  x<  1,  Eq.  (Bl)  gives  a  relative  suppression  of  /,  equal 
'  to  3  X 10'"  per  quantum,  for  the  dimensions  given  above. 


All  the  above  numbers  are  not,  possibly,  as  good  as  may 
be  desired,  but  it  should  be  pointed  out  that  we  have  treated 
the  wont  case,  in  which  all  of  the  flux  density  within  the  loop 
is  assumed  to  cross  the  junction,  for  example  as  depicted  in 
Fig.  10(a).  An  overlap  geometry  as  shown  in  Fig.  10(b) 
would  largely  shield  the  junction  from  the  magnetic  field. 
Also,  one  could  design  shielding  elements  about  the  junc¬ 
tion.  One  may  wish  to  use  the  SQUID  as  an  essentially  dc 
ammeter  with  extreme  precision,  in  which  case  a  very  large 
number  of  quanta  would  be  desired.  For  such  applications, 
one  could  scale  up  the  effective  loop  radius,  as  well  as  pro¬ 
vide  for  shielding.  Thus  it  appears  that  the  nonlinearity  asso¬ 
ciated  with  the  field  dependence  of  the  critical  current  is  a 
problem  which  can  be  handled  readily. 


A  second  source  of  nonlinearity  is  the  field  dependence 
of  the  penetration  depth.  Ginzburg-Landau  theory 


K(K  +  V 8) 

8(*  +  V2)  S') 


(B4) 


For  niobium  at  4  K,  the  Ginzburg-Landau  parameter  k  is  1, 
and  the  critical  field  Bc  is  0. 16  T.  A  field  of  3.6X10'*  Twill 
change /i  by  1  part  in  10*.  It  would  change  the  effective  area 
and  the  number  of  enclosed  flux  quanta  by  twice  that 
amount  This  field  corresponds  to  more  than  5  x  10*  flux 
quanta  for  a  loop  of  1  mm  effective  radius,  and  3  X 10*  quanta 
for  a  100-pm  radius.  Thus,  for  the  usual  fields  expected  for 
A/D  conversion,  the  loop  area  and  inductance  would 
change  by  entirely  negligible  amounts.  For  s  precision  dc 
ammeter,  the  effect  could  be  made  negligible  by  scaling  up 
the  dimensions. 
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ATTENUATION  IN  SUPERCONDUCTING  STRIPLINES + 


Richard  L.  Kautz 


ABSTRACT 


Huturnuti  of  che  Q  of  stripline  resooacora 
yield  value*  for  the  attenuation  of  Nb-Nb.O.  -Pb 
ecrlplinae  typical  of  choee  ueed  in  superconducting 
microcircuits.  Ac  4  K  the  attenuation  between  50  and 
500  MBs  la  proportional  to  frequency  and  probably 
raaulcs  froa  dielectric  loaaee.  Near  che  tranaition 
ceaperature  of  Pb,  Che  actenuation  begin*  Co  show  the 
frequency-squared  dependence  associated  with  super¬ 
conducting  loeses. 


I.  INTRODUCTION 


Superconducting  scripline  resonators  have  been 
used  by  several  authors  tu  examine  the  surface 
impedance  Z  of  superconductors.  Measurement  of  the 
resonant  frequency  has  been  used  to  determine  che 
imaginary  part  of  Z.  and  hence  che  superconducting 
penetration  depch  ^  4  and  measurement  of  che  Q  has 
bean  used  to  determine  che  real  part  of  Z  The 

prasenc  study  applies  the  resonant-cavity5 technique 
to  measure  che  attenuation  of  a  Nb-Nb.O.-Pb  scrip- 
line. 

II.  SAMPLE  PREPARATION 


The  striplines  to  be  considered  consist  of  a 
niobium  groundplana,  a  dielectric  layer  of  anodic 
Nb.Oc,  and  a  Pb-Au  alloy  scrip.  In  fabrication,  a 
400  nm  layer  of  Nb  is  first  deposited  on  a  silicon 
wafer  using  e-beaa  evaporation  in  ultra-high  vacuua. 

The  Nb  is  Chen  anodized  under  conditions  similar  to 
chose  of  Hlckmoct  and  Hiatt®.  The  electrolyte  con¬ 
sists  of  18  gm  of  ammonium  pencaborace  dissolved  in 
100  m)  of  ethylene  glycol.  Anodization  is  carried 
out  at  70°  C  at  a  current  density  of  1  mA/cm*  until 
che  cell  voltage  reaches  50  V.  This  voltage  is  then  „ 
maintained  uncil  the  current  has  fallen  to  0.05  mA/ca~. 
Assuming  a  dielectric  constant4*  7  of  29,  capacitance 
measurements  yield  a  dielectric  thickness  of  104  nm. 

The  striplinc  is  completed  by  evaporation  of  Pb  and 
Au  Co  form  a  pair  of  Pb-Au  alloy  scrips  patterned  by 
photoresist  lift-off.  Each  scrip  is  600  nm  thick, 

50  urn  wide,  and  70  cm  long.  Use  of  chis  balanced 
pair  of  striplines,  shown  in  Fig.  1,  avoids  the 
necessity  of  contacting  the  groundplane.  Since  the 
75  pa  separation  between  strips  is  much  greater  than 
the  104  nm  dielectric  thickness,  che  two  strips  are 
not  coupled.  • 

III.  EXPERIMENT 

The  experimental  arrangement  for  studying  the 
stripline  samples  consists  of  a  sweep  oscillator 
connected  to  one  end  of  the  scripline  via  a  50  li 
coaxial  line  and  an  oscilloscope  connected  to  che 
other  end  also  via  50  u  coax.  The  difference  in 
characteristic  impedance  between  the  scrlpllne  (a 
few  tenths  of  an  ohm)  and  che  coaxial  lines  pro¬ 
duces  a  large  reflection  coefficient  at  che  ends 
of  the  scrlpllne  and  turns  it  lnco  a  resonant 


Manuscript  received  September  28,  1978. 

''Partially  supported  by  the  Office  of  Naval  Research 
under  contract  number  X00014-77-F-0048,  dated  May  1, 
1977.  Contribution  of  the  U.  S.  Government,  not 
subject  to  copyright. 

•NRC/NBS  Postdoctoral  Fellow,  Electromagnetic 
Technology  Division,  Nacional  Sureau  of  Standards, 
Boulder,  Colorado  80503 


cavity.  In  this  configuration,  transmission 
through  the  stripline  peaks  at  frequencies  for 
which  an  integral  number  of  half  wavelengths 
equals  che  stripline  length.  Measuring  che 
resonant  frequency  and  stripline  length  thus 
yields  the  phase  velocity  of  the  stripline. 

In  order  to  accurately  measure  che  cavity 
Q,  it  is  essential  chat  the  external  circuit 
does  not  significantly  load  the  scrlpllne.  To 
insure  this  condition,  capacitors  are  Inserted 
in  series  between  the  stripline  and  the  50  0 
coaxial  lines.  An  equivalent  circuit  of  che 
experimental  arrangement  then  cakes  che  form 
shown  in  Fig.  2.  The  capacitors  are  chip 
capacitors  of  a  few  pF  placed  directly 
against  che  contact  pads  on  the  sample.  By 
changing  the  size  of  these  capacitors  one 
can  determine  whether  or  noc  external  loading 
is  affecting  the  cavity  Q. 


If  che  stripline  is  characterized  by  a 
propagation  constant  y  and  a  characteristic 
Impedance  Z  ,  then  the  slnusoidal-sceady- 
scate  solution  of  che  equivalent  circuit 
shown  in  Fig.  2  is 

Vvs  ■  Rl  ~U1wCl  [:co5h  ^ 


! _ fo _ 

\  R^  +  l/lu,CL 


*■ 


*  1/lu.C, 


sinn 
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Fig.  1.  Geometry  of  the  stripline  samples.  The 
pattern  Is  1.4  cm  on  a  side. 
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Pig.  2.  Equivalent  circuit  of  the  experi¬ 
mental  arrangement.  The  reale tore  S.  replace 
the  SO  a  coax  end  SO  Q  load  resistors  of  the 
eveep  generator  and  oacllloacope.  The  balanced 
pals  of  strlpllnaa  nap  be  thought  of  ee  simple 
transmission  line  with  the  eaaa  propagation 
constant  ea  a  single  atrip  and  a  characteristic 
impedance  twice  that  of  a  single  atrip. 


where  l  is  the  length  of  the  atrlpllne.  This 
equation  defines  the  aaplltuda,  position,  and 
shape  of  the  resonances  observed  as  the 
oscillator  T,  la  swept  In  frequency.  For 
sufficiently  snail  losses  the  persneters  y 
and  Zg  can  be  approximated  by 


Y  “  a  ♦  lu  ^LC 

t 

(2) 

zo-  /l7c  , 

O) 

from  the  first  power  of  frequency  to  some¬ 
thing  approaching  the  second  power. 

While  interpretation  of  this  data  is  far  from 
certain,  sons  tentative  conclusions  can  be  reached 
The  possibility  that  radiation  losses  determine 
the  Q  can  be  elialnaced  on  the  basis  of  estimates 
of  the  radiation  from  various  stripline  dis¬ 
continuities9'  *  which  yield  Q's  about  five  orders 
of  nagnltuda  above  those  observed.  Of  the  ra- 
aalalng  two  types  of  loss,  dielectric  and  con¬ 
ductor,  we  believe  the  first  applies  at  4  K  and 
that  the  second  begins  to  appear  as  the  T  of  the 
Pb-Au  alloy  is  approached.  c 

attenuation  due  to  dielectric  losses  can  ba 
expressed  In  earns  of  the  loss  tangent  e"/t'. 


If  the  loss  tangent  la  independent  of  frequency, 
ea  is  coamon  for  many  materials,  than  a  is 
proportional  to  u  as  observed  here  at  4  K.  The 
attenuation  due  to  conductor  losses  can  be  ex¬ 
pressed  in  terns  of  the  surface  resistance 

“  “  2wZ  (RS,Pb  +  RS,»b)  '  (6) 

o 


where  a  is  the  attenuation  of  the  line  and  L  and  C 
are  the  aarlaa  inductance  and  shunt  capacitance  par 
unit  length.  The  quantltiea  L  and  C  can  ba  obtained 
experimentally  from  a  measurement  of  ths  phase 
velocity 

Vj  •  l/SZc  ,  (4) 

combined  with  a  direct  capacitance  measurement  aada 
using  a  low-frequency  bridge.  A  capacitance  aaasura- 
ment  aada  at  low  frequencies  is  sufficient  in  the 
present  case  as  the  dielectric  constant  of  Nb,0, 
proves  to  ba  tha  same  at  1011  Ha  as  at  107  1  3 

H*7.  Tha  one  parameter  of  Eq.  1  which  remains  to 
be  determined  la  tha  attenuation.  Provided  the 
external  loading  is  email,  a  fixes  the  widths  of 
tha  resonances  pradlcted  by  Eq.  1  and  a  least- 
squares  fit  to  ths  observed  line  shape  provides 
ao  accurate  measurement  of  a.  The  unloaded  Q 
is  related  to  a  by  Q  -  n*/2a i  where  n  la  the 
resonance  number. 


IV.  RESULTS 


Attenuation  measurenanta  were  aada  for  four 
strlpllnaa  for  which  tha  Pb-Au  alloys  of  tha  atrip 
vara  0,  4,  8,  and  162  Au  by  weight.  The  attenuation 
of  each  sample  was  measured  over  frequencies  between 
SO  and  500  Mflt  and  at  csmparaturaa  of  4  K  and  9S2 
of  T  for  tba  Pb-Au  alloy.  Data  for  tha  162  Au 
sample,  shown  in  Pig.  3,  typify  tha  results.  At 
4  K  tha  attenuation  of  all  samples  was  nearly 
identical  and  appeared  to  ba  proportional  to 
frequency.  As  tha  temperature  approached  T  „ 
(which  ranged  from  7.3  K  for  pure  Pb  to  7 . 0CK?B 
for  162  Au) ,  the  attenuation  began  to  increase 
rapidly  and  tha  frequency  dependence  changed 


FREQUENCY.  WHt 

Fig.  3.  Attenuation  of  tha  162  Au  sample. 
Circles  and  squares  indicate  experimental 
points  at  4  K  and  952  of  T  for  the  Pb-Au 
alloy  raapaetivaly.  Solid  lines  are  fits 
to  tha  data  based  on  tha  indicated  frequency 
dependencies. 
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where  w  la  tha  width  of  eba  scrip line.  For 
superconductors,  R_  la  vary  nearly  pro¬ 
portional  to  w2  ana  lacraaaaa  rapidly  aa  tha 
caaperaturs  approachaa  T  .  Thus.  It  la 
reasonable  to  auppoaa  that  tha  attaauatloo 
obaarvad  aaar  T  _  raaulta  from  a  com¬ 
bination  ol  losSAs0 In  tha  dlalactr 1c  and 
tha  Pb-Au  atrip.  Loaaaa  In  tha  Mb  grouad- 
plana  remain  snail  sine a  taaparaturaa  naar 

T  _  «  9.2  K  ara  navar  reached, 
c.ab 

An  experiment  that  confirms  the  assumption 
that  dlalactrlc  loaaaa  doalaata  at  4  X  lavolvaa 
two  samples  which  wara  praparad  dlffarantly 
only  la  that  tha  final  anodization  voltages 
wara  SO  V  and  100  V.  If  tha  two  aaaplaa  ara 
asamad  ldantical  axe  apt  for  dlalactrlc  chick- 
naaa.  chan  tha  only  paraaacara  co  changa  In 
aquae Iona  5  and  6  ara  V,  and  Z  .  In  aaaaura- 
aanta  of  tha  two  aaaplaa  V,  changed  by  a 
factor  of  1/.86  and  Z  changed  by  a  factor 
of  1/.S7  while  the  attenuation  changed  by  a 
factor  of  0.82.  Thua,  the  attenuation  at 
4  K  acaled  nearly  aa  1/V,,  a a  expected  for 
dlalactrlc  loaaaa. 


In  analyzing  tha  attenuation  data  to 
determine  tha  loaa  tangent  of  Nb,0,  and  che 
aurface  reeiatance  of  Pb,  wa  fit  the  data  at 
4  K  to  tha  font  a  -  c.u  jad  at  9SZ  of  I 
to  tha  for*  a  -  c.u+  t.ui  .  The  resultant’™ 
fitting  curvaa  for  tha  i6Z  Au  aaapla  ara  ahown 
In  Fig.  3.  Tha  values  of  c"/c'  ao  obtained 
ara  given  In  tha  following  cable. 
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Fig.  4.  Surface  raaiacance  of  the  Pb-Au 
alloy  aa  a  function  of  bulk  naan  fraa  path. 
Theoretical  raaulta  froa  Mactls-Bardeen 
theory  are  Indicated  by  the  aolld  line 
(local  Halt)  and  daahed  line  (excreaa 
anrmalouc  Unit,  dlffuaa  acactarlng.) 


Table  I.  Loaa  tangent  for  KbjOj 


Saapls 

T  ■  4  K 

T-  0.95  T 

c.Fb 

0Z  Au 

2. a  x  10"3 

3.1  x  10*3 

4Z  Au 

2.14  x  10-3 

2.7  x  10’3 

8Z  Au 

2.12  x  10-3 

2.6  x  iO-3 

16Z  Au 

2.08  x  10-3 

3.3  x  10“3 

Since  the  loaa  tangent  ahould  ba  Independent 
of  tha  Pb-Au  alloy,  che  ecatter  obtained  for 
the  higher  taaparatura  la  Indicative  of  an 
error;  aost  probably  due  to  the  lnsenaltlvlty 
of  the  fit  to  Cite  way  loaaaa  are  apportioned 
between  C.  and  C^,  but  perhapa  due  to 
experimental  error*  or  aa  error  In  Inter¬ 
pretation.  It  la  lntereatlag  to  note  that 
a  frequency- Independent  loaa  tangent  has 
alao  bean  obaarvad  by  Fuachillo  at  al. 
for  anodic  Mb-0.  at  room  taaparatura  and 
frequencies  above  S  z  1<J  Hz.  fuachillo  at 
al.  attribute  this  behavior  to  uncorralatad 
hopping  conduction  for  which  c"  la  pro¬ 
portional  to  caaperature.  If  chair  rooa 
taaparatura  value  la  extrapolated  to  4 Jt 
one  obtains  a  loss  tangent  of  1.4  x  10  , 

in  rough  agreeaent  with  our  aeasuraaents. 

The  surface  resistance  of  the  Pb-Au  alloy 
at  500  MHz  obtained  froa  our  fitting  procedure 
la  shown  In  Fig.  4  as  a  function  of  tha  bulk 
naan-free  path  l  for  normal  electrons.  To 

estlaate  t  we  assuae  that11 


oa/l  -  1015  Q_l  a”2  ,  (7) 

where  o  Is  the  bulk  conductivity  at  T  .  o  Is 
derived11  In  turn  froa  tha  chln-flla  conductivity 
by  estlaatlng  thickness  effects  under  tha  as-,, 
suaptlon  of  diffuse  scattering  ac  tha  surface1  . 

Also  shown  in  Fig.  4  arc  theoretical  results 
for  R.  at  500  MHz  and  a  reduced  ceaperature  of 
0.9S  obtained  by  applying  Mattls-Bardaen  theory 
to  Pb  In  the  local  Unit  (solid  line)  and  the 
excreaa  anomalous  Halt  with  diffuse  scattering 
(dashed  Line)1  .  Because  tha  Plppard  coherence 
length  E0  la  80  na  and  tha  observed  penetration 
depth  1  at  .95  T  ranged  froa  160  na  (for  OZ 
Au)  to  250  na  (f§fro  16Z  Au),  neither  tha  local 
Halt  (which  requires  l  <<  X  or  E  «  X)  or  the 
extraae  anomalous  Halt  (which  requires  1  »»  X 
and  Eq  >>  x)  can  be  legiclaately  applied  In 
the  present  case.  However,  tha  rough  estimates 
for  R,  so  obtained  confirm  that  the  observed 
losses  are  of  the  right  order  of  aagnitude  co 
be  due  co  superconducting  Pb. 

C0MCL0SI0N 

At  4  K  che  attenuation  observed  In  the  saaple 
striplines  can  be  explained  as  losses  In  the  Nb,0, 
dielectric  and  yield  a  loss  tangent  of  about 
2  x  10  .  Extrapolation  of  the  dielectric  losses  to 

'10  CHz  yields  an  attenuation  of  10  dB/a,  leading  to 
losses  acceptable  for  aost  nlcroclrcult  applications 
vhers  strlplins  lengths  are  1  ca  or  less. 


so 
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The  properties  of  nomul-eute  and  »upertonducting  tUiplines  are  calculated  aa  a  function  of  miniaturisacion.  For 
normal  conductors  the  Reuter-Sondheuiier  theory  is  applied  in  order  to  account  for  the  effects  of  finite  film 
thickness  and  mean  free  path.  For  superconductors  the  Mattie  Bardeen  theory  is  used  in  order  to  include  effects 
due  to  the  energy  gap.  Calculations  for  three  example  conductors,  copper  at  29$  K  and  4.2  K  and  niobium  at  4.2 
K.  examine  the  attenuation,  dispersion,  and  characteristic  impedance  of  sthplines  as  a  function  of  frequency  and 
dielectric  thickness.  Simulations  of  pulse  transmission  ate  used  to  evaluate  the  utility  of  the  example  sthplines  for 
high-speed  digital  applications. 

Key  words:  Copper:  niobium:  sthpline:  superconductivity;  surface  impede  nee. 


I.  introduction 


Striplines  interconnect  the  active  devices  of  both  semicon¬ 
ductor  and  Josephson-junction  high-speed  integrated  cir¬ 
cuits.  The  desirability  of  high  device  densities,  particularly 
for  digital  circuits,  provides  motivation  for  the  miniaturiza¬ 
tion  of  all  circuit  components,  including  striplines.  At 
present,  device  densities  for  semiconductor  logic  circuits 
having  100  ps  propagation  delays  are  limited  to  about  25 
gates  per  chip,  simply  because  each  gate  consumes  40  mV 
of  power  [1]‘.  In  contrast.  Josephson-junction  circuits  of 
comparable  speed  require  only  about  one  thousandth  as 
much  power  [2]  and  correspondingly  higher  device  densities 
are  anticipated  [3].  While  power  disaipation  presently  limits 
the  scale  of  integration  for  high-speed  semiconductor  cir¬ 
cuits.  it  has  also  been  noted  that  the  normal-metal  intercon¬ 
nections  used  in  such  circuits  degrade  much  more  rapidly 
with  miniaturization  than  the  superconducting  interconnec¬ 
tions  of  Josephson-junction  circuits  [4].  Here  we  present 
detailed  calculations  comparing  the  properties  of  normal- 
state  and  superconducting  striplines. 

The  stripline  geometry  to  be  considered  is  shown  in  figure 
1.  For  striplines  typical  of  semiconductor  circuits,  the 
conductor  and  dielectric  thickness,  d  and  j,  are  1  and  200 
Jim  respectively  [5],  while  for  superconducting  circuits  these 
dimensions  are  both  around  0.2  jam  (6).  Thus,  although  the 
conductors  are  of  roughly  the  same  thickness  in  both  tech* 
noiogies.  the  dielectric  thicknesses  differ  by  three  orders  of 

t  Partially  wpported  bv  the  Office  ot  Naval  Raaraich  under  contract  number 
\00fM4-7T-F-0048.  dated  *a*  1.  19T7. 
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1  Figures  in  brackets  indicate  literature  references  at  the  end  of  (his  paper. 
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Figure  1.  Stnpline  geometry  in  cmu-eecUon. 


magnitude.  The  thinner  dielectric  of  the  superconducting 
circuit  permits  higher  device  densities  for  two  reasons.  First, 
because  the  characteristic  impedance  of  a  stripline  is.  to  a 
good  approximation,  a  function  only  of  the  dielectric  constant 
and  the  ratio  of  strip  width  to  dielectric  thickness  n/s. 
shrinking  the  dielectric  layer  permits  smaller  strip  widths 
while  maintaining  the  same  impedance  level.  Second,  be¬ 
cause  adjacent  striplines  become  coupled  if  the  separation 
between  them  is  less  than  a  few  dielectric  thicknesses  [?]. 
thinning  the  dielectric  layer  permits  closer  spacing  of  the 
striplines  while  avoiding  crosstalk.  The  actual  strip  widths 
used  in  semiconducting  and  superconducting  circuits  are 
typically  100  Jim  [5]  and  5  jxm  [8]  respectively,  and  differ 
by  less  than  a  factor  of  10*  primarily  because  different 
impedance  levels  are  appropriate  for  the  active  devices  of 
the  two  technologies.  For  present  superconducting  circuits, 
crosstalk  can  probably  be  avoided  if  striplines  are  separated 
by  I  jim.  while  for  semiconducting  circuits  a  separation  of  at 
least  400  jim  is  required.  One  concludes  that  the  mimatun- 
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zatioo  of  amplifies  beyond  their  present  tire  in  semiconduc- 
tor  circuits  can  proceed  without  reduction  of  the  conductor 
thickness  but  not  without  thinning  the  dielectric. 

As  will  be  shown,  both  radiation  and  dielectric  loases 
become  relatively  lesa  important  compared  to  conductor 
losses  as  the  dielectric  thickness  is  reduced.  For  the  purpose 
of  this  study,  it  thus  proves  sufficient  to  focus  on  the  effects 
associated  with  imperfect  conductors,  assuming  a  lossless 
dielectric  and  neglecting  radiation.  The  effect  of  miniaturi¬ 
zation  is  considered  for  striplines  made  from  three  different 
conductors,  copper  at  295  K  and  4.2  K  and  niobium  at  4.2 
K.  These  conductors  represent  respectively,  a  normal  metal 
exhibiting  the  normal  and  anomalous  skin  effects  (at  the 
frequencies  of  interest)  and  a  superconductor. 

To  further  simplify  the  calculations,  we  consider  only  the 
TEM  mode  and  the  stripline  width  is  taken  to  be  much 
greater  than  the  dielectric  thickness,  w  »  s.  such  that 
fringing  fields  can  be  neglected.  These  assumptions  elimi¬ 
nate  dispersion  associated  with  the  discontinuity  in  dielectric 
constant  at  the  air/dielectric  interface,  an  important  source 
of  dispersion  for  w  comparable  to  or  smaller  than  s  [9,  10]. 
One  notes,  however,  that  for  the  TEM  mode  such  dispersion 
can  be  reduced  to  an  arbitrary  magnitude  by  overlaying  the 
stripline  with  a  dielectric  layer  of  sufficient  thickness. 

Under  the  above  assumptions  the  properties  of  a  stripline 
follow  almost  immediately  once  the  surface  impedance  of  the 
conductor  is  calculated.  In  section  2,  expressions  for  the 
surface  impedance  of  both  normal-metals  and  superconduc¬ 
tors  of  finite  thickness  are  reviewed.  These  expressions  are 
used  in  section  3  to  evaluate  the  attenuation,  phase  velocity, 
and  characteristic  impedance  as  a  function  of  dielectric 
thickness  assuming  a  fixed  conductor  thickness  of  1  fzm. 
Section  4  looks  at  the  propagation  of  short  pulses  over  various 
lengths  of  stripline. 

2.  Surface  Impedance 

The  surface  impedance  required  here  is  a  quantity  which 
abstracts  information  about  the  penetration  of  fields  into  an 
infinite  conducting  slab  of  thickness  d.  Taking  the  surfaces 
of  the  conductor  as  the  planes  z  m  0  and  z  m  d.  the  surface 
impedance  is  defined  for  s  sinusoidal  electric  field  £*4*. 
mV1***  and  current  density  J  ,(r.  by 


with  the  boundary  condition  that  the  magnetic  field  be  zero 
at  z  m  d  as  ia  appropriate  for  a  stripline  with  it>  »  s.  The 
real  part  of  Z,,  called  the  surface  resistance,  accounts  for 
stripline  losses  and  the  imaginary  part,  the  surface  reactance, 
contributes  to  the  stripline  inductance. 


Evaluation  of  the  surface  impedance  for  copper  at  295  K 
and  niobium  at  4.2  K  ia  simplified  by  the  fact  that  J  and  E 
are  related  by  the  local  equation 

J  -  o-E.  (2.2) 

where  <r  ia  the  complex  conductivity.  This  local  equation  can 
be  assumed  because  the  mean  free  path/  for  normal  electrons 
ia  short  compared  to  all  other  dimensions  in  the  problem. 
Combining  eq.  (2.2)  with  Maxwell's  equations  yields  for  the 
surface  impedance  [11] 

Z,  *  (uofto/o^coth  (2.3) 

The  surface  impedance  thus  follows  once  the  conductivity  is 
known.  For  copper  <r  is  simply  a  real  constant.  For  supercon¬ 
ducting  niobium.  <r  is  a  complex,  frequency-dependent 
quantity  which  we  take  to  be  of  the  form  given  by  Mattis  and 
Bardeen  [12].  The  Mattis-Bardeen  equation  yields  the  con¬ 
ductivity  of  a  superconductor  at  a  given  frequency  and 
temperature  provided  two  material  parameters  are  known; 
the  energy  gap  parameter  A  and  the  normal  state  conductivity 
at  the  superconducting  transition  temperature  <r„.  A  more 
complete  account  of  the  application  of  Mattis-Bardeen  theory 
to  superconducting  striplines  has  been  given  elsewhere  [13]. 

In  the  derivation  of  eq.  (2.3)  for  the  case  of  infinite 
conductor  thickness,  it  is  observed  that  the  field  amplitudes 
decay  exponentially  with  distance  into  the  conductor.  For  a 
normal  metal  the  characteristic  decay  length  is  the  classical 
skin  depth 

K  -  (oMi.tr/2r »«.  (2.4) 

Because  5r  decreases  with  increasing  frequency  there  exists 
a  frequency  above  which  5 r  <  /  and  a  local  relation  between 
J  and  E  can  no  longer  be  assumed.  At  such  high  frequencies 
the  skin  effect  is  said  to  be  anomalous.  For  a  superconductor, 
the  decay  length  at  frequencies  less  than  the  energy  gap 
frequency,  2SlA.  is  the  penetration  depth  [14] 

\  *  (Xcoth  (A/2kr)/wju(>A<rll)1'*.  (2.5) 

which,  in  contrast  to  6r.  is  frequency  independent.  In  the 
process  of  miniaturizing  striplines,  the  depth  of  field  penetra¬ 
tion  plays  an  increasingly  important  role  as  the  dielectric 
thickness  is  reduced. 

Because  the  mean  free  path  of  copper  at  4.2  K  can  be  long 
compared  to  both  Sf,  and  d.  the  relation  between  J  &nd  E 
assumes  the  non-local  form  [15] 

J(r)  -  ^  j  dipp(p  ■  E(r  +  p))  p~*r~*u.  (2.6) 

where  we  have  assumed  that  the  relaxation  time  r  is  small 
compared  to  the  inverse  frequency,  tor  «  1.  Evaluation  of 
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tilt  aurface  impedance  ia  thia  cat*  requires  an  additionel 
boundary  condition  regarding  the  scattering  of  electrons  from 
the  conductor  aurface.  Of  the  taro  simple  lhnita.  diflute 
scattering  tad  tpteular  reflection,  the  farmer  aetata  to 
provide  the  better  agreement  with  experiment  [15]  aad  haa 
beta  aeamaed  ia  the  preitnt  calculations.  Equation  (2.4) 
aad  Maxwell's  equations  applied  to  the  surface  problem  with 
diffuse  scattering  yield  [15]. 


€L 

dr* 


ia/-*  f‘ di'EAs')KUx'  -*)//) 


<2.7) 


where  a  “  */t  /*/  8r*.  Reuter  aad  Sondheimer  [16]  aoived 
thia  equation  in  the  limit  of  infinite  d  and  obtain  for  the 
aurface  impedance 


«»1.  (2.9) 

The  expressions  for  the  sutface  impedance  of  a  aonaai 
metal  given  above  are  valid  either  far/  <<  8*.  d  (eq  (2^3)1 
or  far  arbitrary  /  and  d  »  ■  (eq  (2.8)).  The  numerical 
solution  of  eq  (2.7)  aad  calculation  of  the  aurface  impedance 
far  arbitrary  /,  8^,  and  d  ia  diacueaed  ia  the  appendix  The 
raanlting  general  program  waa  uaed  to  obtain  ail  raoulia  far 
normal  matala  pwatntad  here. 

Although  the  Soldo  do  nor  foil  off  ia  o  praciaoiy  exponential 
way  far  the  eiuunalrtue  akin  effect,  it  remaiaa  poaeihie  to 
deSae  a  characteristic  penetration  length.  In  general,  the 
diataace  over  which  the  field  amplitude  decay*  is  given  by 


<*pL,  R*{Z,\ 


(2.10) 


2. 

kit) 


Ml  +•  ittk(t)/f*) 


r 


(2.8) 


jj[(1  +  t*)  taa-Ht)  —  t]. 


For  a  «  1.  (normal  akin  effect),  thia  axpwaaion  reduces  to 
the  d  ■  »  limit  of  eq  (2.3).  For  a  »  1  (extreme  annmeloua 
altia  effect).  Reuter  aad  Sondheimer  obtain 


Ia  the  local  limit  thia  reduces  to  8r  for  a  normal  metal  and  to 
h  for  a  superconductor.  In  the  extreme  tnnmakwis  limit,  eq 
(2.9)  may  be  uaed  to  obtain 

ouseua 

8.-~3-(Sc*/)"».  (2.11) 

ao  that  the  akin  depth  Sot  l  >  Sr  is  greater  than  Sr. 
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8.6  x  ltr* 

The  material  parameter*  of  copper  and  niobium  uaed  here 
are  preeeoted  in  table  l.  The  conductivity  of  capper  at  room 
temperature  [17]  is  limited  by  phonon  scattering  aad  ail 
specimens  have  nearly  the  aame  value.  At  4.2  K  scattering 
bom  impurities  and  defects  dominate  aad.  although  the 
conductivity  can  be  more  than  10*  timet  that  at  295  K  for 
specially  prepared  sample*,  a  residual  resistance  ratio  of  10* 
ia  common  far  off-the-shelf  wires  [18].  Here  we  aaeume. 
perhaps  optimistically,  that  the  conductivity  at  4.2  K  is  500 
timet  the  room  temperature  value.  The  mean  bee  path  is 
derived  from  the  fact  that  ail  ia  approximately  constant  for  a 
given  material  and  haa  a  value  of  1.54  x  10**  ft*1  a~*  for 
copper  [11].  The  relaxation  time  ia  calculated  bom 

r -!/»,.  (2.12) 

where  the  Fermi  velocity  is  1.6  x  10*  m/s  for  copper  [ll]. 
Material  parameter*  for  superconducting  niobium  are  baaed 
’  on  the  thin-film  measurements  of  Hcttkela  and  Kircher  [19J. 


The  measured  value  of  am.  1.59  x  10*  fl-1  m”\  wo# 
modified  slightly  ao  that  eq  (2.5)  ia  exactly  satisfied. 

The  real  and  imaginary  pans  of  the  aurface  impedance  of 
the  three  example  metals  are  shown  in  figure  2  for  s  metal 
thickness  of  1  fu a.  Because  the  mesa  bee  path  is  much 
shatter  than  the  film  thickness,  the  surface  impedance  of 
copper  at  295  K  ia  described  by  the  local  formula,  eq.  (2.3). 
for  frequencies  below  about  3  x  10“  Ha.  Above  this 
frequency,  the  classical  akin  depth  is  smaller  than  the  mean 
free  path  and  the  akin  effect  ia  anomalous.  At  sufficiently 
low  frequencies  8,  ia  large  compared  to  the  conductor 
thickness  and  eq  (2.3)  reduces  to 

Z,»—+  «iaq*.rf.  l«  d«  &t.  (2.13) 

<Tu  3 

This  equation,  which  results  when  the  electric  field  is 
uniform  across  the  conductor  thickness,  accounts  far  the 
impedance  of  copper  at  295  K  and  below  about  4  x  10*  Ha 
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Ficlrk  2.  The  real  (solid  line f  and  imaginary 
(dashed  line t  pans  of  the  surface  impedance  for  l  pm 
thick  films  of  copper  at  295  K  and  4.2k  and  ntohium 
at  4.2  K. 

Tht  wnU»  for  coppor  at  4.2  K  end  Impaemtm  a heme 
10*  Hi  ««v  be  iMccumr  dm  to  the  ommirnm  o l  telmsSinen 
ume  effects. 


shown  in  figure  2.  At  higher  frequencies  Sr  becomes  small 
compared  to  d  and  eq  2.3  reduces  to 

Z,  -  1 1  +  i)  ^  I  «  6.  «  d.  12. 14) 

In  this  limit  the  electric  Held  decays  exponentially  to  a  small 
value  befote  teaching  the  back  surface  of  the  conductor. 

At  4.2  K  the  mean  free  path  of  copper  is  sufficiently  long 
that  5r  is  comparable  to  l  even  at  2  X  10*  Hz.  Thus,  the 
anomalous  skin  effect  is  exhibited  over  the  entire  frequency 
range  of  figure  2.  The  relaxation  time  is  sufficiently  long  at 

4.2  K  that  (■ ir  is  unity  at  1010  Hz  (compared  to  7  x  10**  Hz 
at  room  temperature!.  Pippard  [15]  has  shown,  however,  that 
in  the  extreme  anomalous  limit,  relaxation-time  effects  are 
important  only  when  u n  is  considerably  greater  than  unity 
and  in  the  present  case  need  not  be  considered  for  frequen¬ 
cies  less  than  about  10u  Hz.  Since  the  relaxation  time  has 
been  neglected,  the  results  shown  in  figure  2  for  copper  at 

4.2  K  are  probably  not  accurate  above  this  frequency.  As  at 
room  temperature,  the  surface  impedance  at  4.2  K  depends 
critically  on  whether  the  skin  depth  is  greater  or  less  than 
the  conductor  thickness.  The  break  point  occurs  when  5,  “ 
d  or  about  2  x  10*  Hz.  Below  this  frequency.  Z ,  is  given  by 
an  equation  similar  to  eq  (2. 13)  in  which  the  bulk  conductiv¬ 
ity  <7  is  replaced  by  an  apparent  conductivity  a  that  accounts 
for  the  effect  of  finite  conductor  thickness  on  the  mean  free 
path  [20].  In  the  present  case,  the  apparent  conductivity  is 
0.13  times  the  bulk  conductivity.  At  frequencies  above  2  x 
10-  Hz  we  have  Sr  «  L  d  and  eq  <2.9i  applies.  As  can  be 
easily  verified  from  eq  (2.91,  the  surface  impedance  in  this 
region  depends  only  on  ail.  a  constant  of  the  material,  and 


thus  does  not  depend  on  the  somewhat  arbitrarily  chosen 
residual  resistivity  ratio.  If.  for  example,  the  resistivity  ratio 
had  been  taken  as  5000  rather  than  500.  the  only  significant 
change  in  the  surface  impedance  would  be  a  six-iold  decrease 
in  £«{*,]  below  2  x  10*  Hz.  These  considerations  have  been 
discussed  in  detail  by  Keyes  et  ai.  [4]. 

Below  the  energy  gap  frequency.  7.2  x  10"  Hz  for  Nb  at 

4.2  K.  the  surface  resistance  of  a  superconductor  is  orders  of 
magnitude  smaller  than  that  of  a  normal  metal.  In  a  etnpline 
the  smaller  surface  resistance  results  in  lower  attenuation. 
Also,  the  surface  reactance  below  the  energy  gap  vanes  as 
oi.  making  it  appear  exactly  as  an  inductance.  This  property 
yields  a  stripline  with  very  low  dispetsion.  .Above  the  energy 
gap,  a  superconductor  behaves  like  a  normal  metal  of 
conductivity  am 

3.  Propagation  Constant  and  Characteristic 
Impadanca 

The  propagation  constant  y  and  characteristic  impedance 
Z  „  of  a  stripline  can  be  expressed  in  terms  of  the  series 
impedance  Z  and  shunt  admittance  T  of  a  unit  length  of  line. 


y-V/V.  (3.1) 

Z0  «  y/I/Y.  (3.2) 

where  Z  and  Y  for  iv  »  s  are  in  turn  [11] 

s  2 

Z  =*  la*! „  -  -t-  —  Z,.  (3.3) 

K  If 

Y  »  io*«€,,—  .  (3.4) 

s 


The  first  term  of  Z  is  the  inductive  impedance  associated 
with  the  magnetic  field  in  the  dielectric  region  and  the 
second  term  accounts  for  penetration  of  fields  into  the 
conductor.  Y  is  the  capacitive  admittance  between  the  strip 
and  the  ground  plane.  In  the  following  calculations  a  relative 
dielectric  constant  e  of  4  is  arbitrarily  assumed.  Note  that  y 
is  independent  of  w  while  Z„  is  proportional  to  1/ic.  Thus. 

both  y  and  ~  Z  „  can  be  calculated  without  specifying  the 

stripline  width.  The  power  attenuation  in  decibels  per  length 
a^B  and  the  phase  velocity  v4  are  related  to  the  real  and 
imaginary  parts  of  y  by 

a4B  *  Cita  fi<*{y].  (3.5) 

its  *  ui//mfy],  (3.6) 

where  £*8  *  20  log, or.  The  properties  of  a  stripline  thus 
follow  from  the  surface  impedance  of  the  conductor. 
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The  phmt  velocity,  attenuation.  and  charactariatic  imped* 
ance  of  a  room  temperature  copper  stripiine  are  shown  in 
figure  3  for  a  conductor  thickness  of  1  fim  and  a  dielectric 
thickneee  ranging  from  200  jam  to  0.2  iim.  Aa  the  dielectric 
thickneaa  ia  reduced  toward  0.2  >am.  the  line  becomea  highly 
diaperaive  and  looey  and  develop*  a  charectenetic  impedance 
far  horn  the  real,  frequency* independent  impedance  of  an 
ideal  line.  Aa  a  function  of  frequency  the  degradation  of  a 
fixed  length  of  line  due  to  loee  and  diaperaion  increaaee  with 
increaaing  frequency.  Thia  increaae  in  diaperaion  ia  leaa  than 
obvious  since  the  phaae  velocity  approache*  a  conatant  at 
high  frequencies  but  will  be  explained  preaentlv.  Problem 
with  impedance  matching  are  reduced  at  high  frequenciea  aa 
the  characteriatic  impedance  approache*  a  real  conatant. 

A  better  understand  inf  of  the  result*  shown  in  figure  3  can 
be  obtained  by  reviewing  approximate  expreaaion*  for  y  and 
Z  „  applicable  in  the  four  region*  bounded  by  dotted  line*. 
These  regions  are  defined  by 

I.  d.  Js3«  Sr 

II.  d  «  Sr« 

III.  d 

IV.  j  «  Sr«d 

and  reault  from  the  two  posaible  form*  for  the  surface 
impedance,  eqs  (2.13)  and  (2.141.  combined  with  two 
possibilities  for  the  series  impedance  of  eq  (3.3),  either 
resistive  or  reactive.  Of  the  four  regions.  I,  II.  and  III  are  of 
principal  interest  here  since  condition  IV  is  satisfied  only  for 
the  0.2  itm  dielectric  case  over  a  narrow  range  of  frequen¬ 
cies. 

The  asymptotic  forms  for  the  phase  velocity  in  the  regions 


of  interest 

are 

I-  r* 

lojuia 

.  l 

d* 

1  sd 

(3.7) 
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1  ~  ~ 
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(3.81 

II.  <  « 

1  ~  3 

S 

8  s!dlJ 

III.  r* 

l 

1 

8rl 

vW«o 

2 

Tj 

(3.9) 

where  correction  terms  have  been  indicated  in  brackets.  As 
figure  3  indicates,  approximations  I.  H.  and  III  become  valid 
at  successively  higher  frequencies.  In  the  lowest  frequency 
range  r*  is  proportional  to  -Jw  and  goes  to  zero  at  zero 
frequency.  If.  however,  one  considers  the  phase  shift  that 
results  because  the  low  frequencies  travel  at  a  velocity  less 
than  the  high-frequency  asymptotic  velocity  then  for  a  fixed 
length  of  line  t  this  phase  shift 


FREQUENCY,  Nz  ' 

Figure  3.  Phase  velocity,  attenuation,  anti  real 
(solid  line >  and  imaginary  i dashed  linei  parts  of  the 
characteristic  impedance  of  a  295  K  copper  stnpline 
Utah  l  pm  conductor  thickness. 

Results  ace  shown  lot  a  relative  dielectric  conatant  of  4 
and  dielectric  thicknesses  of  200  pm  <ai.  20  pm.  (b)  2  pm. 
(ck  and  0.2  pm  fd).  Dotted  tinea  divide  the  inpha  of  phaae 
velocity  and  attenuation  into  regions  where  simple  approxi¬ 
mate  equations  hold.  Circles,  tquates.  and  triangles  matt 
points  at  which  ■  i.  ■  yn.  and  *  d.  respectively . 
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t3. 15i 


A* -«/(  — - —)■ 

also  goes  to  aero  at  aero  frequency  as  V7I  Thus,  dispersion 
ai  low  frequencies  is  not  the  problem  that  figure  3  suggests. 
At  somewhat  higher  frequencies  6,  becomes  leas  than  yST 
and  approximation  11  applies.  In  this  region  c«  approaches 
the  asymptotic  value  (1  —  d/3»i/y>,*«,  with  a  correction 
term  which  decreases  as  with  inc teasing  frequency. 
Finally,  as  5,  becomes  less  than  d  as  well  as  j.  r«  approaches 
l/VMeM*  with  a  correction  terra  of  -8r/2j.  Because  this 
term  goes  as  l/v  a  the  phase  shift  Ab  goes  as  v  ta  and 
dispersion  effects  inc  tease  with  m  just  as  at  very  low 
frequencies.  Also  note  that  the  correction  terra  is  proportional 
to  1/a  so  that  thin  dielectrics  enhance  dispersion. 

The  attenuation  can  be  approximated  in  the  regions  of 
interest  by 


I.  (*4M  ■ 


II.  a«  - 


III.  a„  - 


3.101 


(3.11) 


(3.12) 


As  a  function  of  frequency.  a4S  increases  as  y«5  at  very  low 
and  very  high  frequencies  with  a  frequency  independent 
plateau  in  between.  As  a  function  of  dielectric  thickness. 

always  increases  with  decreasing  dielectric  thicknesa. 
either  as  1  /y'soraa  1/a. 

Lastly,  the  characteristic  impedance  is  approximated  by 


(3. 13) 


<3. 14i 


In  the  limit  of  low  frequencies  the  real  and  imaginary  pairs  of 
Z .both  diverge  as  l/v m.  The  significance  of  this  divergence 
relates  to  the  reduction  of  the  transmission  line  to  a  simple 
resistor  at  low  frequencies.  For  br<t.d  the  real  part  of  Z  . 

approaches  its  high  frequency  aeymptote.  -  v  hJ**,.  and 

M 

the  imaginary  pair  goes  to  sero.  The  correction  term  to  the 
real  part  is  3r/2i  to  that  the  dielectric  thickness  directly 
influences  the  ability  to  match  impedances  over  a  wide  range 
of  frequencies. 

The  properties  of  copper  and  niobium  aurpiutet  at  4.2  k 
are  shown  in  figure  4.  again  for  a  conductor  thickness  of  1 
pm  and  dielectric  thicknesses  ranging  from  200  pm  to  0.2 
ixm.  Comparing  copper  at  295  k  and  4.2  k  shows  that  the 
300  fold  increase  in  conductivity  which  might  result  from 
cooling  leads  to  substantial  improvements,  especially  si  low 
frequencies.  These  improvements  becomes  less  signif  icant  at 
frequencies  approaching  3  x  I015  Hi.  where  the  skin  effect 
becomes  anomalous  even  for  copper  at  295  k.  The  superior¬ 
ity  of  a  superconducting  stnpline  at  frequencies  below  2XA 
is  at  once  apparent  from  figure  4.  In  this  region  i  a  is 
independent  of  frequency  .  aa  is  orders  of  magnitude  below 
the  attenuation  for  a  normal  conductor,  and  the  characteristic 
impedance  is  a  real  constant.  Above  7.2  *  10"  Hi  all  of 
these  desirable  properties  disappear  as  the  superconductor 
begins  to  behave  like  a  normal  conductor. 

As  for  the  normal  skin  effect  case  al  293  k.  simple 
approximate  formulas  can  be  derived  for  the  properties  of  a 
stnpline  in  the  extreme  anomalous  limit.  The  regions  <>l 
approximation  are  defined  as  for  the  normal  skin  effect 
except  that  6r  is  replaced  bv  6..  These  regions  csn  be 
located  in  figure  4  bv  using  the  points  marked  with  circles 
1 6,  "*  si.  squares  (5.  »  ylal.  and  tnangies  ifi,  *  di  as 
guides. 

At  frequencies  less  than  2A/5  simple  approximations  also 
exist  for  the  parameters  of  a  superconducting  stnpline  [11. 
13).  The  phase  velocity  is 

I  /  2A  -A  -■ 1 

iy  *  — ;  ■  I  1  f  —  coth  -  I  .  u>  <<  2A  fi.  (3.  lbi 
V >*•«*€..  1  *' 


While  this  equation  is  similar  to  eq  (3.9l  for  a  normal  metal, 
it  differs  critically  in  that  A.  as  opposed  to  5r.  is  independent 
of  frequency.  Thus,  as  s  decreases  ia  also  decreases  but 
remains  frequency  independent.  The  attenuation  is 
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*(coth-  +  ^U. 

\  X  stnh  *d/X/ 


2X  L  /A  -•'* 

■e  — coth  -  I 

i  V 

(3.17) 

u  «  2 A/A. 


where  (7,  it  the  reel  part  o f  the  conductivity.  Becaute  <7|  it 
only  weakly  frequency  dependent,  a  it  nearly  proportional 


'•tJlXIO'.  -I 


Fwi’U  4.  Graphs  aamiogaut  to  figure  3  far  copper 
at  4.2  K  > taUd  timet  except  loop  Hashes  tar  lm(Z9J) 
amd  iip4ww  at  4  .2  k  '  than  Hashes  except  Hats  jar  im 

(Z«l> 

f or  copper,  curl**  iumwi.  saS  vmi^m  ~nirt  th* 
pMM  a  which  A,  ■).  «7.  anrt  4,  •  A .  iNpni**h 


to  lit*.  For  the  ethplines  conaidered  here  X  <<  s.  d and  the 
attenuation  it  proportional  to  l/i.  The  characteriatic  imped* 
ance  it 

Z.  -  -  %/—  (  1  +  — coth  -)  '  to  «  2 A/*-  (3.18) 

which  like  t-«  changes  with  i  but  remains  frequency  inde¬ 
pendent.  The  imaginary  part  of  Z  „  is  negligible.  Thus,  for 
«s  <<  2 AJh.  the  only  serious  deterioration  of  superconduct¬ 
ing  stnpline  properties  with  miniaturization  involves  the 
attenuation. 

Before  discussing  how  the  transmission  of  pulses  is  af¬ 
fected  by  stripiine  properties,  we  briefly  examine  the  two 
sources  of  loss,  dielectric  and  radiation,  which  were  ne¬ 
glected  in  the  above  calculations.  Dielectric  losses  occur 
when  «  has  a  non-zero  imaginary  part  such  that  e  *  t  + 
it".  Assuming  e”  <<  f  and  a  perfect  conductor  (Z,  *  Ol. 
the  attenuation  due  to  dielectric  loss  is 


\ 
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__ _  £* 

c«  Vm»«' 


(3.19) 


For  many  materials  t  and  (  are  approximately  independent 
of  frequency  and  a4S  increases  roughly  as  to.  For  good 
dielectrics  «"/«’  is  typically  10-4  at  10,#  Hz  [21)  yielding,  for 
t  m  4.  an  attenuation  of  0.2  dB/m.  Because  this  attenuation 
is  independent  of  i.  the  conduction  losses  will  always 
dominate  for  small  enough  s.  .As  ran  be  seen  from  figures  3 
and  4.  conduction  losses  at  10“  Hz  are  greater  than  0.2  dB/ 
m  for  all  the  normal-state  lines  considered  and  for  the 
superconducting  lines  with  s'  “  2  Mm  and  0.2  gm.  Because 
0.2  dB/m  represents  a  small  loss  over  distances  typical  of  a 
microcircuit,  neglecting  dielectric  loss  will  not  affect  the 
conclusions  of  section  4. 

A  comparison  of  radiation  and  conduction  losses  is  com¬ 
plicated  by  the  fact  that  radiation  losses  cannot  be  converted 
to  a  loss  per  unit  length,  radiation  being  associated  only  with 
discontinuities  in  the  line.  In  order  to  estimate  the  impor¬ 
tance  of  radiation  we  consider  the  Q  of  a  half-wave  open-end 
stripline  resonator.  Radiation  from  this  structure  has  been 
calculated  bv  a  number  of  authors,  most  recently  by  van  der 
Pauw  [22].  Assuming  w  >>  s  >>  &.  van  der  Pauw's  result 
for  the  O  due  to  radiation  losses  reduces  to 


<>» 


i  it1?*  f  (2«  -  nf 

2  {  3 


(3.20) 


where  »// y  /*«««  is  the  resonant  frequency  and  (  is 

the  length  of  the  cavity.  The  y  due  to  conduction  losses  for 
the  same  resonator  is 
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If  the  characteristic  impedance  is  held  fixed  (i c/s  constant) 
then  Q  K  increases  with  miniaturization  and  Q  c  decreases. 
Davidheiser  [23]  recently  noted  this  fact  and  determined  the 
dielectric  thickness  at  which  the  Q  of  a  superconducting 
resonator  is  maximum.  Since  the  ratio  (Jk/Qc  decreases  with 
increasing  ut0  for  both  the  normal  and  superconducting  cases, 
there  is  a  resonant  frequency  above  which  radiation  losses 
dominate  conduction  losses.  This  breakpoint  frequency  is 
given  in  table  II  for  the  various  striplines  considered  here 
assuming  w/s  m  10.  For  the  0.2  and  2  /im  dielectrics,  the 
breakpoint  occurs  at  frequencies  higher  than  those  of  practi¬ 
cal  interest  and  radiation  losses  can  be  neglected.  For  the  20 
and  200  jim  dielectrics,  radiation  begins  to  become  impor¬ 
tant  at  frequencies  around  1010  Hz.  Since  the  wavelength  is 
1.5  cm  at  10l°  Hz  for  e  *  4.  microcircuit  striplines  could  be 
of  the  right  length  to  approximate  a  half-wavelength  antenna. 


While  the  load  at  the  ends  of  a  microcircuit  stripline  is  not 
infinite,  the  radiation  from  a  matched  or  shorted  end  is  a 
sizable  percentage  of  that  from  an  open  end  [24],  Thus 
radiation  losses  may  well  be  important  at  gigahertz  frequen¬ 
cies  for  dielectric  thickness  greater  than  20  g im. 

A  further  consideration  omitted  from  the  present  calcula¬ 
tions  regards  the  propagation  of  modes  of  higher  order  than 
the  assumed  TEM  mode.  Such  higher-order  modes  become 
possible  at  frequencies  sufficiently  high  that  a  half  wave¬ 
length  is  comparable  to  the  cross-sectional  dimensions  of  the 
stripline  (10).  Assuming  w  >>  s  the  cutoff  frequency  is 

(ur  =  tr/if  vV»««o-  (3.221 

such  that  for  €  =  4  and  w  *  200  gm  the  cutoff  is  at  4  x 
1011  Hz.  Thus,  even  for  striplines  of  large  dimensions  in  the 
context  of  microcircuits,  the  higher-older  modes  are  im|x>r- 
tant  only  at  frequencies  above  those  of  immediate  interest. 


Tails  II.  Frequency  i  Hz)  at  which  radiation  and  conduction  Losses  are  equal  for  a  hoif-waie  ooen-ended  stnpltne  resonator  with  nils  m  10  and  «  ■  4. 


Cu.  295  K 

2.0  x  10“ 

1.3  x  10“ 

8.0  x  10'“ 

5.0  x  10* 

Cu.  4.2  K 

1.7  x  10“ 

8.9  x  10“ 

4.6  x  10“ 

2.4  x  10> 

•  Nb.  4.2  K 

2.6  x  10“ 

1.7  x  10“ 

1.4  x  10» 

3.2  x  10* 

I 

1 

1 

IV.  Puls*  Propagation 

*1 

1  O 

The  evaluation  of  a  stripline  from  the  standpoint  of  pulse 
transmission  follows  from  a  knowledge  of  y  and  Z  „  but  is 
sufficiently  complicated  that  simulations  prove  valuable.  In 
this  section  we  consider  the  propagation  of  Gaussian  pulses 
over  a  stripline  of  length  t  using  the  circuit  shown  in  figure 
5.  The  source/load  impedance  RL  is  matched  to  the  asymp¬ 
totic  high-frequency  value  of  Z  „  in  the  case  of  a  normal-state 
stripline 


'll  \ 


STR IPlIflE 

FiCL'RE  5.  Circuit  for  analysis  of  pulse  transmt ssion . 
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where  the  standard  deviation  T  measures  the  pulse  width, 
and  in  the  superconducting  case  is  matched  to  the  low-  Taklng  th„  Fouri„  tran9forni  obtams 
frequency  impedance 

—  ,,  Vs(co)  -  t2ir)*'*  (4.4) 

Rl  ~  ~  yi  —  (  1  +  — coth  —  I  (4.2) 

’  **»  V  5  A  so  that  the  frequencv  spectrum  is  also  Gaussian  with  a 

standard  deviation  Aui  ”  1/r.  Solving  the  circuit  of  figure  5 
The  degree  to  which  a  voltage  pulse  Vs  is  faithfully  repro-  ,or  the  sinusoidal  „eady-»tate  yields 
duced  across  the  load  is  affected  both  by  y.  in  the  dispersion 

and  attenuation  of  the  line,  and  by  Z  in  reflections  at  the  V  t(o)V  l  xlu>)  m  I  2  cosh  yt 

interface  between  the  line  and  the  source  or  load  due  to  '  ,  ,  -. 

t  Z  R  \  \  - 1 

imperfect  impedance  matching.  +  j  ZS.  _  Ik  |sln|1  | 

The  Gaussian  pulses  considered  are  of  the  form  \  Ri  Z„ 
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Ficvas  6.  Pwpapmiaa  <y  c toiau putt* far*  copp ar  mnpiim t  m  19S  K 
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Palm  ar»  i!m  W  lav  laapka  at  atm  nH  MW.  0.1  ca  ilm  IiOh.  1  « 
i<Mk  ami  10  am  i**l  Cm*m  ataah  amaamall*  aaamda  aa*  *a  aaaa  laapk  tmaa 
ham  mat  haam  paphad.  TW  ta ma  atap a  a  a hatad  b a  m  mamma  l,  “  flap !•>  '  ail 
*4*1  •  1.}  a  IIP  ava  « tha  hapaiaapaaaan  p ama  uriwy. 


so  that  the  pulse  received  at  the  load  take*  the  form 


1 1<»)  -  <2/ tr)1'*  tV.  I  dot  e-1'2'* 
x  Rf  [e~(  2co*h  „„h  yf)  '] 


(4.6) 


Numencal  evaluation  of  the  above  integral  was  used  to 
determine  the  shape  of  a  pulse  after  (ravening  various 
lengths  of  stripline.  Since  R  L  and  Z  „  are  both  proportional  to 
1/u-  these  calculations  do  not  depend  on  the  stripline  width. 

The  results  of  such  simulations  for  copper  stnplines  at  295 
K  and  4.2  K  are  shown  in  figures  6  and  T.  For  a  given 
dielectric  thickness  and  initial  pulse  width,  pulses  are  shown 
as  they  would  appear  after  traversing  stripline*  of  length  0. 

.  ■  0. 1.  1.  and  10  cm.  The  time  origin  for  each  trace  hat  been 
shifted  by  an  amount 


r,  ■  f/t*(»). 


(4.7) 


Ficimt  Gnp/u  vaaiagaau  to  fa*n  6  far  a  cappar  luaa  at  4J  K. 


such  that  a  pulse  traveling  at  the  asymptotic  high-frequency 
phase  velocity  t*(*l  would  alwavs  be  centered  at  the  origin. 
The  results  shown  in  the  left  column  of  figures  6  and  7  are 
for  a  pulse  of  r  *  50  ps  or  roughly  at  the  upper  limit  of 
speed  for  present  digital  circuits.  If  we  arbitrarily  require 
that  a  pulse  be  received  with  at  least  80  percent  of  its  original 
amplitude  to  be  useful  in  digital  applications,  then  we  see 
from  figure  6  that  for  s  *  20  tun.  2  Mm.  and  0.2  tun  useful 
stripiine  lengths  are  less  than  10  cm.  1  cm.  and  0.1  cm. 
respectively.  Thu*,  the  performance  of  a  room-temperature 
copper  stripiine  with  a  dielectric  thickness  of  2  tun  or  less 
can  be  a  limiting  factor  even  for  lengths  typical  of  microcir¬ 
cuit  dimensions.  At  4.2  K  this  situation  is  improved  by  an 
order  of  magnitude  in  that  a  dielectric  thickness  of  onlv 
about  0.2  t*m  yields  lines  of  sufficient  quality  for  digital 
microcircuits. 

Theoretical  arguments  indicate  that  it  may  be  possible  to 
generate  pulses  with  widths  as  short  ss  1  p>  using  supercon¬ 
ducting  devices  [25].  The  propagation  of  such  s  pulse  on 
normal-state  striplines  is  shown  in  the  right  column  of  figures 
6  and  7.  using  a  r  of  0.5  ps  to  give  a  full  width  of  about  1  ps. 
In  going  from  r  ■  50  pi  to  r  *  0.5  ps  one  note*  that  the 
distortion  of  the  shorter  pulse  is  much  the  same  ss  that  of  the 
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longer  after  propagating  only  one  tenth  the  distance.  Thus, 
for  t  »  O.S  pt,  a  copper  line  at  295  K  suitable  for  digital 
Microcircuits  requires  more  than  20  pm  of  dielectric  and  at 
4.2  K  requires  more  than  2  pm. 

Further  insight  into  the  transmission  of  pulses  is  gained  by 
examining  four  limiting  situations  in  which  the  transfer 
function  of  eq  (4.5)  can  be  simplified.  To  do  this  we 
introduce  the  series  resistance  R  and  inductance  L  of  a  unit 
length  of  line  defined  by 

Z  =  R  +  iuL  (4.8) 

and  the  shunt  capacitance  C  per  unit  length  defined  by 


For  the  striplines  considered  here  both  R  and  L  are  in  general 
frequency  dependent.  The  four  limits  to  be  discussed  result 
from  two  different  possibilities  for  Z  in  combination  with  two 
possibilities  for  the  length  of  the  line.  If  the  resistive  part  of 
Z  is  much  larger  than  the  reactive  part  then  L  can  be 
neglected  and  we  describe  the  line  as  being  an  RC  line.  The 
propagation  constant  and  characteristic  impedance  of  the  RC 
line  are 


y  *  (1  +  i)\/ff«C/2. 


Z„  =  <1  -i)VR/2(oC. 


These  forms  are  precisely  those  which  lead  to  the  low- 
frequency  divergences  shown  in  figures  3  and  4  for  normal- 
state  lines.  Because  Re[Z  ,]  «  lm[Z ,]  for  a  superconductor, 
no  similar  RC  behavior  occurs  in  this  case.  When  the 
reactive  part  of  Z  is  much  larger  than  the  resistive  part  the 
line  is  said  to  be  an  LC  line  and  we  have 


y~2VZc  +  ‘ 


z0  =  VT/c  -  i 


Wtc[i+~]. 


<o»R/L.  (4.11) 


2<i>VlC 


such  that  y  is  nearly  pure  imaginary  and  Z  „  is  nearly  real. 
The  length  of  a  line  will  be  called  short  or  long  depending  on 
whether  |  yf\  «  1  or  |  y( |  >>  1. 


For  a  short  RC  line  eq  (4.5)  reduces  to 

* 

VJVs  m  2rT+7r  '  l/RC(t-  <412) 

where  we  have  assumed  R  L  ■  JL/C.  That  is,  in  the  short  RC 
limit  the  transmission  line  can  be  replaced  simply  by  a 
resistor  of  value  €R.  The  complicated  behaviors  of  y  and  Z , 
indicated  in  eq  (4. 10)  thus  combine  to  give  the  expected  low- 
frequency  result.  If  the  A<u  of  a  pulse  satisfies  the  conditions 
of  eq  (4.12),  then  the  pulse  is  received  with  reduced 
amplitude  but  unchanged  shape.  The  reduction  in  amplitude , 
becomes  significant  if  the  line  is  long  enough  that  ( R  is 
comparable  to  RL.  For  copper  at  295  K  the  length  of  line  for 
which  €R  =  R  L  is 

(  _<rds  fp2 


rt. 


<o«R/L  (4.10) 


so  that  the  useable  length  increases  in  proportion  to  5  and  d. 
For  d  —  1  pm  and  e  —  4.  this  length  varies  from  0. 1 1  cm  at 
J  *  0.2  pm  to  11  cm  at  s  =20  pm.  One  concludes  that  at 
least  2  pm  of  dielectric  are  required  to  make  striplines  show 
sufficiently  small  amplitude  reduction.  This  is  the  same  ! 
conclusion  that  was  drawn  from  the  r  =  50  ps  pulses  of  \ 
figure  6  and  for  good  reason.  The  frequencies  below  which 
the  various  normal-state  striplines  can  be  considered  RC  j 
lines  are  listed  in  table  III.  Since  for  r  =  50  ps.  Sir  *  AoV 
2  rr  is  3  x  10*  Hz.  the  295  K  lines  with  s  =  0.2  pm  and  2 
pm  are  on  the  border  between  RC  and  LC  lines.  Moreover, 
we  shall  see  that  the  transfer  function  of  eq  (4.12)  also 
applies  in  the  short  LC  range  so  that,  with  the  exception  of 
those  pulses  displaced  from  the  origin,  all  the  50  ps  pulses 
of  figure  6  show  the  amplitude  reduction  predicted  bv  eq 
(4.12). 

At  4.2  K  one  can  define  a  useable  length  for  copper  lines 
at  low  frequencies  similar  to  that  of  eq  (4.13)  but  with  <r 
replaced  by  or.  This  length  varies  from  7.4  cm  at  s  =  0.2 
pm  to  740  cm  at  j  =20  ^tm.  so  that  the  low-frequencv 
characteristics  of  cooled  copper  lines  should  be  very  good  for 
lengths  typical  of  microcircuits.  As  can  be  seen  from  table 
III,  however,  frequencies  considered  low  are  in  this  case 
considerably  less  than  those  characteristic  of  a  r  =  50  ps 
pulse. 

For  a  long  RC  line  the  frequency  response  is  limited  bv 
the  RC  time  constant  ot  the  line.  The  transfer  function  is 


TaiLI  III.  F  requeney  iHzi  at  which  the  resistive  and  reactive  parts  of  the  senes  impedance  Z  are  equal  for  stnpiines  with  d  m  l  pun.  e  m  1.  and  ranoas 

dielectric  thicknesses. 


Cu.  295  K 
Cu.  4.2  K 


1.5  x  UP 

2.5  x  107 


2.1  x  10* 

3.2  x  [0* 


2.1  x  107 

3.2  x  10* 
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VJVS  m  V375  ex p<  y/RmC/2) 

x  ♦'V'Ec/Sjj  .  <4. 14) 

1/RCt*  «  «  «  Jt/L 

which  includes  both  attorn  lotion  and  phaae  lay.  Impedance 

mismatch  cauaaa  the  -JutUR  attenuation  and  the  ~~  phaae 

4 

lay.  The  propagation  conatant  contributes  the  exponential 
attenuation  and  the  phaae  lay  proportional  to  t.  The  long  RC 
limit  deacribea  the  30  pe  pulae  traveling  the  s  «■  0.2  #im.  t 
■  1cm  line  shown  in  figure  6. 

•  In  the  LC  limit,  impedance  miamatch  ia  present  only  as  a 
second-order  effect  and  can  be  neglected.  The  transfer 
function  in  this  limit  ia 


14.151 

at »  R/L. 

such  that  waves  propagate  with  a  phaae  velocity  of  very 
nearly  l/JES.  An  LC  line  ia  short  provided  t  ia  much  leas 
than  the  wavelength.  (  «  Hat  JCC,  and  in  this  case  eq 
(4.15)  approximates  eq  (4.12).  Thus,  the  short  RC  and  LC 
lines  both  can  be  simply  modeled  with  a  resistor. 

Of  the  pulses  shown  in  figure  6.  all  except  the  50  pe  pulae 
on  the  a  m  0.2  pm  line  ate  described  by  the  transfer  function 
for  an  LC  line,  eq  (4.15).  and  the  impedance  miamaich  can 
be  neglected.  The  propagation  of  these  pulses  can  be 
explained  in  terms  of  the  attenuation  and  dispersion  shown 
in  figure  3.  Because  Av  ia  3  x  10*  and  3  X  lO"  Ha  far  r 
equal  50  and  0.5  pa.  respectively,  the  relevant  forms  for  e® 
and  aM  are  approximately  those  of  region  III.  In  this  region 
atjj  is  proportional  to  Vw*(aee  eq  (3.121)  so  that  in  going 
from  3  X  10*  to  3  x  1011  Ha  the  attenuation  increases  by  a 
factor  of  10.  This  ia  confirmed  in  figure  6  by  the  fact  that  the 
0.5  pa  pulses  show  the  same  attenuation  as  the  50  pa  pulses 
after  traveling  only  V to  the  distance.  Because  atf  is  also 
proportional  to  1/s.  an  increase  in  attenuation  of  a  50  pa 
pulse  on  a  2  pm  line  ia  comparable  to  that  of  a  0.5  pa  pulse 
on  a  20  pm  line.  A  similar  conclusion  can  be  drawn  for  the 
phaae  velocity.  As  eq  (3.9)  indicates,  the  deviation  of 
from  its  high-frequency  value  ia  proportional  to  1  A/ws.  The 
phaae  shift  over  a  given  length  of  line  due  to  this  deviation  is 
thus  proportional  to  -Jaih-  Thus,  broadening  of  a  pulse  and 
.•  displacement  of  its  peak  from  r  ■  i,  change  with  wands  just 
aa  the  attenuation. 

TheLC  approximation  applies  to  all  of  the  pulses  shown  in 
figure  7  for  copper  at  4.2  K.  For  frequencies  typical  of  the 


pulses  shown  tbs  surface  impedance  is  governed  by  the 
axtrenm-anoMioua  equations  with  6.  <  i,  d.  In  due  limit, 
the  anew  — tion  and  phaaa  shift  are  proportional  to  er*ty. 
fnnrliiainna  to  thooo  drawn  for  copper  at  295  K  can 

thus  he  aaade  for  figure  7. 

An  analysis  of  puke  propagation  for  superconducting 
is  shown  in  figure  8  Cor  pukes  with  r  ■  1  and  0.5 
pe.  Far  a  poke  with  r  ■  50  pa  thaw  is  no  sipiificant 
distortion  oven  after  propagating  a  distance  of  10  cm  on  a 
line  with  0.2  fim  of  dielectric.  Poke  degradation  ia  essen¬ 
tially  negligible  Cor  superconducting  linos  until  the  frequen¬ 
cies  involved  begin  to  approach  the  energy-gap  frequency,  at 
which  point  the  stripline  characteristics  dagprade  rapidly.  The 
increase  in  distortion  between  the  1  end  0.5  pe  pukes  of 
figure  8  ia  thus  nearly  aa  peat  aa  that  between  50  and  0.5  ps 
piilsrs  on  a  normal-state  stripline.  The  distortion  evident  in 
figure  8  is  primarily  due  to  dispersion  at  frequencies  some¬ 
what  below  the  energy  gap  as  discussed  elsewhere  in  detail 
[13].  Comparison  of  figures  7  and  8  shows  that  a  supercon¬ 
ducting  line  ia  marginally  superior  to  a  normal-state  line 
even  for  t  *  0.5  p*.  For  longer  pukes  the  superiority  of 
superconducting  lines  is  beyond  question. 


Ficuss  S.  CmpOt  sostspw  mfit wv  6fmr  a  mtkhm  Umm  4JK. 

Tka  Wsa  aero  a  aMM  ka  as  mmmm  i,  -  fit, lOS  Wa>  Sta  ta  m  Saaaaaav 

■a  aataasa  ’at*  >  1.10  a  19.  1.44  «  19.  aaS  1.4*  a  IC  art  ia.  ti.liiwn 
ibMMfO.2.  2.  mi  20  tm.  wynuli. 
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V.  Appendix 

The  numerical  solution  of  eq  (2.7)  and  calculation  of  the 
surface  impedance  for  general  /.  8r,  and  d  follows  a  method 
presented  bv  Mason  and  Gould  [26]  for  the  superconducting 
case.  Applying  the  transformations  of  Mason  and  Gould  to  eq 
(2.7)  yields 

a  t * 

Ex(z)  «  EJJ)  +  ijj  ds'  Ex(s) 

•C,((z  -  d)/l.  (s'  -  *)//)  (A.l) 


Gt(a.  v)  —  ~  +  “  ~  ■De"  +■  («*  —  6u)£i(  |  u  | )] 

3  6 

+  -[»*-  v  -  4)e*r  -  (r*  -  6*)£, («•)].  (A.4) 

6 

The  above  procedure  gives  accurate  results  using  a  small 
number  of  integration  points  except  when  the  depth  of  field 
penetration  is  small  compared  to  d .  In  this  limit,  however,  d 
can  be  taken  as  infinity  and  the  Reuter»Sondheimer  result, 
eq  (2.8),  adopted  in  place  of  the  integral  equation. 


Gt(u.  e)  —  J  dr  J  ds  (t>  -  r) 

where  the  kemal  Ct  can  be  evaluated  in  terms  of  E\.  the 
exponential  integral  [27].  Over  the  domain  of  interest  (u  S 
0,  v  2ul, 


G|(u.  v)  —  —  [(tr*  +  u*  —  10  u  —  6 )r“ 
24 

+  (u4  -  12  u*)£|(  |  u| )] 


+  [(|„s|  _t.s_  10 1  r  |  +  6)e" ' n 

24 


+  (r4  —  12  (  |  v  | )] 

+  ^  (e  —  u)[(u*  +  u  —  4)e“ 

6 

+  (!/*■"  6u)£|l  j  u ) )] 


(A. 2) 


+  3  <  M  +  «■•> 


A  numerical  solution  of  the  above  integral  equation  results 
when  the  integral  is  approximated  by  a  sum  consisting  of  N 
terms  in  which  the  integrand  is  evaluated  at  the  points  * 
s,.  :t,  •••.  :»  By  choosing  s  to  be  each  of  these  s(  in 
succession,  one  obtains  a  set  of  .V  linear  equations  which  can 
be  solved  for  the  unknowns  £,(:<).  In  this  procedure  E  x{d) 
may  be  taken  as  any  non-zero  constant  as  it  merely  sets  the 
scale. 

Once  £,(;)  is  known,  the  integral  over  current  density 
required  for  the  surface  impedance  can  be  obtained  from 


f  ds  Jfis)  -  -  <T  [  ds  £,U)G,((z  -  d)/l.  s/l) 
J  n  ■*  J  0 

-f' dr f' 


(A.3) 


where  for  the  domain  of  interest  lu  sO.r  2  OK 
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A  superconducting  sampler  for  Josephson  logic  circuits 

C.  A.  Hamilton,  F.  L  Lloyd,  R.  L  Peterson,  and  J.  R.  Andrews 

Electromagnetic  Technology  Division.  National  Bureau  of  Standards.  Boulder.  Colorado  90)03 
(Received  9  July  1979;  accepted  for  publication  20  August  1979) 

A  method  is  described  for  automating  a  technique  which  is  used  to  sample  transition  duration 
(rise  time)  in  superconducting  logic  circuits.  The  method  is  based  on  measuring  the  time  at  which 
a  biased  Josephson  junction  switches  under  the  influence  of  an  applied  signal.  The  system 
transition  duration  is  limited  primarily  by  time  jitter  which  is  estimated  to  be  7  ps.  Transition 
durations  of  as  little  as  9  ps  have  been  observed. 

PACS  numbers:  74.50.  +  r 


The  characterization  of  Josephson  logic  circuits  often 
requires  the  measurement  of  very  fast  low-level  signals.  Such 
measurements  are  difficult  to  perform  because  of  the  loss  of 
bandwidth  for  signal  transmission  to  room  temperature. 
Also,  standard  sampling  oscilloscopes  have  transition  dura¬ 
tions  (rise  time)  of  about  20  ps  whereas  high  current  density 
Josephson  devices  switch  in  times  less  than  10  ps.  It  is  there¬ 
fore  desirable  to  provide  for  high-speed  measurements  di¬ 
rectly  on  the  superconducting  logic  chip.  This  may  be  ac¬ 
complished  by  using  a  biased  Josephson  detector  as  an 
amplitude  discriminator  for  the  signal  to  be  measured.' 2 
This  detector  will  switch  to  the  voltage  state  when  the  signal 
plus  the  bias  exceed  the  detector  critical  current.  Since  the 
bias  and  the  critical  current  are  known,  the  signal  level  at  the 
time  the  detector  switches  can  be  determined.  Thus,  by  mea¬ 
suring  the  delay  in  the  detector  output  as  a  function  of  bias, 
rising  edges  on  the  signal  may  be  mapped  out.  This  has  com¬ 
monly  been  done' 1  by  recording  traces  of  the  detector  output 
for  a  number  of  values  of  bias.  The  signal  may  then  be  com¬ 
puted  point  by  point.  This  method  is  somewhat  tedious  and 
may  yield  erroneous  results  if  the  time  base  used  drifts  dur¬ 
ing  the  measurements.  Since  time-base  drift  may  be  as  much 
as  several  ps/min,  it  is  a  serious  problem  for  ultrafast 
measurements. 

Typical  sampling  oscilloscopes  used  to  record  the  de¬ 
tector  waveform  have  an  external  time-base  input  which 
may  be  used  to  specify  the  sampling  time.  By  driving  this 
input  from  the  sampled  voltage,  the  sampling  time  can  be 
locked  to  the  steepest  portion  of  the  detector  output  wave¬ 
form.  The  time-base  input,  therefore,  is  a  measure  of  the 
detector  switching  time.  The  signal  waveform  is  thus  repro¬ 
duced  by  sweeping  the  bias  and  recording  the  bias  versus  the 
time-base  input.  This  method  yields  the  signal  waveform  di¬ 
rectly  and  in  a  time  short  enough  to  minimize  drift  problems. 

Figure  1  shows  in  more  detail  how  this  is  accomplished. 
The  superconducting  integrated  circuit  is  shown  in  the 
dashed  box.  It  consists  of  a  first  junction  J 1  which  serves  to 
produce  a  very  fast  signal  transition  which  is  then  sampled 
by  the  second  detector  junction  J2.  The  sequence  is  as  fol¬ 
lows:  a  tunnel  diode  step  generator  operating  at  60  kHz  ap¬ 
plies  a  repetitive  200-ps  transition  duration  step  to  junction 
J I.  When  itscritical  current  /<.,  is  exceeded,  J 1  switches  to  its 
energy  gap  (2.5  mV)  in  a  time  determined  primarily  by  its 
critical  current  and  capacitance.  This  applies  a  signal  cur¬ 
rent  /,  through  the  8-/2  resistor  to  junction  J  2.  J  2  will  switch 
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when  the  current  /,  plus  the  bias  current  Ib  exceeds  Icl .  The 
output  of  J  2  is  transmitted  to  room  temperature  and  sam¬ 
pled  at  time  t.  The  sampled  voltage  is  compared  with  a  refer¬ 
ence  level,  typically  1  mV,  by  applying  the  reference  as  an 
offset  voltage  to  the  sampling  head  amplifier.  The  difference 
signal  from  the  sampling  head  amplifier  is  integrated  by  in¬ 
verting  amplifier  A 1  and  used  to  drive  the  sampling  time 
base  input.  Thus,  if  the  sampled  voltage  is  below  the  refer¬ 
ence  level,  the  integrator  output  increases,  advancing  the 
sampling  time  to  the  point  where  the  signal  equals  the  refer¬ 
ence.  If  the  sampled  voltage  is  greater  than  the  reference,  the 
sampling  time  is  decreased  in  the  same  way.  This  servo  loop 
thus  causes  the  sampling  time  to  lock  to  the  point  at  which 
the  rising  edge  of  the/  2  transition  crosses  the  reference  level. 
The  time  base  input  is  therefore  a  measure  of  the  time  at 
which  J  2  switches,  i.e.,  the  time  at  which  /,  equals  IcI  -  /„ . 
As  lh  is  slowly  swept  by  the  ramp  generator  (  —  30-s  sweep 
time),  the  X-Y  recorder,  which  plots  -  l„  versus  sampling 
time,  will  trace  out  the  rising  portion  of  the  transition  in  /, . 
Since  the  sampling  rate  (60  kHz)  is  much  greater  than  the 
servo  loop  frequency  response,  noise  and  time  jitter  in  the 
sampling  system  are  averaged.  The  600-ps  delay  line  on  the 
IC  chip  provides  isolation  to  separate  the  /  2  transition  from 
cross  talk  between  input  and  output  lines  at  the  chip  mount. 
The  timing  network  and  analog  switch  pull  the  bias  to  zero 
between  samples  to  reset  the  junctions  to  the  zero  voltage 
state. 

The  advantage  of  this  method  lies  in  the  fact  that  it  can 
measure  transition  durations  which  are  substantially  faster 
than  the  transition  duration  of  either  the  chip  mount  ( —  50 
ps)  or  the  sampling  scope  ( —  20  ps).  The  measurable  transi- 


FIG.  I.  Block  diagram  of  the  transition  duration  measurement  system  uti¬ 
lizing  a  Josephson  junction  as  an  amplitude  discriminator. 
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FIG.  2.  Simulation  of  the  signal  voltage  waveform  (solid  curve)  and  its 
reconstruction  (dashed  curve)  based  on  the  use  of  a  second  junction  as  a 
sampling  device.  The  time  scale  of  the  reconstructed  curve  has  been  arbi¬ 
trarily  shifted  to  bring  the  two  curves  together.  The  simulation  parameters 
are  /, ,  «  2.9  mA.  C,  =  4.8  pF,  /  .  *  1.45  mA,  C.  »  2.4  pF  The  input 
is  a  1-mV/ps  ramp. 


tion  duration  is  limited  primarily  by  the  jitter  in  the  sampling 
system.  The  accuracy  of  the  recorded  waveform  is  also  sub¬ 
ject  to  the  assumption  that  the  detector  junction  waveform  is 
constant  and  that  its  delay  depends  only  on  the  instant  at 
which  its  critical  current  is  exceeded.  These  assumptions  will 
break  down  at  sufficiently  high  speeds.  The  most  significant 
error  results  from  the  fact  that  the  turn-on  delay’  of  7  2  is  a 
function  of  overdrive.  Thus  for  low  bias  values  (signal  near 
maximum),  they  2  transition  has  some  additional  delay.  This 
has  the  effect  of  making  the  recorded  waveform  slightly 
slower  than  the  actual  signal. 

Figure  2  is  the  result  of  a  simulation  of  our  circuit.  The 
solid  curve  is  the  computed  signal  voltage  produced  by  7 1 . 
The  dashed  curve  is  the  computed  result  which  the  above 
described  method  would  produce  in  the  absence  of  jitter.  The 
error  caused  by  variable  tum-on  delay  is  evident  in  the  in¬ 
creased  separation  of  the  two  curves  near  the  signal  maxi¬ 
mum.  The  fine  structure  caused  by  Josephson  oscillations  is 
also  washed  out.  The  critical  current  density  of  both  junc¬ 
tions  is  4200  A/cm1,  giving  them  an  estimated  RyC  time  of 
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FIG.  ).  Photograph  of  the  superconducting  integrated  circuit. 
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FIG.  4.  Recorded  transition  duration  of  the  signal  current  /.  for  a  junction 
with  critical  current  density  of  4200  A/cm’. 

2  ps  (R  v  is  the  junction  normal  state  resistance).*  The  effect 
of  jitter  is  roughly  equivalent  to  passing  the  signal  through  a 
low-nass  filter  of  rise  time  equal  to  the  jitter  width. 

Figure  3  is  a  photograph  of  the  integrated  circuit  used  in 
this  experiment.  It  is  a  seven-level  circuit  which  uses  a  lead 
alloy  process  similar  to  that  described  in  Ref.  S.  The  junction 
areas  are  approximately  40  and  80  /rm:  as  indicated.  The 
microstripline  showing  at  the  left  is  a  small  portion  of  the 
line  used  to  obtain  <M0  ps  of  delay  in  the  output  signal.  The 
parameters  used  in  the  simulation  of  Fig.  2  are  taken  from 
measurements  on  this  circuit. 

Figure  4  is  the  experimental  transition  duration  data 
recorded  for  this  circuit.  The  two  traces,  recorded  about  one 
minute  apart,  are  for  increasing  and  decreasing  bias.  They 
are  different  because  the  sampling  time  drift  decreases  the 
apparent  transition  duration  in  one  case  and  increases  it  in 
the  other.  Assuming  constant  drift,  the  correct  transition 
duration  ( 10%-90%)  will  be  the  average,  or  about  9  ps.  This 
result  represents  the  combination  of  the  signal  transition  du¬ 
ration  estimated  to  be  about  5  ps  and  the  jitter  limited  system 
response  of  about  7  ps. 

A  second  experiment  was  performed  in  which  the 
switching  transition  of  7 1  was  intentionally  rolled  off  by  an 
on-chip  RC  filter.  In  this  case  the  recorded  waveform  was 
the  expected  exponential  rise.  The  time  constant  (3.7  ns) 
matched  the  filter  to  the  accuracy  with  which  the  filter  com¬ 
ponents  were  known  (30%). 

It  is  important  to  note  that  the  very  fast  switching  speed 
reported  here  is  only  one  part  of  the  delay  which  occurs  in 
logic  circuits.  Gheewala*  has  used  a  chain  of  gates  to  mea¬ 
sure  logic  delays  of  as  little  as  13  ps. 

This  work  is  supported  by  the  Office  of  Naval  Research 
under  Contract  No.  N00014-79-F-0020. 
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The  rapid  onset  of  quasiparitcle  tunneling  current  in  superconductor-insulator- 
superconductor  (Josephson)  junctions  at  voltages  above  2A/e  is  being  used  for 
millimeter-wave  heterodyne  mixing.  Junctions  with  a  2-jxm  diameter  and  Rs  =  50  ft 
have  little  capacitive  shunting  at  the  signal  frequency  of  36  GHz.  Because  there  is  no 
series  resistance,  residual  capacitance  can  be  tuned  out.  Double  sideband  conversion 
efficiencies  of  0.32  and  mixer  noise  temperatures  as  low  as  7%,  <  7.  K  =  46  v/k  have 
been  observed.  The  results  are  compared  with  shot-noise-limited  mixer  theory.  Photon- 
assisted  tunneling  effects  are  seen  which  indicate  the  approach  to  photon-noise-limited 
operation. 

PACS  numbers;  07.62.  +•  s,  74.50.  +  r 


The  ideal  nonlinear  element  for  a  classical  heterodyne  In  principle,  a  superconductor-insuiator-superconduc- 

mixer  is  a  switch  which  can  be  driven  between  high-  and  low-  tor  (SIS)  tunnel  junction  approaches  the  ideal  switch  limit 

resistance  states  by  the  local  oscillator  (LO).  At  microwave  more  closely  than  the  super-Schott ky  or  related  devices  be- 

frequencies.  Schottky  diodes  with  /-  V  characteristics  of  the  cause  the  singularities  in  the  density  of  states  on  both  sides  of 

form  I  —  /„(exp(SF)  —  1]  are  often  used  for  this  purpose.  the  junction  cause  an  extremely  sharp  onset  of  qusiparticle 

Conventional  Schottky  diodes  have  5  =  e/kT  =  40  V  at  current  at  the  full  gap  voltage  2J/e.  In  practice,  the S  values 

300  K  and  somewhat  higher  values  at  reduced  temperatures.  thus  far  obtained  from  Pb-In- Au  alloy  junctions  in  the  1 .5- 

The  nonlinearity  in  such  diodes  arises  from  thermally  acti-  4.2  K  temperature  range  are  comparable  to  an  ideal 
vated  conduction  and  disappears  at  temperatures  low  Schottky  at  —  1.3  K.  An  important  feature  of  this  device  is 

enough  that  tunneling  currents  dominate.  The  super-  the  absence  of  series  (spreading)  resistance.  External  micro- 

Schottky  diode1  which  operates  by  quasiparticle  tunneling  wave  circuit  elements  can  thus  be  used  to  resonate  out  the 
between  a  superconductor  and  a  semiconductor  has  junction  capacitance.  Since  there  are  always  practical  limita- 

S=  II  600  V'1  at  T  —  l  K.  Such  a  sharp  comer  on  the  /-  V  '  tionstotheamountofcapacitancewhichcanberesonated.it 
curve  is  advantageous  for  small  signal  recievers  at  short  mil-  is  useful  to  inquire  to  what  extent  capacitance  can  be  made 

limeter  and  submillimeter  wavelengths  because  of  the  diffi-  unimportant  by  adjusting  junction  parameters.  As  the  junc- 

culty  in  obtaining  large  amounts  of  LO  power  at  high  fre-  tion  oxide  thickness  is  reduced.  Rs  decreases  exponentially 

quencies  and  because  LO  noise  problems  are  reduced.  Low-  while  C  increases  linearly.  Consequently,  the  importance  of 
temperature  operation  also  gives  low  mixer  noise.  The  maxi-  junction  capacitance  remains  nearly  constant  if  the  normal 

mum  operating  frequency  of  these  Schottky  devices  is  set  by  state  conductance  per  unit  area,  or,  equivalently,  the  critical 

the  requirement  that  the  junction  capacitance  discharges  Josephson  current  density  7 ,  is  scaled  in  proportion  to  fre- 

through  the  series  (spreading)  resistance  of  the  semiconduc-  quency.  We  have  ojCR  vc=  1  at  36  GHz  for  junctions  with 

tor  twice  each  cycle.  critical  current  density  Jc  =  (O'  A/cm-\  In  order  to  maintain 
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FIG.  I.  Diagram  showing  the  junction  configuration  and  the  arrangement 
of  the  junction  in  the  36-GHz  mixer  block. 


suitable  impedances  for  input  and  output  coupling,  the  junc¬ 
tion  area  must  be  scaled  inversely  with  frequency.  The  need 
for  very  small  junctions  at  high  frequencies  could  be  avoided 
by  using  several  junctions  in  series.  This  would  also  increase 
the  saturation  level  of  the  device.  Since  SIS  junctions  are 
being  actively  developed  for  other  applications,  a  wide  vari¬ 
ety  of  junction  parameters  have  been  obtained.  Junctions 
have  been  reported  with  high  enough  current  density*  ’  to 
give  good  performance  up  to  900  GHz,  beyond  which  mix¬ 
ing  action  will  be  degraded  by  the  breaking  of  superconduct¬ 
ing  pairs. 

The  Pb-In-Au  alloy  junctions  used  in  these  experiments 
were  fabricated  at  NBS,  Boulder,  using  the  photoresist  lift¬ 
off  and  rf  sputter-oxidation  techniques.*  The  junction  geom¬ 
etry  is  shown  in  Fig.  1 .  The  first  electrode  was  evaporated  to 
a  thickness  of  4100  A  on  a  Si  substrate  and  subsequently 
covered  with  a  4500-A-thick  square  of  SiO.  The  SiO  was 
then  removed  in  the  area  of  the  desired  junction  and  the 
junction  oxide  was  grown  by  differential  sputtering  and  oxi¬ 
dation.  Finally,  a  5000-A-thick  counterelectrode  was  added. 
A  0.28  X  1.0  X  19-mm  strip  containing  a  single  junction  and 
its  leads  was  separated  by  cleaving  from  the  40  junctions 
produced  on  each  S-cm-diam.  Si  wafer.  The  junction  areas 
were  varied  from  -2x2  to  - 4 x4^m  to  provide  a  range  of 
impedances. 

The  junction  assembly  was  placed  across  a  full-height 
A'.4-band  microwave  waveguide  in  the  E-field  direction  as 
shown  in  Fig.  I.  The  outer  dimensions  of  the  channel  in 


FIG.  2.  Static  /-Fcurves  are  shown  above  for  a  I.S  K  junction  (a)  without 
and  (b)  with  Put.  Plots  of  IF  amplifier  output  voltage  in  the  frequency  range 
from  30  to  SO  MHz  are  shown  as  a  function  of  junction  bias  voltage.  Curve 
(c)  was  obtained  with  a  30-12  1 .5  K  load  in  place  of  the  mixer;  curve  (d)  with 
a  matched  I.S  K  load  in  front  of  the  mixer:  curve  (e)  with  a  calibrated  36- 
GHz  signal  applied  to  the  mixer  from  a  Klystron  oscillator.  Values  of  mixer 
noise  temperature  were  deduced  from  (c)  and  (d),  and  conversion  efficiency 
from  (e). 


which  the  substrate  was  placed  were  varied  periodically  to 
provide  an  rf  choke.  An  adjustable  stub  3/1/4  in  front  of  the 
junction  and  a  plunger  behind  it  completed  the  low-tempera¬ 
ture  microwave  circuit.  Signal  and  LO  power  from  carefully 
calibrated  36-GHz  Klystrons  were  combined  in  a  10-dB  di¬ 
rectional  coupler  and  introduced  into  the  cryostat  through  a 
section  of  stainless-steel  waveguide.  Three  identical  50-Jf? 
stainless-steel  intermediate-frequency  (IF)  cables  were  in¬ 
stalled  in  the  cryostat.  One  was  terminated  with  a  short  cir¬ 
cuit  to  measure  IF  cable  losses,  one  was  terminated  in  a  50-/2 
cold  load  to  provide  a  noise  source  for  calibrating  the  IF 
amplifier  train,  and  one  was  connected  to  the  mixer.  A  direc¬ 
tional  coupler  was  used  to  inject  a  signal  into  the  mixer  out¬ 
put  to  evaluate  its  coupling  to  the  IF  system.  Measurements 
in  the  IF  frequency  band  of  30-80  MHz  were  made  using  a 
string  of  transistor  amplifiers  with  Tir  =  50  K  followed  by  a 
spectrum  analyzer. 

The  /-  Fcurve  of  an  SIS  junction  at  1 .5  K  is  shown  at  the 
top  of  Fig.  2  with  P LO  s  0  and  also  with  P LO  adjusted  for 
optimum  conversion  efficiency.  The  knee  of  the  /-F  curve 
corresponds  to  S  =  9800  V",  compared  with  7730  for  an 
ideal  super-Schottky  at  the  same  temperature.  The  conver¬ 
sion  efficiency  and  the  noise  in  an  SIS  mixer  are  plotted  on 


TABLE  I.  Propen  in  of  two  SIS  mixers  operated  at  I.S  K.  The  noise  temperatures  and  conversion  efficiencin  are  referred  to  a  single  sideband.  The 
experimental  errors  are  given  as  >0%  confidence  limits. 
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the  same  voltage  scale  at  the  bottom  of  Fig.  2.  ForO<  V  <  1.3 
mV,  the  response  is  dominated  by  hysteretic  Josephson  mix¬ 
ing  which  is  very  noisy.’  In  junctions  with  larger  values  of 
zero  voltage  current,  this  response  (and  noise)  was  as  much 
as  one  order  of  magnitude  larger  than  shown,  but  could  be 
suppressed  with  a  magnetic  field. 

Above  1 .3  mV  mixing  occurs  because  of  the  two  regions 
of  curvature  on  the  quasiparticle  /-  V  curve.  The  noise  in  this 
voltage  range  is  small  and  is  independent  of  magnetic  field. 
A  periodic  modulation  of  the  quasiparticle  mixing  appears 
due  to  photon-assisted  tunneling  because  the  range  of  volt¬ 
age  occupied  by  the  knee  of  the  static  /-F  curve  is  compara¬ 
ble  to  hv/e.  This  is  the  quantum  correction  to  mixing  calcu¬ 
lated  by  Tucker  and  Millea.'  The  properties  of  two  SIS 
mixers  are  shown  in  Table  I.’  The  mixer  performance  has 
been  calculated  from  the  static  /-F  curve  in  each  case  using 
classical  mixer  theory  and  assuming  that  the  capacitance  has 
been  resonated  out  at  the  signal  frequency,  but  that  the  ca¬ 
pacitance  short  circuits  all  harmonics.'  This  is  a  reasonable 
representation  of  our  experiments  because  aiCR  v~  1.  The 
calculated  values  of  conversion  efficiency  are  in  good  agree¬ 
ment  with  the  measurements.  The  calculated  values  of  Fw, 
assuming  only  shot  noise,  are  somewhat  less  than  the  mea¬ 
sured  values. 

The  present  performance  of  the  SIS  mixer  is  sufficiently 
good  to  make  practical  applications  very  attractive.  The  sta¬ 
bility  and  resistance  to  thermal  cycling  of  the  Pb-In- Au  alloy 
junctions  is  a  substantial  advantage  over  the  point-contact 
Josephson  mixer  which  had  higher  noise  (TM  —  54  K)  and 


higher  conversion  efficiency  =  1.3  at  the  same  frequency.' 
The  apparent  ease  with  which  operation  can  be  extended  to 
higher  frequencies  appears  to  be  an  advantage  over  the  GaAs 
super-Schottky  diode.  The  full  benefit  of  this  mixer  will  not 
be  obtained  unless  an  IF  amplifier  with  r,F  —  1  K  can  be 
developed. 

This  work  was  supported  in  pan  by  the  U.S.  Office  of 
Naval  Research. 
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